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SECTION I 


INTRODUCTION AND SUMMARY 


A. INTRODUCTION 

This report presents an eclectic flight control system design 
methodology which is characterized hy highly integrated and synergistic 
use of modern and classical control theory techniques. The classical 
flight control aspects follow conventional practice, e.g., as covered in 
Ref. 10, while the optimal control component of the methodology is simply 
one viable approach among the several that have been suggested (e.g.. 
Refs. 1 -5) . It is the intermix of techniques that promises improved 
effectiveness and more widespread consideration of practical matters in 
the design process. Practicality is viewed herein as satisfaction of 
design qualities which are explicit or implied by traditional flight con- 
trol system specifications and structures. In this context practicality 
indicates compliance with flight control system objectives which have a 
history of success, such as requirements upon closed-loop bandwidth, and 
a desire for low-order controllers. 

The class of problem addressed is in the province of stationary, 
linear, Gaussian, stochastic systems. Modern control theory results for 
linear, stationary stochastic optimal control [i.e., linear, quadratic, 
Gaussian (LQG) theory, e.g., Refs. 6-8] are the primary system synthesis 
tools. Classical control theory results [i.e., transfer functions, fre- 
quency response (Bode) plots, transient response plots, s -plane and Bode 
root locus techniques, and multiloop analysis (e.g.. Refs. 9-1 l)] are the 
primary analytical and interpretative tools. Again, the integrated use 
of these modern and classical techniques is crucial in this design pro- 
cedure. 

One feature typical of flight control systems is that measurements 
provided by many sensors (with the possible exception of accelerometers) 
are practically free from broadband noise. This presents a problem and 
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an opportunity in the modern control framework wherein the controller 
comprises a Kalman filter and regulator combination. The absence of 
broadband measurement noise requires solution of a singular Kalman filter 
problem (treated in Appendix A). The singular filter solutions are of 
lower order than those with full measurement noise. Thus, the ability 
to solve the singular filter problem is a key point in attaining low- 
order controllers, which are desirable for practical solutions. 

The design methodology presented is approached in a case study for- 
mat so that examples of each step are available directly as that step is 
presented. The example applications use a UH-1H helicopter as the object 
of control. Two complete syntheses are available: 

• Longitudinal rate of climb command, airspeed hold 
system for cruise. 

• Longitudinal rate of climb command, groundspeed hold 
system for hover. 

The details of these examples are summarized in appendices, i.e., the 
equations of motion for the UH-1H and its actuation system in Appendix B, 
the hover example in Appendix C, and the 100 kt cruise example in Appen- 
dix D. Appropriate material is drawn from these data bases to illustrate 
the technique and methodological developments in the main body of the 
report. This procedure is presented in a more general fashion so as to 
apply with little change to flight control designs for conventional take- 
off and landing, and vertical/short takeoff and landing aircraft. 

B. PREVIEW OF THE REPORT 

The flight control system design synthesis starts as all syntheses 
must with design requirements. These are treated in Section II both 
generally and specifically. At the outset the analytical synthesis and 
design assessment phases are placed in context with the overall design 
process. The distinguishing features of optimal stochastic control and 
conventional control techniques in accomplishing these phases are then 
drawn. With the differences and distinguishing characteristics between 
optimal and traditional concepts delineated, the next step in an eclectic 
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approach is to establish some connections. This is accomplished using 
a simple idealized control system in which analytical relationships can 
he developed between open-loop function crossover frequency and optimal 
control cost function coefficients. This permits the conventional design 
goal of bandwidth to be related to cost function parameters for an equiva- 
lent optimal system. Design goals are then identified as primary (essen- 
tial for feasibility) and secondary (degree of system quality and via- 
bility relative to alternatives). Control system bandwidth is primary j 
while degree of stability, detailed responses of primary and secondary 
controlled variables, control activity, sensitivity, etc., are examples 
of secondary design goals. In the synthesis methodology presented here 
the primary bandwidth goals are associated with the synthesis procedure 
using LQG optimal controller synthesis techniques. The secondary design 
goals are treated primarily via conventional controls analysis. There 
is of course potential, even essential, interaction required between the 
direct synthesis and analysis procedures to evolve a suitable design com- 
promise, so these roles are not independent. 

Although the theory of LQ£- stochastic opt ima l control is well advanced, 
the necessary engineering art required to permit practical flight control 
solutions as straightforward results from theory application is not well 
established. Section III describes some of the more important problem 
formulation details which we have found can make the difference between 
practical and nonsense solutions. These range from such mundane features 
as actuator command and disturbance representation to the subtle but 
extremely important selection of the control rates as the control features 
to be accounted for in the performance index. 

With design goals established and the problem formulation details in 
hand, the controller synthesis itself is presented in Section IV. This 
is divided into the distinct steps of filter -observer synthesis, regula- 
tor synthesis, and combination to form the controller. At the beginning, 
some important facts about linear stochastic control are reviewed as 
needed to define the synthesis problem and solution. The singular filter 
problems and the implications of requiring a solution which does not re- 
quire a differentiation are considered as an important aspect for flight 
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control system synthesis. These provide the last of the background needed 
in the synthesis steps, so the remainder, which is the hulk of the sec- 
tion, is devoted to specific recommended procedures for the synthesis of 
flight control systems and illustrative examples of their application. 

Some of the syntheses are repeated for different parameter sets to illus- 
trate the sensitivity of the designs to various weightings. 

The design assessment is presented in Section V as an examination and 
expansion of the synthesis results into a composite total picture of sys- 
tem properties, behavior, and design margins. In our optimal synthesis 
method the bandwidths of various loops were central issues, and these 
were ’'automatically" met by direct (or perhaps iterative) exercise of the 
optimal design procedure. The design assessment uses a combination of 
analytical techniques from classical and optimal control theory to expand 
the scope of the results, highlight the dominant properties, and reveal 
the sensitivity of the system's characteristics to uncertainties. In the 
process of the detailed design assessment the optimal synthesis design 
results can be better appreciated in that the relative importance of the 
various terms in the control laws and factors which contribute to parti- 
cular control behavior can be clearly identified. This detailed assess- 
ment and broadened understanding lead directly to simplifications of the 
design which may lead to a simpler and more practical controller and other 
aspects of detailed implementation of the system. Section V begins with 
an outline of the system and controller properties to be examined and a 
systematic checklist-like method for exposing the properties of the system 
and controller. The major content of the section is then devoted to a 
close look at many of the techniques for the design evaluation illustrated 
with the helicopter hovering system example. When all this is done, the 
design is not only assessed but simplified Implementations have also been 
considered, together with their implications on system behavior. 

Various technical details required for completeness or reference are 
contained in the appendices. These include, as mentioned previously, the 
development of the singular Kalman filter solution, aircraft equations 
of motion and numerical data, and data summaries for the example appli- 
cations. 
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C. SOFTWARE 


The extensive computations required to accomplish numerical solutions 
are a necessary adjunct to the methodology presented. Volume II of this 
report contains the user's manual for the software package developed in 
the course of this work. The software integrates the design principles 
from optimal control theory with those from classical control theory, 
permitting the user to design and analyze a control system. The highly 
interactive and modular design philosophy employed permits the user to 
develop a "hands on" appreciation and close connection with each aspect 
of the synthesis and assessment process. In addition, the flexibility of 
the software permits the design to proceed in exactly that sequence of 
steps deemed most appropriate for the particular problems. 

Figure 1 illustrates the block structure which forms the backbone of 
the software package. Each block is a separate executable file which 
performs the specific task noted in the figure. The underlying file sys- 
tem structure allows the various blocks to communicate information to one 
another. At each block, a "problem file" is read and/or written. These 
problem files store all intermediate results, which the user may access 
via the service routine. The service routine selectively reads and for- 
mats to a line printer any user -re quested elements in a given problem 
file. In this way, the user can examine not only the final synthesis 
results but also those intermediate results of interest which might pro- 
vide increased understanding of the results obtained. 

Two types of blocks comprise Fig. 1. The first type (Blocks 1, 
h, 6, 7> 10) implements various aspects of the optimal control design 
process using software adapted from Slater's version of Bryson and Hall's 
OPTSYS (Ref. 29). The output from these blocks is stored on problem files. 
The second type of block (Blocks 2, 8, 9) provides links between the opti- 
mal control design and the classical control analysis techniques. The 
classical control software is organized in a similar fashion, as shown 
in Fig. 2. 

A typical sequence of steps in designing an optimal, control system 
is to begin with Block 1 at the top of Fig. 1 and work clockwise through 
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to Classical 



Figure 1 . Software Structure, Optimal Control 
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Figure 2. Software Structure, Classical Control 








all blocks, making the appropriate indicated excursions to the blocks in 
Fig. 2 to use the classical control analysis techniques. Blocks may be 
repeated to cycle a portion of the design, while some blocks ( such as 
the sensitivity block, for example) may be bypassed altogether. 

A complete description of the use of this package, along with a simple, 
but illustrative, example is given in Vol. II (Ref. 25). 


SECTION II 


FLIGHT CONTROL SYSTEM DESIGN REQUIREMENTS 


A. INTRODUCTION 

The flight control system’s (FCS) basic overall function is to provide 

• Stability. 

• Desired responses to specified inputs. 

• Suppression of the effects of disturbances, component 
variations, and uncertainties. 

• Modification or elimination of vehicle cross- coupling 
effects. 

Each specific function can be achieved by appropriate application of con- 
trol techniques in a manner consistent with overall system requirements. 

It is, of course, necessary to have a keen understanding of what is 
"appropriate" and "consistent." This is best arrived at via the overview 
of the design process given in the next three paragraphs (from Ref. 10). 

The aim of the FCS design is to produce a functional 
system that performs its assigned tasks "satisfactorily." 

The design process leading to this end can be broken into 
several phases that are more or less chronological, yet 
are extensively interrelated and interconnected. A typical 
set of such phrases might include the following. 

1) Establishment of System Purpose and Overall System 
Requirements. At the design stage, system purpose can 
be equated with the mission phase or task definitions. 
Requirements are partially derivable from the functions 
needed to be performed to accomplish these mission 
phases (operational requirements), and less directly 
from the characteristics of the interconnected compo- 
nents and the environment in which they operate (implied 
re quir ement s ) . 

2) Determination of Unalterable Element, Command, and Dis- 
turbance Environment Characteristics. Typically, the 
characteristics of some component parts of the system 
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are not easily changed by the system designer* In aero- 
nautical control and guidance such relatively "unalter- 
able” elements often include the vehicle itself and 
possibly the control surface actuators and some of the 
motion quantity sensors. The "structure" of the com- 
mands and disturbances is not subject to the choice of 
the designer but is instead a direct consequence of the 
mission or task and the environment. The latter may, 
however, be open to some interpretation. 

3) Evolution of Competing Feasible Systems. Usually, 
requirements can be met in more than one way, e.g., 
with different sensed motion quantities and equaliza- 
tion elements which are completely alterable within the 
limits of physical realizability and practicality. Then 
it is possible to evolve competing systems that become 
candidates for selection on the basis of their various 
desirable properties. 

4 ) Competing System Assessment, System Selection. The com- 

peting systems can be compared on a very large number of 
bases which can be divided into two categories: design 

quantities and design qualities . Design quantities 
include the dynamic performance (relative stability, 
accuracy, speed of response or bandwidth, etc.) and the 
physical characteristics (weight, volume, power or energy 
consumption, etc.). Design qualities, the so-called 
"-ilities, " include safety, operational capability, 
reliability, maintainability, cost, etc. An optimum 
system is one that has some "best" combination of all 
these features. 

5) Detailed Study of the Selected System. Once a best system 
has been selected, it is still necessary to validate it 
for all nominal and abnormal operating conditions. The 
components that do not yet exist as hardware must be 
specified, designed, fabricated and tested as components. 

As many of these as possible should be assembled in a 
series of system simulations which culminate in flight 
tests of the complete system in its actual operating 
environment. At each stage of the testing process the 
assumptions that were made in previous phases of the 
design should be checked for validity. If actual con- 
ditions violate the assumptions, a new iteration of the 
design should be begun at the point at which the incor- 
rect assumption was made. 

The above steps in the orderly process of evolving, or synthe- 
sizing, a system that satisfactorily meets all its objectives are 

governed initially by a set of functional requirements stemming 

(see Item 1 above) from operational needs and system integration 
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implications. These latter, implied, requirements depend on 
the final system selected (see Item which is then completely 
specified functionally hy the total set of (functional) require- 
ments. The converting of these requirements on each element of 
the system into usable hardware generally requires a specifica- 
tion that includes not only the functional requirements but also 
applicable specifications relating to hardware design and fabri- 
cation practices, including extreme environmental factors. 

In any aeronautical system design the requirements for sub- 
systems evolve in a pyramidal fashion, and become more numerous 
and detailed as definition of the actual equipment is approached. 
The apex of the pyramid is the mission purpose and definition 
(see Fig. 3). Immediately below this central point are three 
blocks involving considerations that interact strongly in the 
earliest preliminary design stages: mission phases, vehicle 

operating point profile, and guidance possibilities. When the 
mission is realistic, one or more feasible vehicle operating- 
point profiles are joined with one or more guidance possibili- 
ties to enable the overall system to perform through the con- 
stituent phases of the mission. 

It is within this larger context that practical flight control system 
design must be considered. This project, while restricted in scope to the 
analytical/quantitative aspects of design, must remain responsive to, and 
compatible with, this context. Within this scope, it is clear that we are 
mainly concerned with the quantitative aspects of the 

Determination of Unalterable Element Command and 
Disturbance Environment Characteristics 

Evolution of Competing Feasible Systems 

Competing System Assessment, System Selection 

phases of the overall design process. 

Address of these design phases using a non- traditional flight control 
system synthesis technique, such as linear quadratic, Gaussian (LQG) opti- 
mal control, requires that we have the distinguishing features with respect 
to the more traditional techniques firmly in mind. These distinguishing 
features are summarized in Table 1 . As shown by the table the differences 
are confined largely to the Evolution of Competing Systems design phase. 
Conventional control theory considers the entire entourage of inputs and 
disturbances ab. initio, whereas the optimal formulation leaves direct con- 
sideration of transients until a later phase of design. Control point 
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Figure 5. Evolution of Flight Control Requirements 
(from Ref. 10) 
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TABLE 1 


DISTINGUISHING FEATURES OF CONTROL SYSTEM 
ANALYSIS/SYNTHESIS TECHNIQUES 


DESIGN PROCESS PHASES 


ANALYSIS/SYNTHESIS TECHNIQUES 


LQG OPTIMAL STOCHASTIC CONTROL 


CONVENTIONAL CONTROL THEORY 


Determination of Unalter- 
able Element, Command and 
Disturbance Environment 
Characteristics 


Similar linearized, constant coefficient equations of motion and Gaussian 
statistics for random commands/disturbances are needed 

Entire ensemble of discrete command 
input forms 


Measurements : Similar alternatives based on traditional/feasible sensors 


are applicable. 

Control Point Utilization: Deter - 


Evolution of Competing 
Feasible Systems 


mined by solution based upon cost 
function parameters. Alternative 
solutions obtained by varying cost 
function parameters. 


Controller Equalization: Deter- 


mined by solution of optimal filter- 
observer problem. Alternative solu- 
tion obtained by altering meas- 
urement noise and/or disturbance 
environment assumptions. 

Realization; Designer-conducted re- 


configurations of standard regulator 
plus filter-observer realizations. 


Control Point Utilization ; Deter- 
mined by considerations of; typical 
operational utilization (e.g., by 
pilots); design-conducted transfer 
function survey based on multivari- 
able control theory; separation of 
trim and high-frequency control ac- 
tions; controls coordination and/or 
decoupling. 

Controller Equalization ; Determined 
using combinations of; s-plane and 
Bode root loci and conventional fre- 
quency response synthesis procedures 
with multivariable control theory, 
response/error coefficients, and 
indicial (sometimes ramp) responses. 

Realization; Inherent in synthesis 


procedure . 


Transient Response; Same transient response characteristics required. 


Competing System Assess- 
ment, System Selection 


RMS Response ; Same rms response characteristics required. 

System Sensitivity ; Same insensitivity characteristics required. 
RMS State Estimation Error: Addi- 


tional basis for selecting sensor 
complement . 

Economy of Controller Equalization; 


Direct result of sensor complement 
selection and measurement noise 
assumptions. 


Economy of Controller Equalization ; 
Inherent in synthesis procedure. 
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utilization factors at first glance appear to "be quite different — yet 
the considerations used in the conventional approach must ultimately he 
reflected into the optimal if an adequate design is to he achieved. LOG 
optimal control techniques provide an additional basis for selecting a 
sensor complement (via rms state estimation error evaluation) and permit 
an automated controller solution (albeit ordinarily in a form requiring 
considerable manipulation before the flight control system implementation). 
The automated solution feature is achieved at the expense of having to 
specify sets of cost function parameters. We shall show later that this 
can be analogous to selection of control loop bandwidths when conventional 
frequency response synthesis techniques are used. Additional features of 
the LQG optimal control technique are that closed- loop system stability is 
inherent in every controller solution (but acceptable insensitivity to 
component variations and uncertainties is by no means assured), and that 
the closed-loop system can be designed specifically for a given command 
and disturbance input environment. 

The assessments needed at the completion of a trial synthesis to 
evaluate the degree of system adequacy are substantially identical and 
any shortcomings must be corrected by additional iterations. Thus, in the 
larger context, the impact of using one or another control system analysis/ 
synthesis technique should be nil. The two techniques are simply system- 
atic methods for evolving and assessing competing system designs. On the 
other hand, the impact will not be nil if the two techniques are not 
equally effective in evolving practical candidate competing systems. This 
in turn emphasizes the importance of technique in applying LQG optimal 
stochastic control to the flight control problem when one chooses to take 
advantage of this more highly automatic design procedure. 

These aspects of technique fall into four areas: 

• Establishing design goals consistent with operational 
and implied requirements. 

• Formulating the problem (equations of motion and input 
and disturbance environment models) in terms appro- 
priate to control requirements determined by mission 
phase tasks. 


13 


• Selecting appropriate sensor complements. 

• Selecting appropriate cost function parameters. 

These matters of technique are discussed in this and the following section. 

B. CONNECTIONS BETWEEN OPTIMAL AND TRADITIONAL 
CONTROL CONCEPTS FOR A SIMPLE IDEAL SYSTEM 

The last section explored the differences and distinguishing charac- 
teristics between optimal and more traditional concepts^ the present one 
will emphasize their connecting relationships. This will be accomplished 
in the guise of a very simple idealized control system. Simplicity permits 
the development of analytical relationships among bandwidth, open-loop func- 
tion crossover frequency, transient and stochastic system response charac- 
teristics, sensitivity, and optimal control cost function coefficients. 

The idealized system is chosen not only for its analytical simplicity but 
also for its practical significance in the application of frequency response 
synthesis techniques. The system, in fact, serves as a prototype or ideal 
solution which the designer attempts to emulate in far more complex applica- 
tions . 

1 . Relationships Between System Response 
Characteristics and Sensitivity to 
Crossover Frequency and Bandwidth 

A block diagram for a single-loop feedback control system is given in 
Fig. 4 . Limiting forms of the closed-loop system transfer functions in 
terms of the magnitude of the open-loop function, |g(s)|, are given in 
Table 2. The key point of these limiting forms is that three regions 
exist. In the first, when |g[ » 1 feedback control is fully active, and 
the system output follows the system input closely with little error. The 
second regime is at the other extreme, j G | « 1, where feedback control 
is essentially not operating and the system error and input are nearly 
the same. The third frequency regime, characterized by |g| of the order 
of 1 , is the region where feedback relationships are most complex in that 
the controller and controlled element dynamics interact to create the 
dominant closed- loop system modes. In a typical closed-loop system which 



G(s) 



Open-Loop Transfer Function 


Cf(s) 


Xfs) 

E(s) 


Closed-Loop Transfer Functions 


Response /Command Relationship 


X(s) = G(s) 
W[sJ 1 + G(s) 


Error / Command Relationship 

E(s) = 1 

W(s) 1 + G(s) 


Figure Prototype Closed-Loop Control System 


TABLE 2. LIMITING FORMS FOR TRANSFER FUNCTIONS 


TRANSFER 

FUNCTION 

|g| 

» 1 

0(1) 

« 1 

X/W 

1 

G 

1 + G 

G 

E/W 

l/G 

1 

1 + G 

1 
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is low pass in character, these three regimes correspond to low frequency, 
very high frequency, and crossover frequency regions, respectively. 

A key point to appreciate in Table 2 is that |g(s) | = 1 provides a 
clear demarcation between the limiting form approximations. Further, 
while the limiting forms represent idealizations of the response charac- 
teristics desired for the controlled variable, x, and for the error vari- 
able, e, these particular idealizations are essentially independent of the 
details of G(s) in the regimes where the limiting forms apply. 

On the other hand, the form of G(s) in the crossover region about 
| G( s ) | = 1, which constitutes a transition between the limiting cases, is 
of great importance. It is here that the battle between high performance 
and stability is joined. The ideal form for G(s) in this crossover region 
approximates G(s) = 0 ) c /s, where the gain, o^., coincides with the crossover 
frequency (i.e., JoDc/sl . = l). A basic design principle in classical 

frequency response design is to adjust G(s) so that; 


G(s) 

= o^/s 

for 

s = jat 


|g( ja>) | 

» 1 

for 

0) « m c 

(D 

|G(jco)| 

« 1 

for 

CD » GD C 



The closed-loop consequences of this type of adjustment include good 
response following of commands and good error reduction over an input fre- 
quency range less than cj^., with adequate closed-loop system stability. 
These characteristics will be obtained regardless of the detailed nature 
of G(s) in the limiting form regimes. This prescription can be general- 
ized in a straightforward way for multivariable control applications 
(Refs. 10 and 11). It can be modified to take more precise account of 

/ — — crO 

phase margins significantly less than it/2 rad by adding an e term to 
the approximation for G(s) in the crossover band (Refs. 10 and 12). The 
optimal controllers with which we shall be concerned here tend inherently 
to be highly stable, so our simple form G(s) = o^/s remains an appropriate 
ideal. 
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The dynamics, open- and closed-loop frequency response, and indicia! 
response characteristics of the idealized system are summarized in Fig. 5. 

The closed- loop frequency response functions in Fig. 5 show that 
crossover frequency coincides with the classical definition of signal 
handwidth. That is, bandwidth is associated with the —5*01 dB point of 
the frequency response. Notice that for this particular choice for G(s) 
the error suppression bandwidth and controlled variable response band- 
width are both equal to If 2$G(joi) «/2 when |g(joo) | = 1, these 

bandwidths are unequal, but both are still approximated by when ade- 
quate phase (stability) margin is provided for in G(s). [That is, when 
|g( joj) | = 1, (0 - cpj^) > 4G(«jo>) > (cpM - it/ 2 ) 111 order to insure a phase 
margin of at least cp^ radians. ] 

The transient response characteristics for this prototype system are 
illustrated by the indicial responses in Fig. 5 . These are: 


1 — exp(— co c t) 

(2) 

exp(— (u c t) 

(3) 


The responses to a unit initial condition on the controlled variable (con- 
sidered to be the state) are: 

x = exp(— oo c t) = — e (4) 

The transient responses are smooth, free from overshoot, and are charac- 
terized by a response time, 3/(^ c . The steady-state error response to the 
step input is zero. 

Consider next the root -me an- square (rms) response characteristics for 
this prototype system. The block diagram for the system is in Fig. 6. 

The stochastic input is generated by low-pass filtering unit white noise. 
Normalized rms controlled variable (o^/o^-) and error (crg/o^) responses 
are plotted in Fig. 7 as a function of the normalized stochastic input 
half-power frequency (o^/oi c ). For stochastic input half -power frequencies 
which are less than the crossover frequency. Fig. 7 indicates that 
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Amplitude Ratio (dB) 



X(s) _ G _ tJc 
W(s) l + G s + oj c 


E(s) _ I _ s 
W(s) l + G s + ou c 



Bode Root Locus 


-cue 

- 0 - 



s- plane Root Locus 


a) System Transfer Function 


b) System Dynamics .* Root Loci 



log cu — ► 


c) Open and Closed- Loop Frequency 
Response 



d) Closed- Loop System 
Indicia/ Responses 


Figure 5. Characteristics of Idealized System 
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attenuation of the controlled variable rms response with respect to the 
rms input does not exceed 1 1^/2 , while attenuation of the rms error 
response with respect to the rms input always exceeds the lesser of 1 j^/~2 
and . This demonstrates that stochastic system performance for the 

prototype is consistent with desired design properties. 

Because the open-loop system phase lag is confined to — 90 deg, this 
ideal system is stable for all values of gain. Its response characteris- 
tics will, however, be affected by changes in system features. These are 
conveniently treated using first-order sensitivity factors which relate 
shifts in open-loop parameters to their consequences as closed-loop pole 
shifts. In general, the variation d}^ in a closed- loop pole (Xi), due to 
changes in open-loop gain (dK), open-loop zeros (dzi) and open-loop poles 
(dpi) is given by (Ref. 1 0) : 


. , n ntt-n 

= sj — + £ sLdz i+ £ sj^dPi (5) 

d=i J j=i 

The S 1 quantities are the first-order sensitivity factors. The subscript 
and superscript notation indicates that a differential increment in the 
open-loop parameter (defined by the subscript) results in a differential 
increment of the ith closed-loop root (denoted in the superscript) which 
is equal to the sensitivity factor times the open-loop parametric varia- 
tion. Provision is made for an excess of m poles over zeros in the open- 
loop transfer function. The gain sensitivity is given by, 


[dG/dsJ 


s=X- 


( 6 ) 


which for this simple ideal system is — gd c . This leads to the obvious 
result; 


dXi 


— 01 - 


d K 
K 


( 7 ) 
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The crossover frequency is thus the gain sensitivity factor which relates 
proportional changes in open-loop gain to closed- loop root shifts. 

At this point we have established connections among open- loop func- 
tion crossover frequency, closed- loop system transient responses, stochas- 
tic input responses, and performance and system sensitivity to parameter 
variations. In all cases the bandwidth or crossover frequency is the sig- 
nificant (in this case solitary) design variable. These results can be 
summarized by the design aphorism T, with stability assured, everything else 
is related to bandwidth. " This is of course the complete story for the 
case at hand. 

2. Connections Between Desired Open-Loop Function 
Crossover Frequency and Optimal Control System 
Cost Function Coefficients 

Given the goal that the resulting system exhibit the favorable proper- 
ties enumerated above, let us now turn to the task of designing a controller 
C(s), using linear-quadratic-Gaussian optimal control procedures. For the 
system to be optimal requires that a cost function, such as 

J = / [lR e2 + *rU 2 ] dt ( 8 ) 

•'O 

be minimized. Assume that the controlled element component of the system 
Y c (s) has a form in the frequency region near cq, which can be approximated 
by: 


Y c 



n > 0 


(9) 


The desired controller properties are to be such that the total open-loop 
function G(s) = C(s)Y c (s) exhibits the three properties enumerated above, 
i.e., |g(s) | » 1, |g(s) | = |cu c /s|, and |g(s) « 1 at frequencies below, 

near, and above the crossover frequency, o^., respectively. The optimal 
control procedure can accomplish this only if the weighting factors in the 
performance index or cost function can be selected to reflect these desires. 
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Fortunately, this problem has been solved, and it has been shown in Ref. 13 
that if the cost function coefficients, q R and r R , are chosen to satisfy 


^R 



( 10 ) 


then the open- loop function G(s) will tend to have the three desired prop 
erties. Figure 8 presents the approximate asymptotes for the G(s) ampli- 
tude ratio frequency response resulting when q R and r R satisfy Eq. 10. 
Notice that crossover frequency, 0 i c , for the low-frequency asymptote is 
maintained regardless of the controlled element characteristics in the 
frequency region near ao c . 

The actual frequency response can be expected to be close to its 
asymptotes, and controlled variable and error response bandwidths are 



Figure 8. Relationships Between Crossover Frequency 
and Cost Function Coefficients 
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c|0 



approximately equal to the open-loop function crossover frequency. Thus, 
we have here the basis for specifying q^ and r^ in a cost function that 
induces an optimal control system which satisfies a given bandwidth 
requirement. This observation is the cornerstone of the method presented 
in later sections for design of practical, optimal flight control systems. 

C. DESIGN GOALS 

Design goals are deduced from mission tasks by considering what is 
required to perform the task (e.g., the landing approach task requires 
speed and path control). Design goals can be identified as primary (task 
cannot be accomplished satisfactorily if not satisfied) or secondary (degree 
to which minimum standards of task performance are exceeded) . That primary 
design criteria be satisfied is essential for feasibility, whereas secon- 
dary design criteria satisfaction is more a matter of system quality and 
viability relative to the alternatives. The design goals elaborated below 
are not unique or exhaustive. They merely represent one recommendation 
which has been tested successfully. 

1 . Elementary Mission Task Requirements 

Experience and mission task analysis (e.g.. Ref. 10) have demonstrated 
that certain modes of flight control system operation are required on vir- 
tually all rotorcraft missions. In our demonstration here of practical 
optimal rotorcraft flight control system design we shall focus on these 
more or less standard modes of operation and disregard special modes pecu- 
liar to less universal mission tasks. (The same approach and philosophy 
for determination of requirements and for design development carry over 
for these special modes, of course. ) 

Table 3 lists standard operational modes of interest for rotorcraft. 
These are grouped mainly by control axis, although no distinction between 
the pitch and vertical axes is made because of the significant differences 
in appropriate control technique at hover and at cruise. Entries in 
Table 3 are ordered to roughly reflect increasing level of flight control 
system responsibility (in distinction to pilot responsibility) in the over- 
all conduct of these mission tasks. 
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TABLE 3 


STANDARD FLIGHT CONTROL SYSTEM 
MODES FOR ROTORCRAFT 


Longitudinal (pitch and/or vertical control axes) 

Damping augmentation 
Rate eommand/attitude hold 
Attitude command 
Airspeed hold 
Translational rate command 
Rate-of-climb command 
Altitude hold/path track 

Hover (longitudinal and vertical position hold) 

Lateral (roll control axis) 

Damping augmentation 
Rate eommand/attitude hold 
Attitude command 

Heading command /ho Id (at mid and high speeds) 
Translational rate command 
Path track 

Directional (yaw control axis) 

Damping augmentation 

Rate command/attitude hold 

Heading command/hold (at low speeds) 



2. Primary Design Criteria 


Primary design criteria involve two considerations: 

• Identification of outer-loop error variables per- 
tinent to a mission task (e.g., for "follow com- 
manded rate of climb, " rate of climb error is the 
pertinent outer-loop error variable to control to 
zero) . 

• Specification of the bandwidths over which outer- 
loop errors and other key variables (e.g., pitch 
attitude) are to be controlled. 

The first item is fundamental to accomplishing the mission task. The 
second is, in effect, a catch-all specification of the performance required. 
This was demonstrated previously in presenting the relationships among 
bandwidth, open-loop crossover frequency, closed-loop system dynamic steady- 
state, and stochastic performance. Stability per se is not identified here 
as a primary design criterion because the design will be evolved using an 
optimal control formulation which assures stability. Thus we return again 
to the aphorism "with stability assured, all else relates to bandwidth." 

The other key variables in the second item require further explana- 
tion. If bandwidths are specified for outer-loop error control only, then 
the bandwidths (or crossover frequencies) provided by the optimal control 
design algorithm for all other variables will be at the minimum levels 
required to achieve the outer-loop error bandwidth objective. Modes par- 
ticipating significantly in the outer-loop error response may be separated 
from the outer loop * s dominant mode frequency by as little as a factor of 
two, while non -participating modes may be changed hardly at all. Fre- 
quently a factor of two separation is not adequate. For example, it is 
usually inadequate for path error and attitude modal separation in view 
of customary flying qualities requirements related to motion harmony and 
the partitioning of attitude and path control frequency bands. It can also 
be inappropriate when the effective actuator bandwidth is fixed because of 
an unalterable hardware selection or when significantly higher actuator 
bandwidth is required for other reasons. In either event, modal separa- 
tion may be at the designer's choice merely by specifying target bandwidths 
for control of these other key variables. 
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Selection of outer-loop error variables is subject to fundamental 
limitations. The number of outer -loop error variables which may be con- 
trolled independently cannot exceed the number of independent controlled 
element control points. Further, these outer-loop error variables must 
be controllable via the independent control points. Selection of the 
outer-loop error variables must be consistent with these two fundamental 
limitations. 

Outer-loop error and other key variables for helicopter flight control 
are summarized in Table 4. Entries in this table correspond to the Stan- 
dard Flight Control System Modes for Rotorcraft listed in Table 3- The 
outer-loop error variable for each mode is the right-most entry in each 
string not appearing under an integral. Variables to the left of the 
outer-loop error in each string are the "other key variables." In several 
instances the integral of the outer-loop error appears at the extreme right 
of the string. These integrals are included to assure zero steady-state 
error to a step command input. Integral terms are required only when a 
command input is explicit in the outer-loop error. 

Tables 5 and 6 recommend bandwidth requirements for the outer-loop and 
other key variables. In several Instances alternative recommendations are 
given for precision and other (non -precision) tasks. The recommendations 
generally are in the form of a lower bound on bandwidth. These lower bound 
values have been found to be reasonably adequate (albeit lower than optimum) 
in practice. Larger bandwidths generally correspond to increased perfor- 
mance (faster response, smaller rms error). However, larger bandwidth 
systems, because they are more sensitive to parasitic nonlinearities and 
to the influence of higher frequency modes, may produce "twitchy" system 
response characteristics to which pilots object. Furthermore, increases 
in bandwidth tend to equate with increased control authority and rate 
requirements. Therefore, extreme caution should be exercised when requir- 
ing more than twice the bandwidth lower bounds listed. 

Tables 5 and 6 do not include bandwidth recommendations for the inte- 
grals of outer-loop error variables. These low frequency Integral modes 
cause "tails" in the transient response. We have found that the bandwidth 
for the integral should be 0.2 to 0.^ of the bandwidth for the outer-loop 
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TABLE k 


OUTER-LOOP ERROR AND OTHER KEY VARIABLES a 
Pitch Control Axis 


Actuator : 

& B> 

S B 



Damping augmentation: 


5 B> 

<1 


Rate command/attitude 
hold; 


b B> 


/ qe dt 

Attitude command; 


^B> 


9 e , / 0 e dt 

Airspeed hold; 

8 B> 

§B> 


8* u AS e > S u AS e dt 

Translational rate 
command: 

8 B> 

b B’ 

A 

0, x e , / x e dt 

Hover (longitudinal 
position hold) : 

^B> 

5 B' 


0, i* 3 , x e , / x e dt 

Vertical Control Axis 





Actuator ; 

8 C; 

8 C 



Damping augmentation; 

&C, 

8 C> 

h 


Rate-of-climb command; 

8 C> 

8 c> 


/ h e dt 

Altitude hold/path track; 

8 Cj 

8 C> 


h e , / h e dt 


a See Appendix B for variable definitions. 

^Used only when adjustments to the closed-loop damping are 
required. 

(continued on following page) 
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TABLE 4 (CONCLUDED) 


Lateral Control Axis 
Actuator: 

Damping augmentation; 

Rate coimnand/attitude 
hold; 

At t i t ude c ommand ; 

Heading command; 

Translational rate 
command; 

Path track/hover (lateral 
position command) ; 

Directional Control Axis 
Actuator; 

Damping augmentation; 

Rate command/attitude 
hold; 

Heading command; 

b Used only when adjustments 
required. 


b A’ b A 
&A, 5 A , P 

5a, 5 a , p e , / P e dt 

5 a , 5 Aj P b , cp e , /cp e dt 

5 a , S a , P fe , q>, ijr e , f* e dt 

5 A , & A , P b , cp, Y e , f Y e dt 

5 A , 5 a , p b , cp, Y b , Y e , / Y e dt 

5r, 5r 
5r, 5 r , r 

5r, 5 r , r e , f r e dt 
Sr, 5r, r b , f e , / t e dt 

the closed-loop damping are 
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a These are minimum values. Actuator bandwidths are usually considerably larger in typical instal- 
lations. Actuator bandwidths may be specified rather than being design parameters. 



error variable. A good compromise or starting value is 0.3* Smaller 
values tend to produce ’’tails" in transient responses which have smaller 
magnitude but longer duration, and vice versa for larger values. 

Next, consider the secondary design criteria. 

3. Secondary Design Criteria 

The secondary design requirements involve a great many considerations. 
These may be classified into two categories: verification that primary 

design criteria are satisfied, and investigation of system characteristics 
determining the quality of performance. Secondary criteria in both cate- 
gories are applied after the fact (or at least after the first iteration) 
as far as the optimal control synthesis is concerned. Therefore, they may 
be thought of as design assessment or evaluation criteria. 

Verification that primary design criteria are satisfied is necessary 
because the relationships between bandwidth, crossover frequency and the 
cost function coefficients given previously are approximate. While we 
have found these approximate relationships to be robust, verification that 
desired bandwidths and crossover frequencies are actually obtained is recom- 
mended. 

The second category includes the numerous items which have been found 
effective in exposing the strengths and weaknesses of control system 
designs. A complete list of system design features and properties which 
should be considered in this context is given in Table 44. The system 
features listed there, and exemplified in Section V, include ; 

• System stability characteristics and margins 

• Responses 

Primary controlled variables 

Secondary controlled variables 

Control activity 

• Sensitivity 

Key modes 

Parasitic nonlinearities 

• Gain levels 

• Sensor /equalization economy 
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The investigation of system characteristics to determine the quality of 
performance and to quantify these system features pr imar ily involve, hut 
are not confined to, such items as: 

• Closed-loop frequency response 

• Steady- state error response 

• Transient response 

• RMS response 

• Actuator displacement and rate limits and activity 

• System insensitivity 

Closed-loop frequency responses are used to determine that physical 
responses of the system are free from resonant peaks, that the response 
hand is flat, etc. Steady-state error characteristics can he determined 
hy examining the very low frequency closed-loop transfer function for the 
outer- loop error response to the command input as s — 0. Transient 
responses can he inspected for overshoot, slow "tails," smoothness, etc. 

RMS responses of the errors, controls and other key variables can be com- 
pared against specifications and error budgets. In particular, three times 
the rms actuator displacement (ct x ) and rms actuator rate ( 03 -) can be com- 
pared with the available displacement and rate capability to insure that 
the fraction of time these capabilities are expected to be exceeded is less 
than 0.26 percent. Control activity in terms of positive -going zero cros- 
sings of the trim setting can be calculated using c x / 2 jta x . 

System sensitivity characteristics can be examined from several points 
of view. The optimal control solution itself facilitates computation of 
the closed-loop regulator pole first-order sensitivity to plant coeffi- 
cient changes and to regulator gain changes. The dual relationships 
between the regulator and filter solutions can be used to evaluate the 
closed- loop filter pole first-order sensitivity to plant model coefficient 
changes and to filter gain changes. First-order sensitivity effects are 
also available via the partial fraction expansions of the closed-loop trans- 
fer functions used in computing transient responses. The modal response 
coefficients of the partial fraction expansions are gain sensitivities 
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(Ref. 10). Finally, sensitivity to use of the designed controller at 
flight conditions other than the design case, to use with more comprehen- 
sive models of plant dynamics, and to use of simplified approximate ver- 
sions of the designed controller can he evaluated. 

D. APPLICATION OF THE DESIGN 

CRITERIA IN SYNTHESIS 

Primary and secondary classes of the design checklist items for rotor- 
craft flight control systems have heen presented above. In the synthesis 
routine the primary requirements are to he addressed directly via the LQG 
optimal controller synthesis procedure. This is accomplished first hy 
formulating the problem so that appropriate outer -loop error quantities 
are included, and second by making use of fundamental relationships between 
response bandwidth and cost function coefficient values. 

Design evaluation for satisfaction of the secondary criteria is after 
the fact of the optimal design. To the extent that any changes are needed, 
secondary requirements are reflected indirectly in the next design cycle 
via the LQG optimal controller synthesis procedure. Connections between 
secondary requirements and the cost function coefficient changes required 
tend to be straightforward. (For example, if activity of one control point 
is excessive, increase the cost function coefficient family weighting that 
control and those variables which are affected sensitively by that control, 
by a common factor. ) In addressing secondary requirements, primary require- 
ments must remain satisfied. Thus primary requirements act as constraints. 
If both sets of requirements cannot be satisfied simultaneously, then it 
will be necessary to modify the sensor array, provide additional controlled 
element capability or adjust bandwidth requirements downward. If simul- 
taneous satisfaction still cannot be obtained after maximum compromise, 
then it is not feasible to accomplish the mission task with the available 
resources. 

The single most important point to appreciate in connection with this 
design checklist is that only a very few design requirements are addressed 
directly in applying the LQG optimal controller synthesis procedure. Yet, 
because many of the secondary requirements are closely connected with 


32 



■bandwidth, these requirements tend to he satisfied as well. Satisfac- 
tion of the remaining design requirements is still an artful process which 
requires a full appreciation of the requirements and of their cause-effect 
connection with the dynamical physics and mathematics of the particular 
application. This is the case even though LQS optimal controller synthe- 
sis techniques are used. Therefore, the ultimate suitability of flight 
control system designs will remain keenly dependent upon the designer's 
experience. Nevertheless, modern synthesis techniques can contribute very 
significantly and effectively by providing candidate flight control system 
designs for complex tasks and plants, which are responsive to primary design 
goals. 
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SECTION III 


PROBLEM FORMULATION 


A. GENERAL 

Synthesis of practical optimal flight control systems depends upon 
appropriate formulation of the problem to an extent which we find surpris- 
ing. Careful attention to actuator representation, integrated errors, 
etc., can make the difference between practical and nonsense solutions. 

In formulating the plant, i.e., the model of the fixed or unalterable 
elements of the system, we adopt the viewpoint that we are dealing with a 
closed, stochastic system. This means, for example, that all significant 
elements of the command and disturbance input environment must be included, 
and that all variables which may be used in the cost function must be pro- 
vided. The subsections which follow attempt to s ummar ize practices we 
have found effective. Figure 9 (page 35) provides the generalized form 
into which the flight control synthesis problems are ultimately cast. 

Figure 9 can be partitioned into several subsections, as explained below. 

• Stochastic Command and Disturbance Input . The subsystem 
having matrices with S subscripts (Fig. 9a, below) is 
used to generate the co mman d and disturbance vector, yg, 
by providing shaping for the white process noise vector, w. 


*s = F s x s + fs" 

(11) 

y s = H s x g + t s 

(12) 


ts 



Figure 9a. Stochastic Command and Disturbance Subsystem 


3 ^ 





Test Points 



Figure 9- Formulation of Problem 












• Controlled Element State Variables . The subsystem having 

matrices "with C subscripts (Fig. 9^) is used to generate the 
controlled element state variables, Xq. Provision is made for 
excitation of the controlled element via control variables, u, 
commands and disturbances, yg, and white process noise, w. 

*o = r c x c + G c u + c c y s + r c w 03) 



Controlled 

Element 

States 

Xc 


Figure 9b. Controlled Element Subsystem 


• Measurements. The subsystem having matrices with M sub- 

scripts (Fig. 9c) is used to model output variables, y M , and 
sensor outputs, z. Sensor outputs may be linear combinations 
of the commands and disturbances, yg, the controlled element 
state, x c , and the white measurement noise, w. 


Mm - F M y s + % x c 

X M y M + v + tv+ -^M^OL 


z = 


Commands and 
Disturbances 
*s 

Controlled 
Element States 
X c 


( 1*0 

(15) 


^0L 



Measurements 

Z 


Figure 9c. Measurement Subsystem 
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• Other Output Variables . Variables of interest which are 

not the objective of measurement (i.e., which are not in y^) 
are included in y Q . The subsystem having matrices with 
0 subscripts (Fig. 9*3-) provides these variables. 


yo = H o x c + F oys 


(16) 



Other Outputs 

yo 


Figure 9d. Other Output Subsystem 


• Outputs . The outputs, y (Fig. 9e), consist of the variables 
which are the objects of measurement, y^, plus the other out- 
put variables, y 0 - Notice that all variables potentially of 
interest in the cost function must be included at the outset 
in formulation. 


y = 



(IT) 



Outputs 

y 


Figure 9e. Outputs Subsystem 


• Loop Opening Points . The input, y M , to every sensor (i.e., 
the kinematic portion of every measurement, z) may be inter- 
rupted individually by setting elements of the identity matrix, 
I M , to zero (Fig. 9f ) • This feature provides the ability to 
examine closed-loop system transfer functions with some or all 
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feedback loops open at the measurement point. This central 
feature permits the calculation of opt ima l system properties 
in a form consonant with the criteria and considerations used 
to evaluate conventionally designed systems. 


Outputs 

y*i 



Measurements 

Z 


Figure 9f* Loop Opening Subsystem 


• Test Points . Three types of test input points are provided, 
shown in Fig. 9g> The test input vector, tg, enables compu- 
tation of closed- and open-loop system transfer functions 
between every command and disturbance input to the controlled 
element, yg, and every state, Xg, and output, y. The test 
input vector, ty, enables computation of closed- and open- 
loop system transfer functions between every measurement noise 
input, v, and every state, xg, and output, y. (in cases where 


Shaping 

Filter 

States 

Xs 


Controlled 

Element 

States 

X C 


Test 

Measurement Test 



Figure 9g. Test Point. Subsystem 

the outer -loop function is regulation to a constant value, or 
where effective inner-loop error regulation properties are of 
interest, the test input vector, ty, can be used to supply 
surrogate command inputs.) The test input vector, top* enables 
computation of open-loop system transfer functions when outer- 
loop error feedbacks are interrupted. Transfer functions 
between error points corresponding to every command component 
of yg and every output variable y can be obtained. 
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• Cost Function. The cost function is in the form 


J c = E[y' Q E y + u'E e u] 


(18) 


where Qg and Eg are diagonal matrices. The fact that 
and R F matrices are diagonal results in a cost function 
integrand which is a weighted sum of squares of the ele- 
ments of y and u. This restriction stems from our belief 
that the cost function should be interpreted as an index 
of performance, and to this end should consist of a weighted 
sum of squares of variables indicative of performance for 
each particular application. This set of variables must be 
contained in the output and control vectors, y and u, by 
virtue of their respective definitions. 

The separability principal is used in obtaining the optimal 
controller solution. This leads to the use of two separate 
cost functions, J F and J F : 


J R 



(y'Q y + u’R u) dt 


( 19 ) 


J F = E[(x c -x c ) , (x c — x c ) + (x s — x s ) f (x s -x g )] 


Minimization of these functions permits the separate solution 
of the optimal regulator and filter problems, as an alterna- 
tive to conib ined solution using Jq. 

The process noise vector, w, and measurement noise vector, v, are 
independent, zero-mean Gaussian white noise processes. Their intensities 
are given by the diagonal Q and R matrices, respectively, i.e., 


E[w] 

= E[v] = 0 

(20) 

c Pww' 

- \ 

(21) 

<Pw r 

\ 

= R 

\ 

(22) 


To this point, the matters dealing with optimal control are more or 
less standard. Here, however, we need to recognize that many measure- 
ments in flight control systems are, for all practical purposes, noise 
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free. Thus, some or all of the measurement noise components may have 
zero intensity. This, in turn, implies an ability to solve the singular 
Kalman filter problem. 

In the following paragraphs the definitions of the partitions of the 
problem formulation shown in Fig. 9 are specialized for the flight con- 
trol system application. 

B. RECOMMENDED FORMULATION FOR 
FLIGHT CONTROL SYSTEMS 

1. Stochastic Commands and Disturbances 

Stochastic commands may be represented by first- or second-order shaped 
white noise. First-order shaping can be used to represent maneuver require- 
ments, guidance source inaccuracies, or inadvertent pilot inputs. One form 
of the second-order shaping filter equations is: 


X S1 = 

-CDgi X S1 + 


(23) 

X S2 = 

- 0^2 x S2 + 

\/2chg2 ^2 

(24) 

ysi = 

X S1 + X S2 


( 25 ) 

^w-| = 

2 

X S1 ' 

^2 = °XS2 

( 26 ) 


To obtain first-order shaping merely omit the equations for x S2 and Q^. 

The output, y g1 , is the sum of two independent first-order shaped 
processes. One process might represent maneuvering co mman ds and have 
moderate level and low bandwidth. The second process might represent 
guidance source inaccuracies, and have low level and moderate bandwidths. 
Relatively little numerical data are available to guide selection of level 
and bandwidth. Some data are contained in Ref. l4 and further suggestions 
appear in Appendix B. Stochastic disturbances are mainly due to gust 
effects. Gusts may be modeled by the Dryden form of the power spectral 


1+0 
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densities (Ref. 1 5 ) . These power spectral densities correspond to first- 
and second-order shaped white noise. The longitudinal (u g ) and effective 
rolling (Pg) gust components are first-order forms; 


*s 

= 

-^ x s + v^ w 

( 27 ) 

ys 

= 

x s 

(28) 

Qw 

= 

ai s , x s = u g or p g 

( 29 ) 

Side and normal gust components are second-order forms ; 


X S1 

= 

-CDgXgi + W 

(30) 

X S2 

= 

"“S X S2 + x sl 

(31) 

y 

= 

7572 * S1 - ( 7? - i)x S2 

(32) 

Qw 

= 

°ys ' y s = v g or w g 


In each case represents the characteristic frequency of the particular 

gust component. Second-order shaping may be approximated by first-order 
shaping for simplicity. The form is ; 

X S 

= 

-1.59i)-a^x s + y2(l.594)a^ w 

(33) 

ys 

= 

x s 

(34) 

Qw 

= 

4 S > = v g or w g 

(35) 


This approximation preserves the half -power frequency, 1.594ocg, and the 
rms value, oy^, of the second-order power spectral density. 




Additional background and numerical data can be found in pages ^17- 
46l of Ref. 15, pages 65^-660 of Ref. 10, and in Appendix B. 

2. Actuators 

Actuation dynamics, activity, displacement limits and rate limits 
almost always provide the main limitations on automatic flight control 
system performance and frequently contribute to performance limitations 
in manual control. For this reason inclusion of an actuation dynamic 
model at some point in the optimal flight control system design cycle is 
essential. Fortunately, actuation dynamics need not be modeled in detail. 
Simple approximations capture the essential effects. 

Consider the NASA VSTOLAND UH-1H helicopter actuation system (Ref. 16 ) 
in Fig. 10. An approximation adequate for design purposes is shown in 
Fig. 11. The approximation in Fig. 11 results in a reduction in order for 
the actuation system model of three. The effective actuation lag is com- 
puted by the method given in Appendix B. We have found it to be good prac- 
tice to represent the effective actuation lag by its open-loop element, 
i/s, and supply the actuator feedback gain via the optimal design process. 
This practice results in accurate representation of the mean square (boost 
actuator) control deflection and control rate in the cost function. 


SERIES ACTUATOR 



Figure 10. Actuation System Block Diagram 
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SERIES ACTUATOR 



Figure 1 1 . Approximation to Actuation System 

Other quantities in the actuation system are frequently subject to 
physical limitations. Series actuator displacement and parallel actuator 
rate are usually in this category, since both quantities are ordinarily 
limited for safety-of-flight reasons. These quantities can either be 
evaluated directly using variables existing in the plant formulation, or 
by including additional state and output equations as described in Appen- 
dix B. 

3 . Airframe 

The airframe can be represented by uncoupled longitudinal and lateral- 
directional sets of equations. Longitudinal equations may be ■written con- 
veniently in terms of states u, w, q, 0, S-g, S c , and the lateral- 

directional equations in terms of states v, p, r, cp, & A , 5p, The 

last state in each set represents the stabilizer bar state. Crosscoupling 
between longitudinal and lateral-directional systems has been found to be 
relatively unimportant in connection with flight control systems design. 
Specific sets of state equations for the UH-1H used in our examples are 
given in Appendix B. 

Airframe models which are more, or less, detailed than the one described 
above may be appropriate for a particular application. The problem most 
often faced is that of obtaining a simplified model for flight control 




system design purposes from a very high-order, detailed model of the 
airframe. We shall assume that the high-order model is linear and sta- 
tionary. Simplification of this model is accomplished by making low- 
frequency approximations to very high-frequency modes. This process is 
called "aggregation, M e.g.. Ref. 17. "High" frequency is interpreted 
relative to the usual range of flight control system modes, typically 
the hand 0.0 to 50.0 rad/sec. There are at least two methods for develop- 
ing the simplified model. In the first (e.g.. Refs. 18 and 19), a modal 
representation of the system is constructed in terms of its eigenvalues 
and eigenvectors. Derivatives of modal coordinates having eigenvalues 
with modulus larger than 50.0 are approximated by zero. The approximated 
representation is reconstituted to the original coordinates after algebraic 
equations* of the system are eliminated by substitution. These algebraic 
equations become additional output equations for the reduced-order model. 
Dependencies of the measurement and original output equations upon the 
variables of these additional output equations must be eliminated by sub- 
stitution. A compact exposition of this method is given on pages 3^-36 of 
Ref. 18. 

The second method is in the spirit of the first but does not require 
construction of a modal representation of the system. It is frequently 
the case that modal subsystems which couple only weakly with the rest of 
the system can be identified by inspection of the equations of motion, 
(inspection may be by comparing eigenvalues of the subsystem with eigen- 
values for the total system, for example. If subsystem eigenvalues are 
nearly equal to corresponding total system eigenvalues, then coupling is 
weak.) Weakly coupled subsystems having eigenvalues larger than 50*0 rad/ 
sec may be replaced in the equations of motion by approximate algebraic 
equations which represent the static (low-frequency) characteristics of the 
subsystem. This approximation has been applied in Ref. 20 and in Appen- 
dices C through H of Ref. 21^ for example. 


*These algebraic equations result from those equations for which deriva- 
tives of modal coordinates are set to zero. 
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k. Kinematic! 


It is essential that all kinematic relationships which are in differ- 
ential equation form be included in the plant. These may include not 
only obvious kinematic relationships such as 


h = u sin ©o — w cos © Q + Vp Q 0 

(56) 

\Jr = r /cos ©Q 

(5T) 


but also error integrals which may be required in the cost function such 
as 


h c - h 

(38) 

*c “ * 

(39) 


5- Measurements 

Equations for the measurements assumed available must be accurate with 
respect to sensor location and orientation effects. Sensor noise charac- 
teristics must be represented with reasonable fidelity. Most flight con- 
trol system measurements are free from broadband (i.e., "measurement'') 
noise as a practical matter. Accelerometer measurements are a possible 
exception, however, especially when installed in a vibratory environment. 
Some flight control measurements, such as glide slope error for example, 
have correlated noise (e.g.. Ref. 22 , Appendix B) . It is essential for 
proper formulation of the problem that correlated noise be represented by 
shaped process noise rather than by white measurement noise. 

A list of the conventional sensors which have been found useful in 
flight control systems is given in Table 7- A summary of special consid- 
erations in using, orienting and locating certain of these sensors follows 



TABLE 7- CONVENTIONAL FLIGHT CONTROL SENSORS 


Longitudinal Axes 


q. 

0 

e 

h 

h 

R 

R 


Body axis pitching velocity 
Pitch attitude rate 
Pitch attitude 

Rate-of-climb (barometric or haro-inertial) 
Altitude (barometric or radar) 

Glide slope deviation (angular measure) 

DME or IRS derived slant range 
Derived DME or IRS slant range rate 
IRS ground speed 


Body-fixed longitudinal acceleration 
az Body-fixed normal acceleration 

(a — ag) Aerodynamic angle of attack 

V A g Airspeed 

5g/5 e Longitudinal cyclic pitch or elevator deflection 

5c/&t Collective pitch or power lever deflection 


Lateral 


P 

9 

9 

r 

Uf 

v 

( ^ - X ) 

4 

(P — Pg) 

^A^a 

5 R /8 r 


Body- fixed rolling velocity 
Bank angle rate 
Bank angle 

Body -fixed yawing velocity 
Heading rate 

Heading (magnetic or gyro compass) 

Localizer or VOR deviation (angular measure) 

IRS drift angle 

Body-fixed lateral acceleration 

Aerodynamic angle of sideslip 

Lateral cyclic pitch or aileron deflection 

Tail rotor collective pitch or rudder deflection 



Special Considerations 


ai: Body-fixed longitudinal acceleration is sensitive to orientation 

and location effects in the airframe. Location effects are secondary, 
but near c.g. alignment is preferred. Orientation effects are important 
and are usually dictated by alignment with the flight path in the flight 
regime where close groundspeed or airspeed control is desired. In these 
instances it is often desirable to use "longitudinal acceleration inde- 
pendent of pitch," i.e., let a x = a x — g sin 0 be the effective measure- 
ment, by combining a x and g sin 0 at the outset. This results in an 
approximate measurement of Vy w ^ en both quantities are referenced to/ 
aligned with the reference flight path. 

ay: Lateral acceleration is sensitive to location effects. This is 

frequently used to advantage in flight control systems in order to create 
a surrogate (3 measurement. By locating the accelerometer a distance 

= “ y 5h/ n 5r aliead - the aircraft c.g.,* one obtains a^- and p transfer 
functions for 5^ inputs which are identical over a broad range of frequen- 
cies except for a scale factor. 

a z : Normal acceleration is sensitive to both orientation and location 

effects. Preferred values depend upon the particular application. 

a z may be used to extend the bandwidth of the barometric rate of climb 
measurement which has an inherent lagging characteristic. In this case, 
the accelerometer should be located close to the aircraft c.g. and should 
be oriented so that its sensitive axis is near vertical at trim for the 
flight regime where close rate-of-climb control is desired. This same 
location and orientation is appropriate when the a z measurement is fed 
back to collective pitch in order to alter the apparent mass of the air- 
craft for flight path control. 

Normal acceleration can also be used to create a surrogate a measure- 
ment via proper location. By locating the accelerometer a distance 


*That is, at the instantaneous center of rotation with respect to the 
tail rotor collective control. 
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= Zgg/Mgg ahead of the aircraft c.g.,* one obtains a z and a transfer 
functions for 5p inputs -which are identical over a hroad range of fre- 
quencies except for a scale factor. This use of normal acceleration is 
frequently impractical for rotorcraft application because Zgp may change 
by more than an order of magnitude between hover and maximum speed and is 
usually positive. 

For all of the above uses of normal acceleration, it is necessary to 
wash out, bias out or cancel the gravitational component of the accelera- 
tion measurement. This requirement is not evident in the control design 
problem because linearized perturbation equations are typically used wherein 
the constant gravitational component does not appear. 

q, r, 9, \jr: Airframes requiring large amounts of artificial damping 

or which must execute gross maneuvers must make careful use of superfi- 
cially equivalent measurements, i.e., of q and 0 or of r and {jr. For exam- 
ple, an airframe requiring heavy pitch axis damping will lose that damping 
augmentation in steeply banked turns if 0 is used in lieu of q since 

© = q cos 0 — r sin $ 

On the other hand, use of q in lieu of 0 when the requirement is for ver- 
tical path equalization will result in opposition of steady turns since 

q = © cos $ + y cos 0 sin $ 

(a— ag), (|3 — Pg) : Measurements of aerodynamic angle of attack and 

angle of sideslip are frequently unsatisfactory in rotorcraft applica- 
tions. These measurements, when used, require correction for installation 
and location effects. 


*That is, at the instantaneous center of rotation with respect to the 
longitudinal cyclic pitch control. 
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6. Outputs 


Numerous linear combinations of the state variables in a given appli- 
cation are of interest because of their relevance to the physics of the 
problem, to performance or flying qualities specifications or because of 
their potential use in the cost function. Our recommendation is that all 
measurement kinematic components (y^ in Fig. 9 ) , and candidate measurement 
quantities, all candidate quantities for use in the cost function, all 
quantities which may encounter physical or specification limits be included 
in the output vector, y. Virtually all of the variables in Table 7 may be 
included in the output vector. Additional variables associated with the 
actuation system such as series actuator displacements and parallel trim 
actuator rates must be included when appropriate to the particular appli- 
cation . 

C. EXAMPLE FORMULATION 

Application for design of a longitudinal flight control for the UH-1H 
helicopter at hover is used for illustrative purposes. Pertinent data for 
the primary design objectives are summarized in Table 8. Vector definition 


TABLE 8 

UH-1H LONGITUDINAL FLIGHT CONTROL SYSTEM 


Flight Condition: 


Hover 


Functions : 


Rate-of-climb command 
Groundspeed hold 


Bandwidths : 

rad/sec 

Cyclic (DB) 

25-73 

Collective (DC) 

25.73 

Pitch (TH) 

2.0 

Rate-of-climb error (HDE) 

1.0 

Integral of HDE (HDl) 

0.82 

Groundspeed error (XD) 

0.5 

Integral of XD (XDl) 

0. 1 



for the plant, measurements and outputs are given in Table 9. Note that 
each variable has an associated computer mnemonic, as indicated below 
each vector- With the exception of the state variables, the mnemonics 
are input to the software (Ref. 23) by the user, and are used to produce 
the data from the software in report-ready form. In the case of the 
state variables, however, the mnemonics are internally generated by the 
software. The names X01, X02, ..., X05 are reserved for the shaping filter 
states, while X06, X07, ..., XT 5 are reserved for the controlled element 
states. Most matrix data in the remainder of this report were produced by 
the software just described, and the mnemonic conventions apply through- 
out. 

Throughout this report, it is important for the reader to maintain the 
distinction between variables which are states and variables which are 
measurements, even though the same variable names are used. This confusion 
could arise with the variables q, 0, 5 B and Sq. For instance, we see that 
the state variable 63 (the actual longitudinal cyclic deflection) is defined 
as the integral of the control variable 63, via the equation x = Fx + Gu + rw. 
Once the controller has been defined, though, we will see that the measure- 
ment variable 5 b (which may only exist as a quantity inside a processor) is 
not the integral of the control variable 5 b* We will continue to point out 
this distinction in the text. 

The software and its associated mnemonics do not suffer from these ambi- 
guities. For example, the state variables q, 0, 63 and- £>c 3X6 X09, 

X10, and XII, whereas the measurements q, 0, 5 b and &c are TH, DB and 
DC (see Table B-l). 

Table 10 presents the literal form for the F, G, F, and Hr matrices 
for the example problem. These matrices are composites of the matrices 
which comprise the problem formulation (Fig. 9) as indicated in the table. 
The numerical data for the example problem in terms of these problem for- 
mulation matrices are given in Table 1 1 . The numerical data for the F, G, 

F and Hr matrices are given in Appendix C. 

Throughout the report it will be useful to present factored transfer 
functions. The notation used is explained in Table 12. Table 13 presents 
selected factored transfer functions for the example controlled element. 


50 





TABLE 9- EXAMPLE PLANT DATA 

VECTORS 



Shaping 

Filter States 


X S 

= jUg w g 

[ 


= jXOI X02 

XO?} 


Controlled Element States 


f ^ 

x c = 

i u 

w 

q 

. 

0 

&B 

5c / h e dt 

■^BAR 

/ x dt [ 

= 

Jxo6 

xot 

xo8 

X09 

X10 

X1 1 

XI 2 

X13 

X11+ 

i 

Process Noise 











■w’ = 




w u 

w w 


W 5 B 

w s c v I he it 

¥y 

X BAR 

W / x dt f 

= 

{ PUG 

PWG 

FHC 

PU 

PW 

pq PTH PDB 

PDC 

PHD 

PXB 

PXD } 

Controls 












u’ = 

) Sb 

&C 

i 









= 

{dbd 

DCD 

( 









Shaping Output 











T ^ 

y S = 

{ u g 

w g 

! 









= 

jUG 

WG 

HDC ! 

[ 








Measurements 











z' ^ 

1* 

q 

0 

5 B 

Sc 

^e 

/ h e dt 

X 

/ x dt [ 



= 

|bd 

Q 

TH 

DB 

DC 

HDE 

HDI 

XD 

XDI [ 



Outputs 












y* = 


q 

0 

5 B 

5 C 

^e 

/ h e dt 

X 

/ x dt 

a 

u AS e J 

= 

I HD 

Q 

TH 

DB 

DC 

HDE 

HDI 

XD 

XDI 

AOA 

ASE } 
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TABLE 10. LITERAL FORMS FOR F, G, F AM) H R MATRICES 
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(continued on following page) 



TABLE 10. (Continued) 
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(concluded on following page) 



TABLE 10. (Concluded) 



5b 


a XDI used in hover example ; ASI used in TOO kt cruise. 

^XD used as measurement in hover, output at 100 kt,- ASE used as measurement at 100 kt, output in hover. 
c Used to investigate acceleration measurements in 100 kt case. 





TABLE 11. (Continued) 
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TABLE 11. (Continued) 



11. (Continued) 
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TABLE 11. (Concluded) 
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998 0. 705E-01 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1 AS E 



TABLE 12. NOTATION FOR FACTORED TRANSFER FUNCTIONS 


CASE: UH1H HOVER 122L0NG 31 -JAN- 79 CONTROLLED ELEMENT TF'S 

DENOMINATOR: 

1.0000 — — - — — ■ — — 

( .00000 ) ( .00000 ) ( .00000 ) ( .00000 )\ 

( .38494 ) J 

(( .10539 , .18934 , .19576E-01, .18833 )) j 

(( .26279 , .92717 , .24365 , .89459 )) r 

< .ii863E-oi> 


Title 


High-frequency gain 
Real poles 


Complex poles 
Low-frequency gain 


NUMERATOR: TH/DBD 


Numerator title 


-.16910 — 7 

( .00000 ) ( .00000 ) ( .00000 ) {-.79065E-02) 

( .33300 ) ( .39184 ) 

< .17445e-03> — ■ • 


High-frequency gain 
Real zeros 
Low-frequency gain 


_e_ _ -■ 1691 s5(s- .0079) (s+ • 333)(s+ .39184) 

5 B s u ( s + . 38494) [ s 2 + 2( .10339) ( . 18934 )~s + ( . 1 8934) 2 ] [ s 2 + 2( .26279) ( - 9271 7) s + ( -927i 7) 2 ] 


-. 1691 (s + .533) 

[s 2 + 2(. 10339) ( -18934)s+ (. 8934) 2 ][s 2 + 2(. 262 79) (.9271 7)s + (.9271 7) 2 ] 


Alternately, the notation can be considered shorthand for expressing factored polynomials; 


K 

(z) 

((?, -^2 )) 

<Kzo^> 

corresponds to 

K(s + z)(s 2 + 2£cus + u?) 


Note that the complex poles and zeros have four components; 

((damping ratio, undamped natural frequency, real part, imaginary part)) 
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TABLE 13 


SELECTED CONTROLLED ELEMENT TRANSFER FUNCTIONS 


CASE: UH1H HOVER 122LONG 31-JAN-79 CONTROLLED ELEMENT TF’S 

DENOMINATOR: 

1. 0000 

{ .00000 ) ( .000130 ) { .00000 ) ( .00000 ) 

( .33494 ) 

(< .10339 , .18934 , .19576E-01, .18833 )) 

(( .26279 , .92717 , .24365 , .89459 )) 

< . 11863E-01 > 


NUMERATOR: TH/DBD 

-.16910 

( .00000 ) ( .00000 ) ( .00000 ) (-. 79065E-02) 

( .33300 ) ( .39184 ) 

< . 174 4 5E-0 3 > 


NUMERATOR: TH/DCD 

-. 33000E-02 

( .00000 ) ( .00000 ) ( .00000 ) ( . 10854E-01 ) 

( .33300 ) (-11.270 ) 

< ,13442E-03> 


Shorthand notation is used to express factored poly- 
nomials. For example, 

K(s + z)(s 2 + 2f;cos + o£) 

corresponds to 

K 

(z) 

((£, (Q, CDv/-! _^2 )) 
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SECTION IV 


CONTROLLER SYNTHESIS 

The controller optimal synthesis can be divided into three distinct 
steps: 

• Filter-observer synthesis 

• Regulator synthesis 

• Combination of the above into the controller 

This section starts by reviewing some important linear stochastic control 
facts in order to define the synthesis problem and solution. The solu- 
tion of the singular filter problem, and the implications of requiring a 
solution which does not use differentiation, are considered. This is fol- 
lowed by specific recommended procedures for the synthesis of flight con- 
trol systems and illustrative examples of their application. 

A. STOCHASTIC OPTIMAL CONTROL PACTS 

Stationary plant state and measurement equations and process and 
measurement noise characteristics for the system which is the object of 
optimal control are summarized in Table 14. The notation and definitions 
used follow those customary in the literature (e.g.. Ref. 1 ) . We have 
made specific provision in Table 1^ for noise-free measurements via the 
zg partition of the measurement vector. 

It is well known (e.g., Ref. 6) that the linear stochastic optimal 
control solution with respect to the cost function 

j c = E[y f Q R y + u’ R r u] (^0) 

where the plant outputs are 

y = H R x 
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TABLE 14. LINEAR PROBLEM MODEL 


States 

x = Ex + Gu + Pw , x(O) = x 0 
Measurements 



z-| = H^x + 

.0 

z 2 = HgX + 

Boise Characteristics 

E[w] = 0 , E[v] = 0 

E[w(t.,) w f (t 2 )] = Q6 ( t 2 — t -j ) 

E[v(t-| ) v' (t 2 ) ] = R5(t 2 — t-|) 


can be partitioned into solution of two separate problems. These are 
the (Kalman) filter and regulator problems. The combined solution, shown 
in Fig. 12, results in the linear optimal stochastic controller. 

The filter problem is solved to minimize the mean square estimation 
error via the cost function: 


J F = E[ (x — x) * (x - x)] (41) 

Solution results in the filter gain matrix, K. In formulating the problem 
we have assumed that some or all of the measurements are noise free. Thus, 
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n 



Figure 12. Form of Linear Stochastic Optimal Control Solution 


the measurement noise intensity matrix R may he singular. This, in turn, 
implies difficulty in obtaining the filter solution by computational methods 
which depend on the existence of B \ Appendix A presents an alternative 
solution technique which yields the required partitions, Kq i and K-jg, 
the filter matrix, K, for this singular case. Solution of this singular 
filter problem is a key contribution for optimal flight controllers because 
most measurements in flight control applications are practically noise free. 
It also enables us to obtain lower -order controllers which are necessary 
for practical flight control systems. 

The regulator problem is solved to minimize the cost function: 


J R 



(y' Qr 7 + u’ R r u) 


dt 


(te) 


Solution results in the regulator gain matrix, C. To be consistent with 
the physical constraints on control use in flight control applications 
we shall require R R to be non-singular. The central issue in the regu- 
lator problem is the selection of appropriate Qr and Rr matrices. The 
selection must produce adequate performance yet not waste control resources. 
A methodology for their rational selection is another key contribution of 
this research. 
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As depicted in Fig. 12 the "controller" is defined "by the filter gain 
matrix (k), the regulator gain matrix (c), and the equation; 

u = — Cx (k3) 

In the flight control application we elect additional restrictions. 

We shall require the plant to he stabilizable with respect to the avail- 
able control points, and detectable with respect to those elements of y 
which have non-zero weighting in Qr (Ref. 8) 

The cost function shall be a weighted sum of mean square values of 
variables having operational significance for each application. There- 
fore, Qr and Rr are restricted to be diagonal matrices because the meas- 
urement matrix. Hr, can be used to generate all variables having opera- 
tional significance except for the controls. However, the controls are 
available in their operationally significant form by virtue of their 
definition. Further, we shall require Rr to be positive definite. 

Singular filters have order (number of states) which is less than or 
equal to (n — mg), where n is the number of plant states and mg is the 
number of independent noise-free measurements. When the order is less 
than (n — mg), differentiation of the noise-free measurements may be 
required. Differentiation is undesirable for the flight control applica- 
tion because very low level, broadband noise components of measurements 
are a virtual certainty even if negligible from all other viewpoints. In 
order to make differentiation of the noise-free measurements unattractive 
as an optimal solution we resort to the artifice of augmenting the process 
noise excitation of the plant . The added noise is specifically chosen to 
cause the derivative of every noise -free measurement to contain an inde- 
pendent process noise component. Mathematically, this is equivalent to 
requiring Tg QTg to be nonsingular (Appendix A).* The levels of the 


*This solution for the filter problem having noise -free measurements 
was reported by Uttam and O’Halloran, Ref. 25, pp. 327-332, but is not 
recognized as the optimal filter. 
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augmenting noise can “be made very low, although numerical roundoff in com- 
putation imposes lower bounds. When r 2 QTp is nonsingular, the filter order 
is exactly (n— m^) and the filter can be realized using that number of inte- 
grators. The levels of the augmenting process noise can be manipulated to 
affect the closed-loop filter eigenvalues. This can be used to advantage 
in controlling the range for the closed-loop filter eigenvalues. It is 
also the reason we refer to the singular filter as the "filter-observer." 

Although augmenting the process noise excitation of the plant is an 
artifice, it can be rationalized. The stationary, linearized, perturba- 
tion plant equations we tend to regard as a true model are in fact an 
approximation. Nonlinear, nonstationary effects and small terms in the 
equations are neglected, as is the uncertainty in the values for the coef- 
ficients in terms we retain. One may look upon the augmenting process 
noise as an aggregated representation of these conveniently neglected 
effects. 

When T 2 QFp is nonsingular, the optimal filter-observer has the struc- 
ture shown in Fig. 13. (Figure 13 assumes for simplicity, and without loss 
of generality, that the plant state definition includes the noise-free 
measurements, Z 2 , as the X£ partition. ) 



Figure 13. Optimal Filter -Ob server Structure 
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Given the regulator gains, C, in terms of the same state vector, the 
controller structure and equations can he shown as in Fig. 1^-. The coef- 
ficient matrices in Fig. lA are given by 


Ap - [F n “ K T2 F 21 “ K 1 1 H 1 1 + (G-J -K q2 G2) C F^ 

Bp = [K-ji, (F-j i — K-]2F21 — K-] -|H-] i )K-|2 + F-)2 — K'i2^'22 — ^-1 1^-12^ 



Figure 1^-. Form of Controller Solution 
(Filter-Observer + Regulator) 

This block diagram, and the matrices therein, corresponds to the 
simplest form of the generalized controller equations. This form, includ- 
ing the Ap, Bp, Cp and Dp matrices, is used in the optimal design soft- 
ware. Figure 15a shows this same controller, but here the block diagram 
structure more clearly delineates the filter-observer and regulator 
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a) Generalized Optima / Controller . c) AH Measurements Noise- Free , 



Figure 15. Form of Controller Solution — Special Cases 






















components . Two limiting cases are also of interest. When all the meas- 
urements are noisy the controller reduces to the block diagram of Fig. 15b; 
the filter -ob server reduces to the familiar Kalman filter, while the regu- 
lator is unchanged. With all measurements noise-free, and with fewer meas- 
urements than states, the filter -observer becomes a pure observer, used to 
est ima te those states for which measurements are unavailable (Fig. 15c). 

If the n umb er of measurements equals the number of states, however, an 
observer is not needed at all/ the noise -free measurements constitute alter- 
native states, and are subsequently fed to the usual regulator (Fig. 15<i)* 

The recitation of stochastic optimal control facts given above estab- 
lishes the mechanics for obtaining solutions. The following subsections 
present the techniques for purposeful application and illustrative exam- 
ples. These subsections focus on two key questions: 

• How should one choose the diagonal Qp and Rp 
matrix elements for practical designs? 

• How should one augment the process noise input 
to avoid differentiators in the filter solu- 
tion? 

Answers to these questions complete the data required to synthesize an 
optimal flight control system responsive to primary design goals. All 
other data required are obtained from the physics of the problem. 

B. FILTER-OBSERVER SYNTHESIS 

In defining the shaping and controlled element states, the measure- 
ments, and the outputs in preparation for optimal controller design we 
have attempted to use experience and foresight to make the actual optimal 
control calculations straightforward. For example, we have chosen to 
include the integrals of variables when the indicial error response for 
those variables must approach zero as time approaches infinity. Also, 
care has been taken to insure that all pertinent variables are included 
in the output vector. 

Still, there are certain fine points which remain elusive. Foremost 
among these is the definition of a desirable group of measurements. 
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together with the corollary question of the adequacy of the associated 
sensor array. The synthesis of the filter-observer permits these points 
to be resolved. In particular, the synthesis provides a definitive meas- 
ure of the sensor array performance for a particular measurement set. The 
order of the controller will also be finalized because it is determined by 
the order of the filter -observer. 

It is recommended that the filter -observer synthesis be addressed 
before the regulator synthesis. This is mainly for the reason that once 
the rms state estimation error and filter-observer eigenvalues are in 
acceptable ranges (specified later) the designer need not reconsider the 
filter -observer design in the later stage of optimal closed-loop system 
assessment. (The same cannot be said for the regulator solution.) Fur- 
thermore, controller transfer functions are imme diately available at each 
cycle of the regulator solution when the filter-observer solution is already 
in hand. In this way, the controller tends to proceed in a top-down fashion. 

The filter -observer has three main functions: 

• It defines the quality of the instrumental complex 
used to make measurements, via the rms state estima- 
tion error. This is the primary reason cited above 
for beginning the controller design with the filter- 
observer synthesis. A stage is set upon which several 
alternate sensor arrays may compete, and the best one 
be chosen. 

• It provides the best estimate of the total state 
vector, given the available measurements, to use as 
input to the regulator. 

• It is the only place where the detailed properties of 
the noise and disturbances can affect the controller 
design. 

This final function is perhaps the most important. In essence, the noise 
and disturbance properties directly affect the closed-loop system in that 
the filter-observer poles are closed-loop system poles. Filter -observer 
poles are a function of the process and measurement noise levels. These 
poles are the equivalent of the compensator poles in a conventional design 
in that they adjust the controller behavior as a function of the distur- 
bance environment. This is analogous to the adjustment of washout and 
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complementary filter “break frequencies appropriate to the command, dis- 
turbance and noise environment in conventional design practice. In the 
course of the filter-observer discussion we will see how these same poles 
affect the controller in the command environment. 

Steps in the filter-observer synthesis are summarized in Table 15 
for ready reference. These five steps are elaborated in the following 
paragraphs and are illustrated by examples. 

TABLE 15. FILTER -OBSERVER SYNTHESIS 


Specify process and measurement noise levels 
based on physics (Q and R) 

Augment process noise to get independent noise 
on the derivative of every noise -free measure- 
ment 

Compute closed-loop filter eigenvalues and rms 
estimation error 

Change sensor array if rms estimation error is 
extreme 

Adjust augmenting process noise if eigenvalue 
range is unacceptable 


1 . Physical Components of Process 

and Measurement Noise 

The process noise components which are used to generate shaped command 
and disturbance inputs to the controlled element have levels which are 
determined by the physics of the operating environment. In practice, 
these levels may be definitively established, as in the case of gusts, 
or may merely be estimates, as in the case of some commands. Here, we 
assume the noise levels which are features of the environment are known. 
Measurement noise components should also be based upon the physical situa- 
tion. For example, it is reasonable to assume that rate and attitude 
gyros and actuator position pickoffs provide noise-free measurements. 
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Measurement noise on body -mounted accelerometers may or may not be 
regarded as significant depending upon factors specific to the applica- 
tion. If significant, measurement noise levels may be inferred from the 
threshold and linearity specifications for the particular accelerometer 
instrument. In the case of Microwave Landing System (MLS) guidance, for 
example, the measurement noise (power spectral density) level can be cal- 
culated from either the levels of error budget power spectral densities 
(e.g.. Ref. 25) for the broadband sources, or from measured power spec- 
tral densities. 

2. Augmenting Process Noise Components 

The requirement to differentiate any noise -free measurement can be 
avoided if the derivative of every noise-free measurement contains an 
independent component of process noise.* One means by which this condi- 
tion may be satisfied is to simply include an identity partition in the 
controlled element r c matrix (refer to Table 11). This identity parti- 
tion is in addition to any columns of the Tc matrix which arise from the 
problem physics. A corresponding null partition must be included in the 
shaping r s matrix to insure dimensional compatibility. Subject only to 
the assumption that no command or disturbance input is a noise -free 
measurement, this provision and the ability to specify the levels of the 
augmenting process noise arbitrarily are sufficient to provide an inde- 
pendent component of process noise in the derivative of every noise-free 
measurement. That this is so is verified by the following equations. 


^Alternatively, there will never be a requirement to differentiate 
any measurement containing measurement noise. However, adding measure 
ment noise for the sole purpose of avoiding differentiation of a meas- 
urement is unattractive because it increases the order of the filter- 
observer. The order increases by one for each measurement so treated. 
Practically, this would require an increase in controller complexity. 
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z = ••• + [F M H s r s + H M r c ] w + ••• (^6) 

Inspection of [FjyjHgrg + re] reveals which augmenting process noise 
components must have non-zero levels . This inspection proceeds row-by-row 
for only those rows which are associated with noise -free measurements. 

An array formed of these rows can be constructed.* This array must have 
rank mg, where mg is the number of noise-free measurements. Up to n — mg 
columns of this array may be eliminated without reduction of the array rank. 
Any set of columns which are coefficients of the augmenting process noise 
and which can be eliminated without reducing the array rank, may have cor- 
responding zero levels in the diagonal process noise intensity matrix, Q,. 

The values for the remaining augmenting process noise levels must be greater 
than zero, but are otherwise arbitrary. 

The above approach to augmenting the process noise has been adequate 
for all applications made to date. However, it is clear that augmenta- 
tion of the Fq and Tg matrices by identity and null partitions, respec- 
tively, is not required. Indeed, these augmentation partitions may be 
completely arbitrary subject only to the restriction that T2 have full 
rank. 


3. Filter Eigenvalue 8 and RMS Estimation Error 

The third step, that is, the computation of closed-loop filter- 
observer eigenvalues and rms estimation error, is accomplished using the 
software package (Ref. 23 ) . 


*This array is computed by the software as T2* 
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Changes to Sensor Array 


Next, the rms estimation error is compared with target values for 
closed-loop system performance. Since rms estimation error establishes 
a lower hound on achievable closed-loop system performance, it is essen- 
tial that no target rms performance value he exceeded hy its correspond- 
ing estimation error component and that significant margins exist. If 
target values are exceeded, this is an indication that the sensor array 
is inadequate with respect to variables measured, quality of measurement,* 
location/orientation of sensors, or any combination of these factors. 
Changes to the sensor array with respect to one or more of the above 
characteristics should be made before continuing the controller synthe- 
sis. On the other hand, if rms estimation error is exceedingly small with 
respect to target values, this may be regarded as an invitation to explore 
use of fewer or lower quality* sensors. Use of fewer noise -free measure- 
ments, of course, requires use of a correspondingly higher-order filter- 
observer. This tradeoff aspect must be borne in mind when exploring the 
possibility for reducing the number of sensors. 

5- Adjusting Augmenting Process Noise 

The final step is to inspect the eigenvalue range for the closed-loop 
filter-observer. Eigenvalues with very large modulus (say, greater than 
50. rad/sec) indicate that large gains will be required In filter implemen- 
tation and may indicate excessive dynamic range requirements. Eigenvalues 
with very small modulus (say, less than 0.1 rad/sec) indicate that very 
slow closed-loop response modes arising in the filter-observer may be 
present. These modes may be excited by command or disturbance inputs, or 
as the result of differences between the real-world plant and the plant 
model for which the filter -observer is designed. Some closed-loop filter 
eigenvalues may be altered when some or all of the measurements are noise- 
free. This can be accomplished by adjusting the augmenting process noise 


*Lower quality measurements can only be explored when noise effects are 
explicitly modeled. 



via changes in the appropriate columns of the r and Q matrices. Our 
experience in application is that changes in the augmenting process noise 
levels via the diagonal Q matrix axe usually sufficient to obtained closed- 
loop filter eigenvalues in the desired range. 

6. Filter -Ob server Synthesis Examples 

This subsection provides illustrative examples for the filter-observer 
considerations addressed in previous paragraphs. The application to be 
described in detail is to the longitudinal dynamics of the UH-1H for the 
hover flight condition. Table 16 is a data summary for the final filter- 
observer design. 

Consider the elements of the diagonal process noise intensity matrix, 

Q, in Table 1 6 . The associated mnemonics have P prefixes to denote pro- 
cess noise. The remainder of each mnemonic has been used to designate the 
differential equation each process noise component forces. The first three 
element values are fixed by the physics of the assumed operating environ- 
ment. The values given are , and cru , respectively, because the 

shaping filter gains, y / 2o>^, and -/2u> are contained in Tg. 

The remaining nine elements of the Q matrix diagonal are augmenting process 
noise intensities. Notice that one of these intensities is zero. 

Next, consider the elements of the diagonal measurement noise inten- 
sity matrix, R. The associated mnemonics indicate the measurement in 
which that measurement noise component appears. All elements of the 
diagonal R matrix are zero because all nine measurements are assumed to 
be noise-free in this application. 

Consider the rms estimation error in Table 1 6 . Only two of the 12 rms 
state estimation errors are non-zero. Those two values are for estimates 
of the longitudinal and normal gust velocities. We can compare these rms 
errors to the input rms gust velocities; 

Vg = - 175 %Ug ~ 8 -^°u g 

°Ew g = ^Wg 3-2?= °w g 


°Ug = 2.06 

tfwg = 1-71 
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TABLE 16. UH-1H, HOVER, FINAL FILTER DESIGN 

Q MATRIX DIAGONAL, FILTER R MATRIX DIAGONAL, FILTER RMS STATE EST ERROR, 
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Because the estimation error is a small percentage of the rms gust velo- 
cities, the sensor array is judged to be adequate. In other applications, 
it may be desirable to compute the rms estimation error as reflected in 
the output vector, y, in order to verify that design performance target 
values are not exceeded. 

Consider the filter gain matrices K-] -\ and K-|2* -j does not appear 
in Table 1 6. This is the case because K-j ^ has zero columns when all meas- 
urements are noise free. The matrix in Table 1 6 (and other plant data) 
is used in the structure shown in Fig. 13 to define the filter -observer 
design. 

Next, we focus attention on the filter eigenvalues, in the center of 
Table l6. Previously, we identified a rather arbitrary range of 0.1 to 
50.0 rad/sec for "acceptable" eigenvalues. Clearly, the eigenvalues from 
this example fall within the acceptable range. Concerns with eigenvalues 
on the high end of the range centered around the high gains which might 
appear in the filter as a result. From the K-|2 filter gain matrix we see 
that all gains are at moderate levels, thus corroborating the acceptability 
of the eigenvalue at 15.6 rad/sec. 

On the other hand, the eigenvalues at the low end of the range were 
of concern because, in general, the filter eigenvalues become closed-loop 
system roots. Closed- loop roots with modulus less than 0.1 rad/sec clearly 
are undesirable. For example, they produce long tails in some of the tran- 
sient response characteristics of the system. In fact, even the eigen- 
values at 0.315 rad/sec and 0.333 rad/sec, for the example, could cause 
unacceptable transient response tails. But we shall see presently that, 
for this example and many other systems, these eigenvalues are not an 
important concern. 

This example constitutes a very interesting special case which has 
practical value. Notice that the filter -observer structure shown in Fig. 13 
implies that the controller and closed-loop filter eigenvalues can differ 
only if [G-^ — K12G2] no ^ nu H- * However, the data for and G2 

*This is the case because all filter modes are uncontrollable with 
respect to all elements of the plant control vector. 
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in Tables 1 6 and 17 indicate that [G-j — is null. Therefore, the 

controller and closed-loop filter eigenvalues are the same, and the closed- 
loop filter eigenvalues are independent of the regulator solution. Fur- 
thermore, in this special case the controller eigenvalues are therefore 
guaranteed to be stable. Stability of the controller-alone is usu ally 
desired in traditional flight control system synthesis. It is therefore 
important for us to recognize the causal factors which lead to similar 
optimal controller designs. 

One factor is the application of control through serially coupled 
actuation subsystems. Another factor is noise-free measurement of the 
actuator outputs. Together, these result in G-j being null and G2 being 
sparse. Additionally, we observe that the only non-zero elements of G2 
occur in the actuator state equations. Finally, we have observed in numer- 
ous applications that the coefficients of the actuator output measurements 
in K-|2 are zero as long as no noise-free measurement (other than actuator 
output itself) depends explicitly upon the actuator output states. K12G2 
is null when these two observations apply. We have already observed that 
when [G-| — a nul1 matrix, the filter eigenvalues and the con- 

troller eigenvalues are identical. 

The filter eigenvalues are, in general, closed-loop system poles. They 
do not, however, necessarily affect all aspects of closed-loop behavior. 

For instance, those inputs which drive the plant controls directly through 
a gain matrix result in closed-loop transfer functions in which zeros can- 
cel the closed-loop filter -observer poles exactly. (These inputs are those 
commands wherein error and integral error are noise-free, and have gains 
in the Bp matrix equal to zero.) This reduces, to some extent, the impor- 
tance of the filter eigenvalue location, since the modes they represent are 
not present in some of the command transfer functions of the closed-loop 
system. There remain, however, three areas where the filter eigenvalues 
do play a role : 

• The filter eigenvalues become roots of the controller 
transfer functions. If the controller is used apart 
from the total FCS (e.g., in a checkout phase) the 
eigenvalue locations will partially determine the 
controller's response characteristics. 
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TAB IE 17. CONTROL DISTRIBUTION MATRIX PARTITIONS 
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• The poles of the closed-loop system identified as 
filter eigenvalues will not he cancelled exactly in 
general by closed-loop system zeros for transfer 
functions where the input is a disturbance. Often, 
however, they nearly cancel with zero, and are thus 
not much of a concern. 

• In the next chapter we will consider the behavior 
of the closed-loop system in the face of various 
nonlinearities and component degradation. These 
off-nominal conditions have the effect of opening 

or changing one or more of the feedback loops. This 
will naturally destroy the cancellation of filter 
eigenvalue poles by closed-loop zeros. Considera- 
tion of filter eigenvalue location then becomes quite 
important . 

Turn attention now to the matter of augmenting process noise excita- 
tion of the plant in such a way that differentiation of any noise-free 
measurement is not required. Consider the process noise distribution 
matrix, r, in Table l8. The upper and lower partitions of r are Tq and 
T c , respectively. The left and right partitions of p contain the coeffi- 
cients of the process noise components arising from the problem physics 
and arising from augmentation, respectively. Notice that the upper right 
partition of r is null in accordance with our recommendation for augment- 
ing Tg, and the lower right partition is identity in accordance with our 
recommendation for augmenting Tg. 

Consider the r 2 matrix for this example application In Table 19 . The 
first three columns of r 2 distribute those process noise components which 
have a basis in the problem physics into the derivatives of the noise-free 
measurements. The first two columns cannot contribute to the rank of rp* 
The last nine columns of r 2 distribute the augmenting process noise com- 
ponents. Column 11 cannot contribute to the rank of T2i therefore, the 
EXB component of the process noise intensity can be zero. (This is the 
eleventh element in the diagonal Q matrix.) The remaining nine columns 
of r 2 are linearly independent,- therefore, r 2 has full rank. Eight of 
these nine columns are associated with augmenting process noise compo- 
nents. The levels for these eight components must be greater than zero 
in the diagonal Q matrix. 
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TABLE 18. PROCESS NOISE DISTRIBUTION MATRIX 



000 0.000 0.000 I 0.000 0.000 0.000 0.000 0.000 0.000 1.00 0.000 0.000 l X12 



TAB HE 19- PROCESS NOISE DISTRIBUTION MATRIX FOR DERIVATIVES 
OF NOISE-FREE MEASUREMENTS 
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The final filter -observer design summarized in Table 1 6 is the 
result of three filter -observer design trials. The initial filter- 
observer design is summarized in Table 20. This initial design differs 
from the final only with respect to levels used for the augmenting pro- 
cess noise. Thus, Table 20 can be compared with Table 1 6 to demonstrate 
the substantial influence adjustment of augmenting process noise levels 
may have upon closed-loop filter eigenvalues. Initial design levels are 
larger for most components, and especially so for the u state equation. 

The initial design has a range of eigenvalues which is unacceptable 
because the smallest eigenvalue modulus is less than 0.1 rad/sec. This 
design also is characterized by somewhat larger rms estimation error and 
larger gains in the matrix:. Adjustment of the augmenting process 
noise levels (by trial and error) brings the low modulus eigenvalue well 
into the acceptable range (refer to Table l6). 

Another interesting and practical dimension to filter-observer syn- 
thesis involves use of different sensor complements. The study described 
below explores the results of incorporating measurements of longitudinal 
acceleration independent of pitch and normal acceleration. Application 
in this case is for longitudinal dynamics of the UH-1H at a 100 kt level 
flight condition. The baseline sensor complex includes noise-free meas- 
urement of instantaneous vertical speed, pitch rate, pitch attitude, 
longitudinal cyclic actuator output, main rotor collective actuator out- 
put, rate -of- climb error, integral rate-of-climb error, airspeed error, 
and integral airspeed error. Five additional cases are considered in 
terms of measurements added to the baseline measurement set. The char- 
acteristics of these cases are summarized in Table 21. 

Table 22 summarizes the design parameters, closed-loop filter eigen- 
values, rms estimation error performance and the filter gain matrix for 
the baseline case of the UH-1H at 100 kt. In the process noise inten- 
sity diagonal matrix, Q, the fourth and subsequent elements are augmenting 
process noise components. The intensities of these augmenting components 
are kept constant in this study to the extent that the addition of noise - 
free measurements permits. The baseline filter design is third order,* 
its eigenvalue range is acceptable, and the rms state estimation error 
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TABLE 21. SENSOR COMPLEMENT STUDY 
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a ICRB denotes the instantaneous center of rotation location for longitudinal 
cyclic pitch control inputs, Zgg/Mgg. 


















TABLE 22. BASELINE FILTER DESIGN, UH-1H, 100 KT 

ILTEN ItMS STATE EST ERROR, FILTER 
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is acceptable although estimates of longitudinal gust and inertial velo- 
city are not very accurate. The controller and closed-loop eigenvalues 
are the same because the K 12 gain matrix coefficients of DB and DC are 
zero. (Gi and G 2 are the same as in Table 1 7* ) 

In Table 23 we explore the effect of adding a noise-free normal 
acceleration measurement. The principal effects are to lower the filter 
order from third to second order, and to produce non-zero coefficients 
for DB and DC in the K-j2 filter gain matrix. This latter change means 
that controller eigenvalues will be a function of the regulator gains, 
and the controller eigenvalues will not necessarily be stable. We assert 
that non-zero coefficients of DB and DC in K-\2 arise because the noise- 
free measurement AZP is an explicit function of DB and DC in the measure- 
ment equations. The effect of the additional measurement in reducing rms 
state estimation error is almost nil except for the normal gust state 
which is halved. 

In Table 24 we demonstrate the truth of the above assertion. Here 
we relocate the noise -free measurement of normal acceleration, AZP, from 
the pilot’s station to instantaneous center of rotation for longitudinal 
cyclic pitch control inputs.* The latter location results in elimina- 
tion of the explicit dependency of the AZP measurement equation upon DB 
since its coefficient (Zg-g — M§-g) ^eing zero defines for the 

instantaneous center location. The principal change with respect to 
the filter design of Table 23 is that the coefficients of DB are again 
zero in the K 12 filter gain matrix. 


*This location is frequently regarded as a desirable one for normal 
acceleration measurement in fixed-wing aircraft. This is so because the 
accelerometer measurement is then very nearly proportional to angle of 
attack which, in turn, is ideal for augmentation of short-period dynamic 
characteristics (Ref. 10, pp. 446-453). 
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Data in Table 25 show the effect of using a noisy measurement of 
normal acceleration. The differences by comparison with Table 23 are 
that the filter order is increased from second to third order, and the 
coefficients of DB and DC in the K 12 filter gain matrix are again zero. 
Consequently, the controller eigenvalues are stable and are the same as 
the closed-loop filter eigenvalues. Notice that for the first time the 
K-| 1 filter gain matrix has a number of columns greater than zero. The 
rms state estimation error is not much affected by the normal accelera- 
tion measurement noise. 

Table 26 explores the effect of adding noise-free measurement of 
longitudinal-acceleration-independent-of-pitch, AXP, to the baseline 
measurement set. It is interesting to explore using this additional 
quantity because of its traditional use in flight control technology to 
provide smoothed, broadband airspeed feedback signals via complementary 
filtering. The principal effects are to again reduce the filter order 
from third to second order, and to produce non-zero coefficients of DB 
and DC in the K-jp filter gain matrix. Somewhat to our surprise, there 
is a small reduction in rms estimation error for normal gust velocity 
and almost no reduction for all other states. 

In Table 27 we explore the effects of adding noise-free measurements 
of both normal acceleration and longitudinal-acceleration- independent-of- 
pitch 'to the baseline measurement set. The principal effect is that the 
filter is first order. Coefficients of DB and DC in the K-jp filter gain 
matrix are non-zero. The rms state estimation error is not changed 
significantly with respect to values obtained when either acceleration 
measurement was used separately (refer to Tables 23 and 2 6 ). 

7- Re -Emphasis of Key Points in 

Filter -Ob server Synthesis 

To conclude the design synthesis of the filter-observer, we empha- 
size again the three ways the filter-observer affects the final controller 
design: 
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TABLE 25. EFFECT OF a. z MEASUREMENT NOISE ON FILTER DESIGN, UH-1H, 100 KT 

Q MATRIX DIAGONAL, FILTER R MATRIX DIAGONAL, FILTER R^'IS STATE EGT ERROR, FILTE 
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TABLE 27. EFFECT OF NOISE -FREE (a x - g sin ©) AND MEASUREMENTS ON FILTER DESIGN 


U 2 



O 

o 




r— ze cr. cn 


vo 

r. 

O' 


o 

I 






93 




K12 GAIN MATRIX, FILTER 



• The filter -oh server synthesis process provides a 
means for choosing the best sensor complex from 
among many candidates. As a result, the quality 
of the complex is also quantified. 

• The filter -observer provides the best estimates 

of the total state vector, given the desired meas- 
urements, to use as the regulator input. 

• The filter-observer provides the only direct means 
for the noise and disturbances to affect closed- 
loop system behavior. If the matrix [G-j — KigGg] 
is null, then the controller and closed-loop filter 
poles are identical. The controller poles are then 
also guaranteed to be stable. 

The above considerations lead to the recommendation that filter- 
observer synthesis precede regulator synthesis. Because these are com- 
pletely separable design processes, the filter-observer design can be 
set aside once satisfactory performance is achieved, and work then pro- 
ceeds with the regulator. This is the topic of the next subsection. 

C. REGULATOR SYimiESIS 

The steps involved in synthesis of the regulator portion of the con- 
troller are summarized in Table 28. Although ten steps are listed. 

Steps 3 through 8 are repeated until all bandwidth design goals have been 
addressed in any given application. These ten steps are elaborated in 
the following paragraphs where required, and are illustrated by examples. 

1 . Coat Function Control Weightings 


One element of the cost function diagonal control weighting matrix, 
Rr, should be set initially to unity. This is permissible since one 
coefficient in the cost function may be an arbitrary positive number. 

The remainder of the diagonal elements are chosen in ratios to the first, 
such that limitations of the physical controls, or upon their use, are 
taken into account. These limitations may be the result of rate, author- 
ity, or power limits or possibly the result of flying qualities con- 
straints, e.g., avoid "throttle thrash." Specifically, the remainder 
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TABLE 28. 


REGULATOR SYItfTHESIS 


1 . Set R R ta accommodate limitations upon controls or their 
use. 


2. Order key variable bandwidth requirements in decreasing 
order . 


5- Identify highest bandwidth requirement not yet addressed, 
<%• 

4. Compute transfer functions 

Output; Variable subject to requirement 

Inputs: Variables satisfying prior requirements 

which may serve as surrogate controls 

All loops satisfying prior requirements 

must be closed 


5- Approximate above transfer functions at s = jo^ by 
K/s n+1 . 

6. Calculate (2 n o^) n+ ^ r/K 2 for each approximation (r is 
cost function weighting on input variable). 

7. Set smallest (2 n ac§) n+1 r/K 2 value equal to cost function 
weighting, q R , for variable subject to cjd c requirement. 

8. Compute regulator solution with additional cost function 
weighting, q R . 

9* Repeat from Step j5 until all bandwidth requirements are 
addressed. 


10. Change R R to affect relative use of the available control 
points,- maintain ratios of 0 uR and R R elements supporting 
related control objectives. 
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of the diagonal elements is the ratio of the square of the limit on the 
element having the unity coefficient to the square of the limit on the 
control variable being weighted. 

A particularly desirable alternative is to use the R R values result- 
ing from a previously completed solution at a nearby flight condition. 

The Rr values used at this stage of design are ultimately modified in 
the last step. 

2. Order Key Vaxiable Bandwidth Requirements 

Designation of key variables and their associated control bandwidth 
requirements define the primary design goals for regulator synthesis as 
explained in Section II. The ordering of these key variables and their 
control bandwidth requirements according to decreasing bandwidth provides 
a rational basis for systematic synthesis. This approach results in 
satisfaction of the more stringent, high gain, inner loop, high band- 
width requirements firsts the next most stringent, and so on. Proceed- 
ing in this manner assures that the more stringent of either all prior 
bandwidth requirements or the minimum system augmentation necessary to 
support the immediate bandwidth requirement is met at each stage of 
design. It is also interesting to note the similarity of this procedure 
to the one described on page 664 of Ref. 10 for synthesis via frequency 
response methods. 

3. Computation of Transfer Functions 

Each bandwidth requirement is addressed in turn, starting with the 
highest. Consider the highest bandwidth requirement not yet addressed. 
(We shall designate this bandwidth by since it is approximately equal 
to crossover frequency, as demonstrated in Section II.) Transfer func- 
tions are computed for the variable which is subject to the bandwidth 
requirement. The input variables for these transfer functions are 
selected by the designer from among the controls and all variables 
satisfying prior bandwidth requirements. (The latter variables may be 
thought of as surrogate controls.) If the designer is unfamiliar with 
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the conventional bases for selection (e.g., airspeed is controlled with 
pitch attitude in cruising flight ), all of the above named variables may- 
be included. Later steps will result in an optimal selection. All loops 
satisfying prior bandwidth requirements must be closed when computing 
these transfer functions. 

Approximation of Transfer Functions 

Approximations to the computed transfer functions are required for 
the frequency region near o^. The form of the approximations required 
is K/s n+ ^. There are two viable ways for obtaining these approximations. 
The first, and easiest, method is to simply approximate the transfer 
function by its magnitude asymptote at s = jcn c . The appropriate phase 
approximation is carried in the sign of the gain for the approximation 
and the value of n + 1 . The second method involves finding the tangent 
to the Bode plot at s = jcd c . The phase approximation is again carried in 
the sign of the gain for the approximation and the value of n + 1 . Either 
method results in a pair of values for K and n for each transfer function. 
The values are somewhat dependent upon the approximation method used, but 
the effects of this dependency are minimized by the man ner in which K 
and n are ultimately used. It is our recommendation that one approxima- 
tion method be used throughout the regulator synthesis, however. 

5. Obtaining the Cost Function Coefficient 

The cost function coefficient, q R , for the variable subject to the 
current bandwidth requirement is determined next. (q R denotes an element 
of the diagonal matrix, ) The quantity 

(2 n a|) n+1 r/K 2 

is calculated using the K and n values for each transfer function approxi- 
mation. r in the above quantity represents the cost function weighting 
coefficient for the control or surrogate control variable corresponding 



to the particular transfer function. This may he an element of the 
diagonal R R matrix or an element of Qp determined in a previous step. 

q^ for the current bandwidth requirement is set equal to the small- 
est calculated value for the above quantity. 

This step contains the key to understanding of this approach to regu- 
lator synthesis. The quantity being calculated is an evaluation of the 
effectiveness of each control or surrogate control variable for achiev- 
ing the specified bandwidth. Large values for q R indicate low effective- 
ness and vice versa. The lowest value for q R is picked so as to insure 
that the most effective control (or surrogate control) variable is used 
to achieve the specified bandwidth. 

6. Continuing the Regulator Synthesis 

The next step is computation of the regulator solution using the 
additional cost function weighting coefficient, q R . At this point the 
designer is ready to address the next bandwidth requirement and so returns 
once again to the step for computation of the appropriate transfer func- 
tions. The sequence of steps from that point repeats (at least in prin- 
ciple) until all bandwidth requirements have been addressed. 

As experience is gained with the controlled element and the series 
of regulator design steps, the designer may find that several, or all, 
bandwidth requirements can be addressed simultaneously. The reasons for 
this are that uncoupled or weakly coupled or serial subsystems within the 
controlled element often present situations wherein the effects of loop 
closures satisfying prior bandwidth requirements upon the various trans- 
fer functions are known or easily estimated without recalculating trans- 
fer functions. This is especially the case when the effect of one (or 
more) prior loop closures is nil, or if one (or more) prior loop closures 
is at very high gain. In the latter case, ratios of transfer function 
numerators of higher kinds give close approximations to the actual closed- 
loop transfer functions (refer to Ref. 10 , pp. 163 - 177 , and Ref. 11 , 
pp. 63 - 70 ). 


98 



Techniques using numerators of higher kinds (that is to say, multi- 
variable system transmission zeros) are the core of a regulator synthe- 
sis technique recently advanced by Harvey and Stein (Ref. b -) . Their 
technique is based upon a single parameter relaxation from a singular 
regulator problem solution. The particular singular regulator solution 
selected for the relaxation starting point is derived from the modal 
properties desired for the finite modulus poles of the singular system. 

7 . Adjustment of the Diagonal Rr Matrix 

After all bandwidth requirements have been addressed, it may turn out 
that the relative use of control points is inappropriate. For example, 
the main rotor collective pitch control may be more active in terms of 
zero crossings per second than the longitudinal cyclic pitch control in 
level cruising flight. (Main rotor collective should be a trim control 
in this situation.) Inappropriate control use (required authority, 
required rate or activity too high) is rectified by increasing the cost 
function weighting coefficients in Rp for those controls. In order that 
bandwidth requirements continue to be satisfied, it is necessary to 
increase by the same factor all those elements in Qp which have a ratio 
relationship to each element of Rp being increased. (This ratio rela- 
tionship arises from successive applications of qp = min [(2 n o)§) n ' f "1 r/K].) 
Adjustments made in this manner assure simultaneous satisfaction of band- 
width and control use requirements. 

8. Regulator Synthesis Examples 

This subsection illustrates application of the ten regulator design 
steps in Table 28. Application is to longitudinal dynamics of the UH-1H 
for the hover flight condition. Regulator design function and bandwidth 
requirements are listed in Table 29* 

Step 1 ■ The control weighting cost function coefficient arbitrarily 
selected to be unity is that weighting longitudinal cyclic pitch rate, 
DBD. 

R R (1) = 1.0 
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TABLE 29 


UH-1H LONGITUDINAL FLIGHT SYSTEM FUNCTION 
AND BANDWIDTH REQUIREMENTS' 


Flight Condition: 


Hover 


Functions : 


Rate-of-climb command 


Groundspeed hold 


Bandwidths : 

(Rad/Sec 

Cyclic (DB) 

25 . 7 

Collective (DC) 

25.7 

Pitch (TH) 

2.0 

Rate-of-climb error (HDE) 

1 .0 

Integral of HDE (HDl) 

0.82 

Groundspeed error (XD) 

0.5 

Integral of XD (XDl) 

0.1 


Our expressed desire is to have maximum collective pitch rate, DCD, he 
significantly less than maximum longitudinal cyclic pitch rate to satisfy 
flying qualities requirements. A ratio of ^To to 1 is selected on the 
basis of flying qualities engineering judgment. Consequently, the cost 
function coefficient weighting collective pitch rate is: 



Step 2 . The key variable bandwidth requirements are already in 
order of decreasing bandwidth in Table 29 . 
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Step 3 . Refer to Table 29* The highest bandwidth requirement not 
yet addressed is for the longitudinal cyclic pitch actuator output, DB. 

(The bandwidth requirement for the collective pitch actuator output, DC, 
is equal to that for DB. When this occurs it is immaterial which require- 
ment is addressed first. Indeed, since the UH-1H has no unstable roots, 
and the two actuators are completely independent, both requirements can 
be addressed simultaneously. ) 

Step 4 . Compute transfer functions. The output variable subject to 
the bandwidth requirement of 25-7 rad/sec is DB. The input variables 
which may serve as controls are DBD and DCD. (There is no possibility 
for surrogate controls because no prior bandwidth requirements have been 
satisfied. Consequently, there are no prior loop closures to incorporate.) 
The candidate transfer functions are; 

DB = J__ _DB_ _ Q 

DBD s DCD 

That this is so is readily determined from inspection of the plant equa- 
tions (Table 10). Similar considerations lead to transfer functions for 
collective pitch, DC. 


DC = J_ JC_ n DC 

DCD s DCD DB 

The third transfer function for collective pitch arises from the poss- 
bility that DC, which "satisfies a prior bandwidth requirement, " may be 
considered a surrogate control for DB. The uncoupled nature of the two 
actuators and stability of the open-loop plant result in DB having no 
effect upon DC even when the loops satisfying the bandwidth requirement 
for DB are closed. 

Step 5 . Approximate the transfer functions. No approximation is 
required in this case since the transfer functions are already in the 
form K/s n+1 . The values for K and n are; 
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Transfer Function K, n 

5 B /6 B (DB/DBD) 1.0 0 

5 C /5 C (DC/DCD) 1.0 0 

The zero valued transfer functions need not "be considered. 

Step 6 . Calculation of (2 n o ^) n+1 r/K 2 . For the DB requirement the 
value is: 

[2° (25- 7 ) 2 ] 1 R E (1)/(1.0 ) 2 = 662. 

For the DC requirement the value is; 

[2° (25- 7 ) 2 ] 1 R r (2)/(1.0 ) 2 = 6620. 

Step 7 . Selection of the cost function coefficient. There is only 
. / n 2 \ 1 / p 

one finite value of (2 ouc) r/K for the DB bandwidth requirement. It 
is therefore the smallest value and is set equal to the corresponding 
cost function weighting for DB, 0^(4). 

Qr(4) = 662 . 

The situation for the DC bandwidth requirement is similar. The cost 
function weighting for DC is Q^( 5 ). 

Qr( 5 ) = 6620 . 

Step 8 . Compute the regulator solution with the additional cost 
weighting. The computed regulator gain matrix, K R , for Qr(^) = 662 . and 
%(5) = 6620. is a 2 X 12 matrix of zeros except for two elements: 

K r ( 1,4) = 25.7 rad/sec 

Kr(2,5) = 25.7 rad/sec 
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Step 9 . Repeat the design process until all bandwidth requirements 
are addressed. All requirements are not addressed^ therefore return to 
Step 3 * 

Step 3 « Refer to Table 29 * The highest bandwidth requirement not 
yet addressed is for pitch attitude, = 2.0 rad/sec. 

Step !<- . Compute transfer functions. The output variable subject to 
the bandwidth requirement is pitch attitude, TH. The candidate controls 
and surrogate controls are DBD, BCD, DB and DC. We will restrict our 
attention to the surrogate controls DB and DC for the purpose of illus- 
tration. These transfer functions may be obtained from the appropriate 
ratios of numerators for the open-loop controlled element transfer func- 
tions in Table 30 . (A flight controls engineer would ordinarily reject 
collective pitch, DC, as a candidate for controlling pitch attitude. We 
shall continue to include it here in order to show that the design pro- 
cess automatically rejects use of collective pitch for pitch attitude 
control. ) These transfer functions and the corresponding Bode plots are 
given in Figs. l 4 a and iVb. (The two actuator loops satisfying the prior 
bandwidth requirements are closed, but this does not happen to affect 
the transfer functions in this particular case.) 

Step 3 - Approximate the transfer functions. Two methods for 
approximating transfer functions are illustrated in Fig. l6a for the 
TH/DB transfer function. The solid amplitude ratio curve in Fig. l6a 
is | G- ( j 00) and the dashed curve is composed of the asymptotes for 

|G( joo) | (jg. In the vicinity of the pitch attitude bandwidth requirement, 
u^TH = 2.0 rad/sec, the asymptote is given by -0.169/s 2 . (This can also 
be computed readily from the factored transfer function itself. It is 
accomplished by making low-frequency approximations to the factors with 
root modulus greater than 2.0 rad/sec, and high-frequency approximations 
to the factors with root modulus less than 2.0 rad/sec. Multiplying 
out the result gives an answer directly in the form K/s n+1 ) . 

Approximation of the actual amplitude ratio curve requires the Bode 
plot of Fig. 16 a to expedite approximation. A tangent to the curve is 
drawn at m = 2.0 rad/sec. JkI^ is determined by the intersection 
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TABLE 30. OPEN-LOOP CONTROLLED ELEMENT TRANSFER FUNCTIONS 


CASE: UH1H HOVER 122LONG 31-JAN-79 CONTROLLED ELEMENT TF'S 


DENOMINATOR: 

1.0000 




( .00000 ) 

( .38494 ) 

( .00000 

) ( .00000 ) < 

. 00000 

(( .10339 

.18934 

, . 19576E-01 , .18833 

)) 

(( .26279 , 

< .11863E-01> 

.92717 

, .24365 , .89459 

)) 


NUMERATOR: TH/DBD 


16910 

( .00000 )' 

( .33300 ) 

< .17445E-03> 


( .00000 ) 

( .39184 ) 


{ .00000 ) 


(-. 79065E-02) 


NUMERATOR: TH/DCD 

-. 33000E-02 

( .00000 ) ( .00000 ) ( .00000 ) ( . 10854E-01 ) 

( .33300 ) (-11.270 ) 

< . 13442E-03> 


NUMERATOR: HDE/DBD 


.24729 
( .00000 
(-1.4789 
(( .35965 
<-.22334 


) ( .00000 
) 

, 1.3542 

> 


) ( .00000 ) ( .33300 ) 

.48705 , 1.2636 )) 


NUMERATOR: HDE/DCD 


-9.7980 
( .00000 
(( .18545 
(( .25760 
<-. 17544 


) ( .00000 
, .14356 
, .93210 


) ( .00000 ) 

. 2662 3E-01 , .14107 )) 

.24011 , .90064 )) 


NUMERATOR: HDI/DCD 


-9.7980 
( .00000 
(( .18545 
(( .25760 
<-.17544 


) ( .00000 
, .14356 
, .93210 


) 

. 26623E-01 , ,14107 )) 

.24011 , .90064 )) 


NUMERATOR: XD/DBD 

1. 0606 

( .00000 ) ( .00000 ) ( .00000 ) ( .33300 ) 

( .38296 ) 

(( . 19712E-01 , 2.2663 , .44674E-01, 2.2659 )> 

< .69473 > 


NUMERATOR: XDI/DBD 

1.0606 

( .00000 ) ( .00000 ) ( .33300 ) ( .38296 ) 

(( .19712E-01, 2.2663 , .44674E-01, 2.2659 )) 

< .69473 > 



K/s n+1 

APPROXIMATION 

( 2 n «£ IH ) n+1 W (® 2 

Asymptotic 
K = -0.169, n = 1 .0 

1,484,183 

G(Jco) 

K = -0.316, n = 1.5 

2,851,632 


1691 (s - .0079)(s + . 333)(s + .39184 


(s+ .38k9k-)[s 2 + 2( . 1 6339) ( . i8934)s + ( . 1893U) 2 ] 

x [s 2 + 2(. 26279) (.9271 7) s+ (.9271T) 2 ] 


& ( rad /sec ) 


Figure 16a. Transfer Functions for Selecting Qr( 3) Value for 
Pitch Attitude Control Bandwidth Requirement 
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Phase (deg) Amplitude (dB) 


K/s n+ 1 

APFROXIMA.TIOH 


QdcAk® 


Asymptotic 
K = 0.0372, 


n = 2.0 


19, 604,^68, 000 


2 


_e_ 

5 C 


.01 

20 - 


--.0033(s + . 

(s + .38494) [s 2 +2( 


0108^4) (s + .333)0 + 11.27) 
• 1 6339) (-1 8934) s + (.1893*0 


2 ] 


x [s +2(. 26279) (• 9271 7)s + (.92717) 2 ] 


.10 


uj ( rad /sec) 


1.0 



Figure 1 6b. Transfer Functions for Selecting Qr( 3) Value for 
Pitch Attitude Control Bandwidth Requirement 
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of the tangent with co = 1.0 rad/sec, and — (n+ l) is the slope of the tan- 
gent in 20' s of dB's per decade. 

Either method of approximation may be used, but one method should be 
used exclusively throughout. We shall continue using the asymptotic 
approximation which, for the TH/DC transfer function, is 0. 0372/s ^ as 
shown in Fig. l6b. The asymptotic approximations are summarized below. 

Transfer Function K n 

e/8 B (th/db) -0.169 1 

0/S c (TH/DC) 0.0372 2 

Step 6 . Calculation of (2 n ae§) n+1 r/K 2 . 

For the DB control: 

(2 n o| TH ) n+1 Q R (1)/(K™) 2 = [2 1 (2. ) 2 ] 2 662. /(-0. 169) 2 

= 1,484, 183. 

For the DC control; 

(2 n a| TH ) n+1 Q r (5)/(K^) 2 = [2 2 (2. ) 2 ] 3 6620./(0.0372)2 

= 19,604,^68,000. 

Step 7* Selection of the cost function coefficient. The smallest 
value of (2 n a|) n+1 r/K 2 is set equal to the cost function weighting for 
pitch attitude, Q R (3). 

Q r (3) = 1,484,183. 

From here, all steps are very similar up to the point where all band- 
width requirements have been addressed. The regulator solution that 
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results is for the "completed" design. Table 31 summarizes this regula- 
tor solution. The cost function weightings, closed-loop eigenvalues and 
regulator gain matrix are given in this table. 

We are now in a position to consider how Rp might be modified in 
Step 10. This requires knowledge of the relationships among the control 
and output variables that were operative in determining the cost function 
coefficients. These relationships are summarized in Table 32. This 
table has the following interpretation. Consider either of the last two 
columns, DBD or DCD. The bottom-most Qp entry indicates the outer-loop 
variable which is the ultimate object of control. (The variable name is 
read in the left-most column. ) The Qp entry immediately above indicates 
the control (or surrogate control) variable for the column entry below. 
This interpretation applies on up the column to each Qp or Rp entry in 
turn. It also applies (individually) for each column. Thus, fxdt has 
as a surrogate control x which, in turn, has the surrogate control 0. 
Similarly, / he dt is preceded in the control chain by h e , which is con- 
trolled by 5 C , which is effected finally by 5 C , the "actual" control. 

The "Control Variable Association" differentiation directly reflects 
the origin of the smaller Qp cost weighting at each bandwidth step. Thus, 
Qp(3), which Is the weighting for pitch attitude, 0, derived from the 
1.48 x 10^ found for DB as opposed to the 1 .96 X 10 10 for DC. 

If every entry in any one column is multiplied by a common factor, 
the design bandwidths will tend to be unchanged. (Exceptions will occur 
only if the multiplying factor is such that one or more Op's would change 
columns if the formal procedures of Steps 3 through 9 were followed. ) 

Because the longitudinal speed and pitch degrees of freedom are 
nearly uncoupled from the vertical speed degree of freedom, it is not 
terribly interesting to illustrate the impact of multiplying every entry 
in one column by a constant factor for a hover flight condition and vice 
versa (but to a lesser extent). On the other hand, the 100 kt cruise 
condition is interesting because longitudinal speed, pitch and vertical 
speed couplings are strong. The relationships among control and output 
variables and cost function coefficients for this 100 kt flight condi- 
tion are summarized in Table 33. The cost function coefficients related 
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TABLE 32 


RELATIONSHIPS AMONG CONTROL AM) OUTPUT VARIABLES 


AND COST FUNCTION COEFFICIENTS — HOVER 

Output 

Variable 

Bandwidth 

Requirement 

(rad/sec) 

Control Variable 

Association 

DBD 

DCD 



%(i) = i-oo 

Rr( 2 ) = 10.0 

5 B (DB) 

25-7 

Qr( ^ ) = 662. 


8 C (DC) 

25- 7 


Q r(5) = 6620. 

e (th) 

2.0 

Qr( 3 ) = 1 A8x 10 6 


h e (HDE) 

1.0 


Qr(6) = 69.0 

/hedt (HDI) 

0.82 


%(7) = ^6.6 

x (XD) 

0.50 

Qr(8) = 357 . 


/ x dt (XDl) 

0.10 

Qr( 9 ) = 3.57 



TABLE 33 

RELATIONSHIPS AMONG CONTROL AND OUTPUT VARIABLES 
AND COST FUNCTION COEFFICIENTS — 100 KT 


Bandwidth 

Control Variable 

Association 

Output 

Requirement 




Variable 

(rad/sec) 

DBD 


DCD 



Rr( 1 ) = 1.00 

Br(2) 

11 

b 

0 

i-b 

5b (DB) 

25.7 

Qr(^) = 662 . 



5 C (DC) 

25.7 


%(5) 

= 662 . f 

0 (TH) 

1-5 

Qr(3) = 3.55X 103 



he (HDE) 

0.70 


Or ( 8 ) 

= 5-95 f 

f h e dt 

0.32 


Qr(6) 

= k.77f 


f = 1 . 

, 00 , 10 . 0 , 100 . 
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to height control (the right-most column of Table 35) are multiplied by 
f = 1 .00, 10.0 and 100. 

Summaries of the resulting regulator designs are given for Tables 3^, 
35 > and 3 6 for f = 1.00, 10.0 and 100., respectively. These summaries 
include the cost function coefficients, closed-loop regulator eigen- 
values, regulator gain matrix, and rms responses for states, outputs and 
controls. The gains based on control deflections for the three values 
of f are summarized as : 


2TTT + 1 


-.012 


-.013 

1 

5 b - 

o 

o 

vji 

Ug + 

-. 108 

5T7T + 1 


.022 


“3Ti 

o 

OJ 

'Ll 


-.069 


.008 


-.053 


” 

9-63 

-■315 

h c + 

.009 

u + 

-.206 

w + 

9.86 

-.923 


.026 


— 727 


7.^9 


m 


28.3 


.018 


-.031 


-.718 

62.6 

9 + 

.093 

5 C + 

-.181 

/ h e dt + 

-.867 

l6l.4 


.334 


-.66 6 


-.965 

L J 


s 

20 + 1 


-.008 


.053 

2^5 + 1 

5 C ~ 

-.004 

Ug + 

.026 

2^T + 1 


-.001 


.006 


- 




— 



.132 


-.001 


.079 


-.664 

.082 

he + 

-.001 

u + 

.057 

V + 

.01 

.034 


-.001 


.028 


.139 

__ 


_ 


_ - 


_ 



-17.9 


.019 


.075 


.163 

H 

-12.2 

6 + 

.01 

5b + 

.054 

/ h e dt + 

.049 


- 5.44 


.004 


.026 


.01 


f u AS e 


dt 


(W 


f = i . 
f = 10. 


, f = 100. 




1 12 


218 -1.41 -3.55 0.296S-01 -2.11 17.8 47S. -0.501 
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REGULATOR GAIN MATRIX 



Recall that the equation relating the controls to the states for the regu- 
lator is u = — Cx. In Eqs. 47 and 48 we have moved the state variables Sg 
and 5 q to the left-hand side of the equations for and Sq, respectively. 
This is because the state variables 63 and 6^ always remain integrals of 
the control variables S B and 

The gains across the three values for f, for several variables in each 
equation, are plotted vs. f in Fig. 17- The most noticeable effect, and 
one which is wholly expected, is that the gains in the 5g equation (Eq. 47) 
increase and the gains in the 5^ equation (Eq. 48) decrease as f increases. 
This is so because increasing f effectively increases the penalty on use of 
as a control point. Thus, the regulator is forced to rely more heavily 
on the longitudinal cyclic control point. 

Figure 17 also illustrates, for the variables plotted, that a sort of 
power law might be appropriate for the effect of increasing the 5c cost 
function weighting on the regulator gains. To the extent that such a 
relationship exists, it might be exploited in further refining the method 
for computing desired cost function weightings. 

The highest frequency eigenvalues representing the inner-loop band- 
widths ' associated with the actuators are changed only slightly as f is 
varied. This is to be expected since a relative change in emphasis 
between collective and cyclic controls should not affect the highest fre- 
quency modes. Similarly, the three lowest frequency eigenvalues, 0 . 100 , 
0.11 6, and O.269, related to the shaping filters Ug, Wg, and h c , respec- 
tively, do not vary with f at all. This too is to be expected, because 
the shaping filter modes are not controllable. They are used only to 
shape the white process noise. 

A comparative summary of rms response and control activity is given 
in Table 57* Just as with the gain shifts of Eqs. 47 and 48, this per- 
formance summary shows a very marked shift from use of collective pitch 
to control altitude to use of longitudinal cyclic pitch and pitch rate 
to control altitude as f is increased. (Consider the rms values for DCD 
and DC vis a vis the rms values for DBD, DB and Q for the three values 
of f . ) Yet outer-loop performance is not changed that significantly. 
(Consider the rms values for HDE and TH. ) As f is increased, one real 
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TABLE 37 


COMPARATIVE PERFORMANCE SUMMARY 

(For Closed-Loop System Including 
Filter Observer) 


Multiplying Factor, f 


Variable 

1.00 

10.0 

100 . 

h (HD) 

1 .01 

1.02 

1.01 

q (Q) 

0.238E-02 

0 . 758 E -02 

0. 141E-01 

0 (TH) 

0.373E-02 

0 . 899 E -02 

0. 119E-01 

5b (DB) 

0.491E-01 

0. 101 

0.208 

S C (DC) 

0. 122 

0.455E-01 

0.948E-02 

HDI (/ h e dt) 

0.385 

0.254 

0. 1 12 

az (AZP) 

0.519 

0.575 

0.875 

h e (HDE) 

0.263 

0.253 

0.175 

x (XD) 

1 .to 

3-11 

3.84 

a (AOA) 

0.886E-02 

0.669E-02 

0 . 766 E -02 

5 b (DBD) 

0.129 

0.757 

2.01 

5 C (BCD) 

0.331 

0. 181 

0 . 599 E -01 

DB control 
activity 8. 

0.4i8 

1 . 19 

1.54 

DC control 
activity 3 - 

0.432 

0.633 

1 .01 


x control activity = Qx/( 2ita x )* This is the expected 
number of positive -going zero crossings per unit time. 


root increases in modulus by a factor of 3*7 and the closed-loop short- 
period damping ratio decreases from 0 .56 to 0.28. This shift in relative 
use of controls is accompanied by an increase in control activity for 
both controls. As seen from the bottom of Table 37 , there is an almost 
fourfold increase in positive -going zero crossing rate for the 5$ control, 
as f is increased by two orders of magnitude. Notice also that this par- 
ticular- control activity measure also illustrates the dramatic shift in 
the control burden from Sq to 5 B as the penalty on Sq is increased in 
the cost function. 

Figure 18 shows the rms response of several variables and the two 
control activities plotted as functions of the three values of f. First, 
notice that the slope of the line depends on the control point associated 
with the particular variable^ the 5^ variables show a negative slope, 
while the 5p variables have positive slope. This same feature was ob- 
served in Fig. 17- We might also infer, as before, that some power law 
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exists which approximately defines the effect of increasing the Sq con- 
trol point weighting upon the various rms response levels. Since such a 
relationship does not hold for variables not plotted in Fig. l8, we can- 
not draw any conclusions here. 

Flying qualities dictate that collective pitch be used as a trim 
control in cruising flight. Hence, the rms values of DC and DCD should 
be small, f = 10. seems to produce reasonably small rms values for DC 
and DCD, while the DC control activity is only modestly increased. 

Another example case is summarized below for the 100 kt cruise flight 
condition. This example differs from the previous one in that airspeed 
control replaces pitch attitude control as one of the ultimate outer-loop 
goals. Relationships among the control and output variables and cost 
function coefficients are summarized in Table 38. Table 39 is a summary 
of the resulting regulator design. 

TABLE 38 

RELATIONSHIPS AMONG CONTROL AND OUTPUT VARIABLES 
AND COST FUNCTION COEFFICIENTS 


Output 

Variable 

Bandwidth 

Requirement 

(rad/sec) 

Control Variable 

Association 

DBD 

DCD 



R R 0) = i-oo 

Rr(2) = 10.0 

5 B (DB) 

25-7 

Qr(4) = 662. 


5 (DC) 

25.7 


Qr(5) = 6620. 

9 (TH) 

1-5 

Qr(5) = 354,742. 


h e (HDE) 

0.70 

(See note) 

(See note) 

/ h e dt (HDI) 

0.32 


0 r(6) = 2.84 

u AS e ( ASE ) 

0.10 

Qr(9) = 3.44 


■f u AS e dt ^ ASI ) 

0.02 

Qr(8) = 1.37X10 -3 



Note: Bandwidth requirement is satisfied without additional feed- 

back. This is indicated by KDE/DC = 7-3 (n = —1), and 
HDE/TH = 169. (n = -1). 
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REGULATOR GAIN MATRIX 



An interesting feature in this 100 kt cruise example is that hand- 
width requirements on HDE are met as the result of imposing prior hand- 
width requirements on TH, DB and DC. This is indicated hy the fact that 
the asymptotic approximations to the transfer functions (obtained in 
Steps 4 and 5 ) for HDE/TH and for HDE /DC are pure gains. Thus, 


HDE 

TH 


- 169 . 


K = -169. 
n = -1 


HDE 

DC 


T .30 


K = 7. 30 

n = -1 


n = — 1 indicates a pure gain transfer function approximation. When this 
occurs, no cost function weighting (i.e., a value of zero) should he 
used on the variable under consideration. 

Generally, n < 0 indicates either that the bandwidth requirement is 
met (n = — 1), or that a prior bandwidth requirement has, in effect, used 
the current variable under consideration as an inner-loop feedback 
(n < -1). 

D. CONTROLLER SOLUTION 

Synthesis activities to this point have generated all the data 
required for the controller solution. The data merely need to be com- 
bined in the appropriate way. 

The structure of the controller is shown in Fig. 1^. Equations for 
the coefficient matrices (Ap, Bp, Cp and Dp) are given on page 67 . The 
’’controller solution” is simply an evaluation of these four coefficient 
matrices. These matrices, together with the controller equations in 
Fig. 1^, enable a circuit designer or software specialist to implement 
the controller . 

Controller coefficient matrices for the longitudinal UH-1H are given 
in Table 40 for the hover flight condition and in Table t-1 for the 
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TABLE kO. CONTROLLER COEFFICIENT MATRICES, UH-1H, HOVER 


y = A F y + B F z 
u = C F y + D f z 


AF MATRIX 

12 3 

E 1 3 E01 E 02 


! -0 - 33 3 

0. 906E-09 

13. 363E-13 8 ! 

E 1 3 



i 

2 

-34.3 

-3.922 

-2.37 ! 

E01 




3 

-1 . 60 

-3-76 

-15.0 1 

E02 



_* 



BF MATRIX 
1 

2 

3 

4 5 

6 

7 

8 

9 


HD 

Q 

TH 

DB DC 

HDE 

HDI 

XD 

XD I 

_ 

-0.146E-06 

4.97 

0. 310E-0S 

0. 1422-09 -0. 922E-J7 

0.4 52E-1 5 

0.000 

-3. 1O4E-0B 

0.000 

! 1 
: ei 3 
1 

95.6 

-29. 1 

788. 

-34.3 59.7 

-0. 295E-06 

0. 000 

-9.13 

0.000 

l 2 
! E01 

592. 

-76.5 

-0. 123E+04 

-1.60 381. 

-0. 183E-05 

0. 000 

1.36 

0.000 

1 3 
1 E02 


CF MATRIX 


1 

2 

3 


E 1 3 

E01 

E 02 

1 1 

1 

oo 

UD 

-3.254 

-0.523 

! DBD 




! 2 

0.122 

0.264 

1.07 

1 DCD 


DF MATRIX 
1 

2 

3 

4 

5 

6 

7 

8 

9 


HD 

Q 

TII 

DB 

DC 

HDE 

HDI 

XD 

XD I 


21.5 

538. 

0. 149E+04 

-29.9 

0. 518E-01 

-0.717 

-0. 252 

-26. 6 

f 

-1.89 i 

i 

i 

DBD 

-42.5 

1.30 

79. 8 

i3. 517E-02 

-27. 0 

4.16 

2.16 

-0.436 

-fl. 221E-01! 

2 

DCD 
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TABLE 41 . CONTROLLER COEFFICIENT MATRICES, UH-1H, 100 KT 


y = Ajy + B^z 
u = Cj>y + D f z 

AF MATRIX 


1 

E 1 3 

2 

Eei 

3 

E02 


-0. 333 

-0.614E-08 

i 

O.130E-061 

1 

E13 

B.06 

-a. 183 

I 

1.52 I 

2 

E01 

-125. 

1.17 

l 

-25.0 1 

3 

E02 


BF MATRIX 


1 

2 

-3 

t 

5 

6 

7 

8 

9 


HD 

0 

TH 

DB 

DC 

HDI 

HDE 

AS I 

ASE 


•0. 314E-05 

4.97 

0. 530E-03 

0. 6S7E-06 

-0. 171E-05 

0. 263E-12 

-0. 180E-12 

-0.908E-11 

0. 1 62E-06 1 
l 

i 

E13 

-36. 3 

261. 

0.613E+04 

8.06 

-20. 1 

-0. 14SE-04 

-0. 157E-03 

0. 478E-03 

! 

1.87 ! 

1 

2 

E01 

594. 

-0.419E+04 

-0. 100E+06 

-125. 

325. 

-0.422E-04 

0. 546E-04 

0. 147E-02 

! 

-30.6 ! 

3 

E02 


1 


CF MATRIX 


1 

2 

3 


E13 

E 8 1 

E02 

! 1 

-23.9 

-0.558 

-1.22 

1 DBD 




i 2 

0.934 

-0.313 

0.479 

! DCD 


DF MATRIX 


1 

2 

3 

4 

5 

6 

7 

8 

9 



HD 

0 

TH 

DB 

DC 

HDI 

HDE 

AS I 

ASE 

_ 


28.1 

271. 

-0. 434E+04 

-27.9 

0.627 

-1.34 

-2.99 

-0 . 193E-01 

-2.15 

l 

1 

1 

1 

DBD 

-11.6 

-2.14 

0. 206E+04 

0. 627E-01 

-25.9 

0.324 

0.981 

-0. 786E-02 

0. 281 

! 

1 

2 

DCD 
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100 kt cruise flight condition. It must be appreciated that the con- 
troller function is different in one respect for each flight condition 
(groundspeed control in hover in distinction to airspeed control in 
100 kt cruise). At each flight condition only one controller function 
(flight control system mode) is addressed. Other controller functions 
would be developed in similar form by similar methods. Different meas- 
urements may be used for different controller functions as in the case of 
the two example functions given (proportional -plus-integral groundspeed 
error in distinction to proportional -plus -integral airspeed error). 

Controller transfer functions are summarized in Table 42 for the 
hover flight condition and in Table 43 for the 100 kt cruise flight con- 
dition. The characteristic polynomial roots are the same as the closed- 
loop filter -observer roots for these special cases (refer to Tables 17 
and 22). This is for the reason that the [G-| — K-]2^2l matrix is null for 
these cases as was explained earlier. 

The numerators of these transfer functions are not particularly 
remarkable in most respects. However, they are not exactly what a con- 
ventional frequency response design approach would lead to either. This 
is particularly so with respect to the non-minimum phase zeros at low 
frequencies (e.g., the DCD/TH numerator in Table 42) and at mid frequen- 
cies (e.g., the DBD/TH numerator in Table 43). The non-minimum phase 
zeros at high frequencies do not really affect the design because of low 
loop gain at those frequencies. 

Another observation is that every measurement is fed back to every 
control. This results in eighteen controller transfer functions. While 
there are some controller transfer function paths for which the gain is 
negligible at all frequencies (e.g., the DCD/XDI numerator in Table 42), 
not many have this property. This is in contrast to the eight to ten 
non-zero controller transfer functions an experienced conventional fre- 
quency response designer would use. More will be said in the next chapter 
concerning controller simplification. 

The controller transfer functions are third order (not counting inte- 
grals of errors). This order is certainly acceptable as a practical 
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TABLE V3. CONTROLLER TRANSFER FTJNCTIONSj UH-1H, 100 KT 


CASE: UH1H 100KT 128LONG 2-FEB-79 CONTROLLER TRANSFER FUNCTIONS 
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I 


matter for flight control system implementation. It is estimated that 
an experienced conventional frequency response designer would use two 
first-order filters. 

This completes explanation of the controller synthesis procedures 
and examples. The next section explains and illustrates procedures for 
examining more of the system and controller properties and for conduct- 
ing the controller design assessment. These key activities are necessary 
to better understand what the optimal synthesis hath wrought, as well as 
to assure that other quantitative requirements and qualitative require- 
ments, which we are not able to incorporate directly into the problem 
formulation and solution procedure, are satisfied. 
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SECTION V 


DESIGN ASSESSMENT 


At this stage of the design process a controller has been synthe- 
sized that consists of a filter -oh server and a regulator. This section 
discusses the next step in the design process — examination and expan- 
sion of the synthesis results into a composite, total picture of system 
properties, behavior, and design margins . The goal here is to present 
the techniques needed to evaluate a preliminary FCS design. By exercis- 
ing these techniques, the analyst ■will be able to: 

• Determine whether design goals have been met. 

• Define quantitatively all the important properties of 
the system, e.g., 

— expand the scope of the results (e.g., to include 
transient characteristics). 

— highlight the system 1 s dominant properties. 

— reveal the sensitivity of dominant properties 
to uncertainties. 

• Broaden the understanding of the optimal synthesis 
design results, for example; 

— relative importance of the various terms In 
the control laws. 

— factors which contribute to particular con- 
trol behavior (e.g., see Section IV discussion 
of results when G-] — K-|2G2 = 0) . 

— implications for future design strategy. 

• Explore simplifications of the design which may lead 
to a simpler, more practical, controller. 

• Investigate the implications of controller architecture 
on results. 

• Generalize design for other applications. 

A marriage of evaluation techniques from classical and optimal control 
theory will provide the basis for the design analysis. 
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Table 44- presents a more specific outline of the system and controller 
properties that are examined in this chapter. All of the assessment goals 
listed above axe addressed by one or more of the table entries, although 
there is not always a direct correspondence between the two. The design 
goal list addresses rather fundamental questions, whereas the table (and 
the remainder of this section) presents a systematic checklist -like 
method for answering these questions while simultaneously defining and 
evaluating a given FCS design. 

To this end, notice that the topics addressed are divided into two 
major categories: properties of the system, and properties of the con- 

troller. Within each category a definite hierarchy exists, so that the 
evaluation proceeds from the simple and qualitative to the detailed and 
quantitative. Along the way, certain features may be deemed unimportant 
or uninteresting and are thus excluded from more thorough examination. 

In this way, the need for microscopic analysis of each aspect of the 
design is eliminated, and the number of assessment steps is correspond- 
ingly decreased. 

The remainder of this chapter is devoted to a close look at many of 
the techniques for design evaluation listed in Table J+4. We begin with 
an exploration of the stability characteristics of the system; stability 
is typically the most fundamental consideration in any feedback system. 
Certain metrics from classical control theory such as gain and phase 
margins are computed. These provide a quantitative rank-ordering of all 
variables in terms of their tendency towaxd instability. Then, for each 
variable, the mode which is the first to become unstable as the closed- 
loop gain increases is identified as the crucial mode. 

Next, the response characteristics of the system axe assessed. This 
task is divided into three areas: response of the primary variables, 

response of the secondary variables, and response of the controls. Cer- 
tain behavior is expected in a variable with satisfactory response: the 

transient response to commands must be smooth and well-behaved; the 
response to disturbances should be minimal; the control activity must 
be within reasonable physical limits; and command following should be 
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TABLE Ml-. DESIGN ASSESSMENT 
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Commonality of elements/ Common controller elements for each Controller gains over many opera- 

gain settings for different operational mode tional modes 

operational modes 







good. For both the primary and secondary variables the modes that domi- 
nate the response character are identified. 

The final consideration within the system properties category is the 
sensitivity of the design. Here, too, there are three areas of interest. 
The first is the sensitivity of the key response modes (i.e., the cru- 
cial and dominant modes) to off-nominal conditions, which is quantified 
with the use of first-order gain and parameter sensitivities. The second 
is the svstem sensitivity to parasitic nonlinearities such as different 
thresholds for the various feedbacks. Here, open-loop frequency response 
plots provide the data for assessment. Finally, the ability of the 
design to degrade gracefully is evaluated, based on its stability and 
response characteristics when loops are systematically opened. 

Consideration of properties peculiar to the controller follows. 

There are two major considerations: the gain levels and the sensor com- 

plex. Questions concerning control saturation and minimum increment of 
control are addressed within the context of gain levels. Under the 
heading of sensor complex fall many concerns which affect the controller 
mechanization. These include; equalization economy, which refers to 
the desire for minimum sensor /signal compensation/control channel fea- 
tures in the system implementation,- and commonality of elements, a desire 
for the minimum number of gain and/or element changes over the range of 
operating modes. Both economy of equalization and commonality of ele- 
ments are evaluated via variation of system architecture. 

Limited design assessments were undertaken during the course of the 
filter-observer and the regulator synthesis. Those aspects will be only 
briefly referred to at the appropriate points as the discussion winds 
its way through Table 44 . Certain other features which are enumerated 
in the table are beyond the scope of this report^ these too will be dis- 
cussed only briefly, without benefit of examples. In all other cases, 
however, the discussion of the evaluation techniques will be accompanied 
by examples from the UH-1H hover FCS design, presented in full in Appen- 
dix C. 
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A. STABILITY 


As stated previously, an understanding of the stability characteris- 
tics of a system is probably the most fundamental of all design assess- 
ment aspects. Interest lies not only in determining whether or not a 
given system is stable, but also in what causes the system to be stable 
and what modes of the system are most likely to become unstable when dif- 
ferences appear between the actual and analytical systems. Such infor- 
mation is needed to define management safety and to evaluate the per- 
formance of a design in off-nominal conditions. 

The most obvious place to begin the stability assessment is with the 
closed-loop system denominator. All closed-loop roots should lie in the 
left-half complex plane and indeed this is guaranteed for an LQG optimal 
controller if the open-loop plant is detectable and stabilizable and if 
a positive definite performance index was used in controller synthesis. 

In fact, these conditions assure uniform asymptotic stability (i.e., 
x(t) 0 as t oo ? see Ref. 24). One measure of the degree of stability 
is the damping of the closed-loop roots. Table 45 contains an example 
closed-loop transfer function (g/ 0 c )- Notice that not only are the 
denominator roots all stable, as assured because the uniform asymptotic 
stability criteria are obeyed here, but that the minimum damping ratio 
present is 0.674. 

Some additional observations can be made as to the system's dynamic 
behavior by further examination of the 9/0 c transfer function. Many of 
the poles cancel exactly (or nearly so) with zeros for this transfer 
function. These modes will, therefore, not appear in 0 response to e c 
commands. The resulting approximation to 0/0 c shown at the bottom of 
Table 45 reflects the closed-loop actuator characteristics (at 25-73) 
and the well -damped short-period mode represented by the quadratic. 

At this point, one might be tempted to conclude the stability assess- 
ment, by correctly observing that the closed-loop denominator has per- 
mitted identification of all the important aspects of stability. Yet 
there is still much which needs to be said about a system's stability 
beyond simple closed-loop root considerations. In particular, we are 
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TABLE k-5 

CLOSED-LOOP e/0 c TRANSFER FUNCTION 
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interested in two things: general margins and the modes that would he 

the first to go unstable as likely system parameters are changed. 

The first concern is the particular margins of stability for each 
loop. There are two common margins: the gain margin, which is the dif- 

ference between the gain at the crossover frequency and the gain at the 
—l 8 o deg phase point,* and the phase margin, which is the difference 
between the phase at crossover and — 180 deg. Figure 19 is an open-loop 
frequency response plot for 9 / 9 c , where x, / xdt, and 0 loops are open. 
The gain margin of 21 dB and phase margin of 56 deg reveal the excellent 
stability margins for this loop. 

In addition, it is instructive to convert the phase margin, cp^, to 
an equivalent t in a pure time delay, e jayr . In physical terms this 
delay, if added to the system, would be just sufficient to destabilize 
the system. The delay itself can be thought of as the low-frequency 
measure of a large number of high-frequency leads and lags which are not 
accounted for in our system model. If these are given by 


then 


t = £ T i “ 2 T i 


Ghi( s ) 


ni(TiS + 1 ) 
+ 1 ) 


That is, the delay t is the sum of the high-frequency lags minus the 
sum of the high-frequency leads. For the 9/9 c loop presented in Fig. 19 
the allowable t is computed as 

56 / 57-5 

o> c 2.3 

= 0.1125 sec 


This can be considered a delay or lag margin. 
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CASE: UH1H HOVER XD,XDI AND TH LOOPS OPENED 
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Figure 19- Frequency Response for Open-Loop e/0 c 
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Consider now the second concern^ namely, given the gain and phase 
margins for a loop, which mode will become unstable first as the gain mar- 
gin approaches zero or added lags create a time delay which approaches t? 
This question can be answered with the use of an s-plane root locus plot 
or a Bode root locus plot. Both are shown for the 0/0 c example in Fig. 20. 
The root locus is a plot of the open- loop poles and zeros (symbols X 
and 0 , respectively), and the locus of closed-loop roots. In Fig. 20a 
we see the short-period root locus heading for the right-half plane as 
the closed-loop gain increases. Thus, it is the short-period mode which 
will become unstable first when the gain margin is zero. This is con- 
firmed by the Bode root locus plot of Fig. 20b. Here the Bode plot and 
root locus plot are combined. The important open-loop poles and zeros 
are called out, as well as the closed-loop roots (the solid rectangular 
symbols). The crosses which depict the second-order closed-loop roots 
are shown moving into the right-half plane (at u> = 10 rad/sec) for the 
short-period mode, thus identifying it as the "crucial" mode for the 0/0 c 
loop. 

This completes the stability assessment for the 0/0 c loop. The same 
procedure should be followed for every other loop in the design. That 
is, examine the closed- loop roots, compute the stability margins, and 
identify the crucial modes. The results of this process for the h/h c > 
/h e dt/h c , x/x c , fx dt/x c loops are presented in Appendix C, Figs. C-8, 

C-9, C— 10, and C— 1 1 , respectively. Notice that the h/h c loop has a gain 
margin G]y[ = As a practical matter, this is an impossibility and simply 
indicates that one or more modes (and, indeed, one which is the crucial 
mode) has been neglected. Although the crucial mode for this loop cannot 
be identified, it is always possible to compute an effective t from phase 
margin and crossover frequency which represents a crucial value for the 
high-frequency leads and lags. 

Many other techniques are available in the literature for stability 
assessment, for example the Routh and Hurwitz tests (Ref. 2 6 ) and the 
Nyquist stability criterion (Ref. 27)- The methods used here have been 
selected to give the best mixture of meaningful results from a combined 
optimal/conventional standpoint. 
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B. RESPONSE 


Turn attention now to the system's response characteristics. For 
this aspect of assessment the system's output variables are divided into 
three classes: 

• Primary controlled variables perform the explicit 
commanded functions of the FCS. For the UH-1H 
hover example these functions are rate-of-climb 
command and groundspeed hold, so the primary varia- 
bles are h and x as outputs or he and x e as error 
quantities. 

• Secondary variables are those vehicle output varia- 
bles which remain after identifying primary varia- 
bles. The UH-1H hover secondary variables include 
0 , q, a, and u ASe - 

• Controls, 83 and 5 C , for example, comprise the final 
category. 

The next three sections discuss the assessment of response characteris- 
tics for each category of variables. 

1 . Primary Controlled Variables 

The most fundamental single indication of controlled variable 
response behavior is the response/command bandwidth. Among other things, 
the bandwidth; 

• Indicates the range of command frequencies over 
which feedback acts to provide good command follow- 
ing and error suppression. 

• Indicates the frequency range over which disturbances 
are suppressed. 

• Specifies the speed of response, since the bandwidth 
is inversely proportional to the indicial response 
rise time. 

Because bandwidth properties play such an important role, appropriate 
bandwidth values were made the key design requirements. Further, the 
output variables (to be included in the regulator cost function) were 
chosen to insure desirable low-pass properties and particular time 
response behavior (e.g., including / h e dt guarantees that the indicial 
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response in h is such that h e approaches zero in the steady state). Thus, 
"bandwidth" will have its traditional meaning in a low-pass filter con- 
text. 

To be certain that bandwidth goals were met, we examine opened-loop 
frequency response plots such as presented in Fig. 19- Here, for the 
0/0c loop, the goal was a bandwidth of 2.0 rad/sec. The stretch of 
—20 dB/decade slope on either side of the 0 dB crossover frequency (coc) 
permits us to approximate the bandwidth by ay>. Similar bandwidth assess- 
ment is indicated in Figs. C-8 through C— 1 1 for the primary variables. 

We have previously emphasized that the bandwidth numbers partition 
the response into regions of action — regions where feedback is useful 
in improving the response and regions where feedback effects are nil. 

The details of the system behavior in these regions are most conveniently 
addressed in the frequency domain, with the closed-loop response /command 
and error/command characteristics being of particular interest. For 
example. Figs. 21 and 22 present closed-loop Bode plots for the h/h c and 
h e /h c transfer functions, respectively. 

Several essential aspects of the h/h c response are apparent from 
Fig. 21 f There is good following of command input magnitudes at all fre- 
quencies below the bandwidth of approximately 1.4 rad/sec, as evidenced 
by the unity gain. Beyond that, the response to commands is smoothed 
and attenuated (integrated once) out to 25*7 rad/sec, and even more 
attenuated (integrated twice) thereafter (not shown). In addition, the 
phase shift up to the bandwidth Is nearly linear with frequency. This 
closed-loop system response is dominated by only one lead and one quad- 
ratic lag out to 2 6 rad/sec. In fact, the response can be approximated 
by 1.4/(.7 s+ 1) for all frequencies below 26 rad/sec. Thus, the response 
character of what at first glance seems to be a rather complicated trans- 
fer function can be reduced to the simple form K/(T-|S+ l). 

Similarly, Fig. 22 provides much information regarding the nature 
of the h e response to h c . Here again the frequency response has three 
distinct regions: for co < 0.1146 rad/sec it behaves as a zero position 

error system; for 0.1146 rad/sec < o> < 0.91 rad/sec it behaves as a zero 
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Figure 21 . h/hc Closed-Loop Frequency Response 




velocity error system; and beyond the bandwidth of approximately 1 .4 rad/ 
sec there is no error suppression (|h e /hc| = 1 ) . 

Both the h/h c and h e /hc Bode plots lead us to suspect that the mode 
that dominates the h response is, not surprisingly, the path mode 
(£ = 0.923, 00 = 0.90T)* It is this mode that provides the general shape 
for both plots. Indeed, Fig. 23 confirms this suspicion with a Bode root 
locus plot of h/hc, where the lie and / lie dt loops have been opened. The 
dominant mode is closest to the point of crossover, which approximates 
the bandwidth due to the stretch of —20 dB/decade slope. 

The discovery of the dominant mode for the h response leads naturally 
to an examination of the transient time response behavior of this varia- 
ble. Plots of both the indie ial and ramp responses of h to h c input are 
shown in Fig. 24. The path mode does, of course, dominate the response, 
as it did in the frequency response plots of Figs. 21 and 22. However, 
the initial part of the indicial response is quite first-order -like, even 
though the path mode is second-order. This is not surprising, since the 
almost -critical damping of this second-order mode (£ = 0.923) allowed us 
to approximate the frequency response by 1.4/(.7s+l) out to 2 6 rad/sec. 
The transfer function for h/h c contains a lead term (the path mode) as 
well as the quadratic lag, and the lead part of the combination provides 
the overshoot in the indicial response shape. Thus, initially the 
response is first order, with a time constant near then, after the 
first 0.4 sec the response is that of the heavily damped path mode, with 
the 12 percent overshoot attributed to the lead at 0-53 rad/sec. These 
time domain characteristics reflect the dominant features of the frequency 
response (Fig. 21). 

A similar comparison can be drawn between the h e /h c frequency response 
plot of Fig. 22 and the h time response to an h c ramp input shown in 
Fig. 24. Remember the three frequency regions which were present in 
Fig. 22 — a region of zero position error, a region of zero velocity 
error, and a region where the error was passed with unity gain. Those 
same divisions can be made. for the ramp response in Fig. 24. The h 
response at low frequencies looks like that of an equivalent first-order 
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system, 1.4 (T-,s-M), where T-| =0.7 sec. At mid frequencies the response 
exhibits the characteristics of a zero velocity error system, and begins 
to approach the £i c input ramp. But the response never quite reaches the 
command because the zero velocity error behavior is only active in mid 
frequencies. Thus, there is ultimately a non-zero offset of the response 
from the command. This is illustrated by the ramp delayed by T 2 seconds. 
Again, the characteristics of the h/h c ramp response are consistent with 
those we would predict from the h e /hc frequency in the frequency domain. 

The command -following properties of the primary variables can also 
be examined using the rms levels of the states, outputs and controls in 
response to process and measurement noise. These are presented in Table 46. 
Notice, for example, that the rms levels of h c and h are approximately 
equal (h c = 1.0 ft/sec, h = 1.01 ft/sec), which corroborates the command 
following behavior gleaned from both the frequency domain (Fig. 21) and the 
time domain (Fig. 24) . The ability of this system to suppress error can 
also be illustrated by comparing the rms level of £l c and h e (£i c = 1.0 ft/ 
sec), h e = 0.259 ft/sec). The rms level of error is approximately one- 
fourth of the rms command input, which is expected in light of our pre- 
vious observations regarding h e /h c in the frequency domain. 

Considerations of the error suppression character of the rate-of- 
climb loop lead naturally to an investigation of the system response to 
gust disturbances. As a general rule of thumb the design should be indif- 
ferent to gust disturbances, in order for the control system to perform 
satisfactorily in actual flight. As stated at the outset, the behavior 
of the system in the presence of disturbances is fundamentally a function 
of the bandwidth. As such, the gust suppression question has already been 
addressed in principle. An example of the h. behavior in response to lon- 
gitudinal and normal gusts is presented in Fig. 25- Gust suppression for 
this particular variable is excellent. 

This concludes a detailed example evaluation of the primary variable 
response. The example used for the rate-of-climb variable showed it to 
have satisfactory response characteristics in every response. To com- 
plete the story, a parallel evaluation of the groundspeed (x) variable 
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CASE! UH1H HOVER 122L0NG 31-JAN-79 ALL LOOPS CLOSED 
FLOATING UNDERFLOW FC=036035 


DENOMINATOR! 

1.0000 

( .10195 ) ( .31507 1 ( .33300 ) ( .33595 

< .50520 ) ( 15.614 > l 25.716 > <■ 25.726 

(( .92313 7 .90723 7 . 83743 7 .34832 )) 

(( .67330 7 2.7741 7 1.8692 7 2.0498 )> 

119.79 > 

FLOATING UNDERFLOW FC-036035 


NUMERATOR ! HD/ UG FILE NAME? HDUG.2CL 

NEW FILE 

-.98613E-01 
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( .48797 ) < 25.726 ) < 42.347 ) 
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(( .69689 , 2.7215 7 1.8966 7 1.9519 )) 
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h_ 
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w g (s + 15. 6) (s + 25.72) [s 2 + 2 (. 923) (.91 )S + ( .91 ) 2 ] 


Figure 25. Time Responses of h to Gust Disturbances 
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is necessary. Such an investigation was performed, and some of the results 
appear in Appendix C. Again, satisfactory response was observed, and the 
groundspeed mode at O.505 rad/sec provides the dominant behavior. 

In the process of these examinations much has been learned about the 
details of the optimal design. Also, the evaluation has shown that no 
iterations are needed as far as primary responses are concerned. 

2. Secondary Variables 

Techniques that were appropriate for the assessment of primary variable 
responses are again called upon at this stage. The bandwidth is still the 
single most revealing measure of the secondary variable responses. Refer- 
ring back to Fig. 19 * which depicts the open-loop frequency response for 
0/q c , we see that the bandwidth (equal to about o^, due to the stretch of 
—20 dB /decade slope) is 2.3 rad/sec, quite close to the design goal of 
2.0 rad/sec. While secondary variable bandwidth can be used to divide the 
frequency response into regions of feedback effectiveness and ineffective- 
ness (similar to the primary variable assessment), the expected forms for 
response/command and error /command frequency responses are changed. This 
is because the secondary variables are usually inner loops and do not have 
explicit command input points. For example, in the hover FCS we cannot 
command pitch attitude (0); the closed-loop design requires us to co m mand 
groundspeed (x) instead. If, however, the groundspeed loop should fail 
(due to component degradation, say) it would then be imperative that 0/0 c 
exhibit good command following, error and gust suppression and transient 
response behavior, since the system would have to be commanded by 0 C . This 
0 C system is also a possibility as a primary system in its own right. We 
will have occasion to discuss both these events in more detail later. 

For now, we will rely on the Bode root locus plot of Fig. 20 to show 
that the short-period mode, which is closest to (%, dominates the 0 re- 
sponse. This is hardly unexpected, since the 0 and q closures provided the 
improved damping and stiffening for the short-period mode (£ S p!q L = 0.26279* 
“sploL = °- 92797; £ sp l CL = °- 6 ^ 8 ' “spl CL = 2. 77^1). Reference 10, pages 
80-87 also indicates a graphical method for determining the dominant mode 
by plotting the modal response coefficients as time vectors. 
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With the dominant mode exposed, we proceed to an investigation of the 
transient response characteristics. Again, the absence of an explicit 
closed-loop e c command point makes indicial and ramp responses to commands 
meaningless. The real concern in the evaluation of secondary variable 
time responses is for motion harmony and favorable (or negligible) coupling 
between primary and secondary variables. For this reason. Fig. 2 6 pre- 
sents the responses of h, 5c> and §B to a 1.0 ft/sec h c step input. 

It is clear that 6 and h are almost totally uncoupled, since an h c input 
elicits very little response from 9 and a quite noticeable response from h. 
This phenomenon can actually be explained on the basis of control activity. 
Since mainly controls 0 while 5q primarily controls h, and the controls 
are almost uncoupled, we expect the 0 and h responses to be essentially 
uncoupled as well. This point is addressed in the next subsection. 

The same uncoupled behavior of h and 0 can be seen in their respective 
rms responses to process and measurement noise. Table b6 gives the rms 0 
responses as less than 0.1 deg for a 1.0 ft/sec rms h c level. This further 
confirms the high degree of h-0 decoupling. 

There is also interest in the response of secondary variables to gust 
inputs. Figure 27 reveals that the 0 inner loop is quite indifferent to 
gust disturbances. This is desirable, although not really imperative, 
again because 0 is not a controlled variable. But if the x loop should 
fail or be switched out, it is reassuring to know that this inner loop 
exhibits good gust suppression. 

The secondary variable 0 has demonstrated satisfactory response charac- 
ter during the several stages of its assessment. One important difference 
in the methods used to evaluate primary and secondary variable responses 
arose in connection with command following and error suppression. Unlike 
the primary variables, the secondary variables have no meaningful command 
points. We must therefore depend on open-loop transfer function charac- 
teristics to divulge the dominant mode(s). In addition, the concerns in 
the time domain center around motion harmony and the extent of inter- 
variable coupling (or decoupling) rather than rapid and accurate response 
to commands. 
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Figure 26. Transient Responses of Several Variables to h c Step Input 
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Figure 27. Closed-Loop Transient Response of 9 to Gust Inputs 
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3 . Control Activity 


The third area of response assessment has an emphasis different from 
the first two areas. Since the control actuation bandwidths are at high 
frequency ( 25.7 rad/sec for the UH-1H example), the command following, 
error suppression and time history characteristics are fairly uninterest- - 
ing. Instead, the evaluation of control activity centers around the fol- 
lowing two considerations ; the level of control activity required by the 
controller and the extent of control cross-coupling. 

The question of appropriate control activity levels was first addressed 
during the regulator synthesis. As the final step, a measure of control 
activity (positive zero -crossing rate) was used to ascertain whether or not 
the cost function weights on and 5 C in % were appropriate. As it turned 
out, the main rotor collective control point was too active, and thus the 
cost function weighting on 5c was increased by an order of magnitude. This 
10:1 ratio of SgrBg discourages use of the 5g control point in the regula- 
tor synthesis, and positive-going zero crossing rate for 5 q showed the 
desired decrease. 

There are several other measures of control activity level which can 
be used at this point in the response assessment. These are illustrated 
in Table 47. The first two measures are the rms levels of the control posi- 
tion and control rate. These values can be compared to the limiting control 
capabilities of the particular aircraft. For example, the UH-1H longitu- 
dinal cyclic has a 13.0 in. maximum travel, while the main rotor collective 
maximum travel is 10.7 in. We can see that the rms levels of S-g and 5 q are 
well below the corresponding maximum capabilities. Thus, if the assumed 
inputs are realistic estimates the actual control activity should be moder- 
ate. 

Table 47 also shows the number of axis crossings per second (N 0 ), the 
number of positive -going axis crossings per second [N;g(o)], and the number 
of positive exceedances of 10 percent full travel per second [2N L ( . 1 th) ] . 

The first two indicate the level of monitoring activity required of the 
pilot in supervisory or, for that matter, the manual activity needed if 
the pilot were to take over from the automatic system. These are 
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TABLE 47- CONTROL ACTIVITY 


Axis crossings per second 


No 



If CTx 


Positive exceedance of x 
per second 


N L (x) 


No e -x2/2 c x 
2 6 


oji = rms of control rate 
ti x = rms of control position 



Longitudinal 
Cyclic, DB 

Main Rotor 
Collective, 

Ox (in.) 

0.0494 

0.0921 

(in. /sec) 

0.0466 

0.255 

N 0 ( 1/sec) 

0.300 

0.881 

N L (0) (1/sec) 

0.150 

0.441 

N l ( . 1th) a ( 1 /sec) 

5.6 x 1 0 -9 

8.8 x 10 -4 


DCT 


a Positive exceedance of 10 percent full travel per 
second (for DB, 10 percent full travel = 1.30 in., 
for DC, 10 percent full travel = 1.07 in.) 


considerations which come under the heading of "handling qualities." The 
Np( . 1 th) measure can be used in a probabilistic sense to determine the 
average time between exceedances of 10 percent full travel [ 1 /2Np( . 1 th) ] . 

Up to this point, the Investigation of control response has not 
addressed the extent of control cross-coupling. Decoupled controls are 
desirable in that they permit control of variables from a single control 
point. This simplifies the controller both from an implementation and a 
performance standpoint. There are many quantities that reveal the nature 
of control cross-coupling^ two will be presented here. The controller 
transfer functions that involve control position and rate are given in 
Table 48 . The transfer function approximations are divided into low, mid 
and high frequency regions. It will be recalled that DB and DC are the 
measurement vector components corresponding to best estimates of 5g and 
Sq. As such, they are not the integrals of the control variables (DBD) 
and 5 q (DCD). Accordingly, to avoid confusion, the controller transfer 
functions in Table 48 are referred to as DBD/DB, etc. 

The low, mid and high frequency approximations for DBD/DB, and the mid 
and high frequency approximations for DCD /DC are essentially pure gains. 



This is the expected result of designing the actuators as high gain inner 
loops. The low frequency "washout present in DCD/dC r i.e., 27s/(s + 15.6) 
in Table 48] is undesirable in a physical installation (e.g., for actuator 
rigidity considerations) and may be considered an anomaly of the optimal 
design process. 

The cross-coupled transfer functions, and 5^/5g_, illustrate a 

second point, i.e., that the controls are essentially decoupled for the 
frequency region above 100 rad/sec. In this frequency band, both transfer 
functions revert to low gain quantities which may be neglected. In an 
upcoming section we will discuss the advantage of this control decoupling, 
namely, it provides a means of simplifying controller implementation. 

A further confirmation of the cross-coupled control characteristics is 
found by referring bach to the indicial transient response plots in Fig. 2 6 . 
The time responses of five variables to a 1.0 ft/sec h c step input were 
first used to examine the harmony of the primary and secondary variables. 
These time responses also illustrate the degree of control decoupling. For 
a step h c we see a corresponding response, in the positive direction. 
Meanwhile, no 5g response is evident. From this we can conclude that Sq 
provides the h control, and that 6g and 5 q are essentially uncoupled. 

Similar time responses to, say, an x c would show corresponding responses 
for x, '0 and 5g, illustrating that 5g is used in x and 0 control and again 
that 6g and 5 q are not significantly cross -coupled. 

This concludes our discussion of the design's response character. We 
have divided the preceding discussion into three areas; primary controlled 
variable responses, secondary variable responses and control responses. 

The primary variables were characterized by good command following, error 
suppression, indifference to gusts, and smooth, well-damped, minimum-tailed 
transient responses. For each primary variable we identified the mode that 
dominates that variable’s response. The assessment of the secondary varia- 
ble's responses was mainly concerned with motion harmony and coupling among 
the primary and secondary variables, and identifying that mode for each 
variable which dominates the response. Finally, the evaluation of the con- 
trol activity centered on the levels of control activity and the extent of 
control cross-coupling. 
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TABLE 48. CONTROLLER TRANSFER FUNCTIONS 
AND APPROXIMATIONS 


CASE: UH1H HOVER 122LONG 2-FEB-79 CONTROLLER TRANSFER FUNCTIONS 


DENOMINATOR: 


1 . 0000 

( .31537 ) { .333*3 ) ( 15.614 ) 

< 1.6383 > 


NUMERATOR: DBD/ DB 

-29.879 

{ .18186 ) ( .333*3 ) ( 15.423 ) 

<-27.917 > 


NUMERATOR; DCD/ DB 
. 51744E-02 

(-.11656 ) ( .33303 ) (-2064.8 ) 

< .41471 > 


NUMERATOR: DSD/ DC 

. 5178 IE -31 

( .29977 ) ( .33300 ) (-4127.1 ) 

<-21.332 > 


NUMERATOR: DCD/ DC 

-27.008 

. 94 888E-0 1 ) ( .33300 ) ( . 33773 ) 

< .28821 > 


LOW FREQUENCY 
APPROXIMATION 
(0.54 <. go < 15*6 
rad/sec) 


MID FREQUENCY 
APPROXIMATION 

( 15.6 < ai < 2064 . 

rad/sec) 


HIGH FREQUENCY 
APPROXIMATION 
(2064. <, cd) 


DBD -29-879(s + 15.45) 

DB (s +15-6) 

= -29.553 -29.879 -29.874 


DBD .052(s -4127- 1) .052(s -4127. 1) 

DC (s + 15.6) s 


= -13.757 


— 2 1 4 . 6/s 


.052 


DCD 

DB 


. 0052 ( s - 2064 . 8 ) 

(s + 15 . 8 ) 


. 0052 (s - 2064.8) 

s 


-.685 


-10.74/s 


.0052 


DCD -27.008.5 

DC (s + 15.77) 

= - 1 . 73 s 27.008 27.008 
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C. SENSITIVITY 


The design assessment, in addressing questions of stability and re- 
sponse, has exposed performance characteristics only for the system as 
modeled and only within the intended operating environment. It is equally 
important, however, to analyze the effectiveness of the design in various 
off-nominal situations. Thus, we wish to examine the system’s sensitivity 
to airframe variations, sensor thresholds, major component failure and 
the like. Following Table kh, these examinations are divided into three 
areas ; 

• Sensitivity of key response modes. 

• Sensitivity to parasitic nonlinearities. 

• Graceful degradation. 

The next sections discuss each area in turn. 

1 . Key Response Modes 

The assessment of system stability and response led to identification 
of some key response modes for each variable; the crucial mode, which 
becomes unstable first^ and the dominant mode, which shapes most of a 
variable's response. These may or may not be the same mode. Here we are 
concerned with the sensitivity of these key modes to off-nominal condi- 
tions, component tolerances and variations, and state equation uncertain- 
ties, nonlinearities and simplifications. It is possible, of course, to 
perform a sensitivity analysis on all modes of each variable. But the very 
fact that we are able to identify key modes means that other modes are 
less important and also less sensitive. Thus, their assessment may be 
deferred until a detailed design evaluation is undertaken. 

Because "tolerances," "variations," "off-nominal conditions," etc., 
all tend to be small, we can enlist first-order gain and parameter sensi- 
tivities, which derive from small perturbation considerations, to quantify 
the key response mode sensitivities. First-order gain sensitivities com- 
pute incremental changes in closed-loop pole locations due to incremental 
changes in open-loop gains. This gain sensitivity is equal to the 
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negative of the modal response coefficient for the particular closed-loop 
root, and is a measure of the slope along a root locus. As such, it is 
meaningful to plot the sensitivities on a root locus, as shown in Fig. 28. 
Here the first-order sensitivity of the short -period mode to incremental 
Kq changes is plotted on the e/e c root locus. Remember that the short- 
period mode is the critical mode for the 0 loop. Also listed at the bot- 
tom of Fig. 28 are the modal response coefficients for every other root 
of 0/e c . In general, the incremental gain sensitivities for these other 
roots are quite small, reinforcing our assertion that only the key response 
modes have significant sensitivities. 

First-order gain sensitivities are useful in evaluating the effect of 
changing the closed-loop gain upon the location of the closed- loop roots. 

If, however, we wish to assess the effect of changing parameters of the 
controlled element upon the root locations, we can use parameter sensiti- 
vities. The software (Ref. 23) can be used to compute the sensitivity of 
each regulator eigenvalue to incremental changes in the F, G and C matrices. 
This reflects the sensitivity of the closed-loop roots to any uncertainties 
in the various stability derivatives, for example. Figure 27 illustrates 
the parameter sensitivity of the path mode for the h variable. The lower 
portion of the table is the entire matrix of sensitivities of the path 
mode (the dominant mode for h) to changes in the components of the F matrix. 
These sensitivities are normalized by the path mode eigenvalue, and are 
expressed as a magnitude and phase (in degrees). We can plot element f( 6,5) 
(outlined in the sensitivity matrix) on a root locus for h e /h c , to indicate 
the direction in which the path mode would move for an incremental change 
in the stability derivative. (Table B-2 in Appendix B shows that the 
element F(6, 5) is the 1% derivative.) Each non-zero element in this sen- 
sitivity matrix could be plotted as a vector on the root plot in the com- 
plex plane to produce a graphical representation of the parameter sensitivi- 
ties of the system design. This provides a means to account for undertain 
specification of the controlled element, and the effect of this uncertainty 
upon the closed-loop roots. 
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Figure 28. First-Order Gain Sensitivities for 8 Loop 
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NORMALIZED p MATRIX SENSITIVITY TO EIGENVALUE : -0,017 0.349 

1 2 3 4 5 6 7 8 9 1 0 11 1 2 

X01 X02 X03 X06 X 07 X 08 X09 X10 Xll X12 X13 X14 


! 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 J 1 

I 90.0 90.0 90. 0 90.0 90.0 90.0 90.0 90. 0 90.0 90.0 90.0 90.0 | X01 

I I 

! 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 , ? 

I 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 i X02 


I 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1 3 

: 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 I X03 


i 0.000 0.000 0.000 0. 150E-02 0. 49SE-02 0.298E-01 0.749 1,03 0.629E-02 0.000 0.131 0.000 1 4 

! 90.0 . 90.0 90.0 -16. 3 1 G7. -51.2 -3.57 -94. B 161. 90.0 -1 78. 90.0 I Xi'f. 


t 0.000 0.000 0.000 0.664E-02 0. ll'S-01 0. 191E-01 0.790E-02 0.476E-01 0.135E-01 0.000 0.603E-02 0.000 , 5 

! 90.0 90.0 90.0 -129. -126. -164. -116. 152. -131. 90.0 69.2 93.0 ! X07 


0. 662E-03 0.000 0.462 0.000 \ 6 

-1 29. 90.0 -108. 90.0 J X0B 


0.000 0.000 0.000 0.000 | 7 

90.0 90. 0 90.0 90.0 r X 09 


! I 

1 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 J 8 

I 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 I X10 


! 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 I 9 

I 90.0 90.0 90.0 9J.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 90.0 t Xll 


I 0.000 0.000 0.000 0.4O4E-02 0.308E-01, 0.000 0 . 6 1 4E-02 0.000 0.000 0. 000 0.000 0.000 ’■ 10 

I 90.0 90.0 90.0 -120. 62.7 90.0 -108. 90.0 90.0 90.0 90.0 90.0 I X12 



I 0.000 0.000 0.000 0.000 0.000 0.425E-01 0.000 0.000 0.000 0.000 0.5S3E-02 0.000 "l 11 

! 90.0 90.0 90.0 90.0 90.0 19.4 90.0 90.0 90.0 90.0 72.6 90.0 I X13 


I 0.000 0.000 0.000 0 . 355E-01 0.11GE-02 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1 12 

! 90.0 90.0 90.0 161. 164. 90.0 90.0 90.0 90.0 90.0 90.0 90.0 i X14 


Figure 29- Parameter Sensitivity. of Path Mode to 
Changes in for h 


160 



We have discussed some possible techniques for evaluating first-order 
sensitivities in order to assess the behavior of the key response modes 
due to various off-nominal conditions. In particular, gain sensitivities 
are useful in computing incremental changes in closed-loop root location 
due to changes in closed-loop gains, while parameter sensitivities reveal 
the effect of controlled element parameter variations upon closed- loop 
root location. Only two key modes were presented as examples,- in the 
course of the assessment it is necessary to examine the sensitivity char- 
acteristics of all key modes. 

2 . Parasitic Nonlinearities 

In constructing a linear set of state and output equations for optimal 
controller synthesis, various nonlinearities were neglected. The assess- 
ment of key response mode sensitivity can be useful to quantify the first- 
order effects of some of these neglected nonlinearities. In a detailed 
design evaluation, however, it is necessary to have a more complete under- 
standing of these parasitic nonlinearities and their effect upon system 
performance. Because this is a detailed design task, we will give only an 
outline of the assessment procedure here. 

One basic problem arises when certain types of nonlinearities produce 
limit cycle behavior. The most common types of parasitic nonlinearities 
are threshold -like phenomena in the sensor and actuation/load elements. 
These can cause problems in two ways; 

• If one of the control points has a higher threshold 
than the others, then for some small region in the 
control space the feedback to that control point is 
not active. For the UH-1H, suppose that the thresh- 
old for 6g is larger than that for So Then for some 
input levels, only the Sq control is active. 

• If each of the sensors has a different threshold, 
then in some region of the control space certain 
sensors will not be sending useful signals. Thus, 
although all control points may be active, several 
of the sensors may have outputs equal to zero and 
the system is simply incomplete, possibly even unsta- 
ble in the threshold region. 
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As an example, suppose that the h e and / he dt sensors have higher 
thresholds than all other UH-1H sensors. Examining the root locus and 
transfer function for the 0 closure (Fig. 30), we see that, for some gain 
values, the 6 closure has right-half plane roots when the h e and / h e dt 
loops are open due to thresholds. This would lead to low-amplitude limit 
cycle behavior. 

It is therefore imperative to examine possible combinations of open 
and closed loops in order to identify those which exhibit oscillatory 
behavior. The use of Bode plots In conjunction with describing function 
techniques (Ref. 26) serves to estimate the resulting limit cycle ampli- 
tudes. Of course, the details of the limit cycle itself are dependent on 
the nature of the nonlinearity. We have given an example of the problems 
that occur due to differential control point or sensor thresholds. There 
are two additional categories of parasitic nonlinearities that should also 
be considered in practical flight control systems: hysteresis, due to 

cable/pushrod-Coulumb friction interations^ and valve friction operation 
within the actuator's closed loops. These are discussed in greater detail 
in Ref. 28. 

3. Graceful Degradation 

Above we have examined the problem of design sensitivity on a small 
scale (that is, the sensitivity to relatively small effects such as off- 
nominal conditions and parasitic nonlinearities). Here we consider major 
component failure, to ascertain whether the system performance will degrade 
gracefully. This is a major concern in the operating environment, and 
should be a consideration even within the preliminary design phase. 

As one might expect, the techniques used for assessing graceful degra- 
dation depend on the controller architecture. Here we consider two pos- 
sible controller mechanizations, which comprise two limiting cases. The 
first is a single thread system, wherein the controller gains are fixed 
and satisfactory performance is obtained only when all elements are work- 
ing correctly. These are typical of the majority of flight control sys- 
tems today. However, when fail-safe or fail-operational requirements are 
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present, modern FCS mechanizations use two-, three-, or four-fold 
("threads”) redundancy. 

Considerations for a single-thread design's graceful degradation are 
quite straightforward, and in fact very similar to techniques used in 
determining the design's sensitivity to parasitic nonlinearities. We 
noted that for some small regions in the control space certain control 
points or sensors could he inoperative due to differential trehsolds. 

In an extension of this idea, we can identify some large regions of the 
control space wherein control points or sensors are imperative due to 
component failure. Thus, the evaluation of the effects of differential 
thresholds and component failure in single -thread installations can be 
conducted in parallel. For the former, the differential thresholds and 
the control region affected are usually small, whereas for the latter the 
differential thresholds are infinite and the control region affected is 
quite large. The analysis is identical, however. 

The real drawback in assessing a single-thread architecture is the 
enormity of the task. It is necessary to make an exhaustive catalog of 
likely failure possibilities and examine each in detail. For a multiply - 
redundant system, the necessary work is far greater still, if fail- 
operative performance is to be guaranteed. 

At the other extreme from the single -thread system is a system employ- 
ing an adaptive filter -observer controller architecture. Such a control- 
ler is automatically reconfigured by an onboard computer if one or more 
components should fail. It is the computer’s job to detect the failure 
and recompute the necessary controller gains, similar to the controller 
synthesis discussed in Chapter IV. For this type of mechanization the 
workload in analyzing graceful degradation is substantially reduced. The 
system is easily assessed using the software developed in this work,* the 
task consists of systematically deleting sensors and control points, and 
using the optimal controller software to synthesize the resulting con- 
troller gains. Each of these "degraded" designs must be subjected to 
stability and response assessment, in order to quantify the effects of 
degradation. This approach does not take into account the software and 
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hardware which must detect failures and recompute controller gains, nor 
does it account for transients in the process. The method will, however, 
provide an upper hound on degraded system performance. The matter of 
system behavior in the transient region during the controller recomputa- 
tion is perhaps best evaluated using the actual controller. 

Concluding Remarks on System 

Property Assessment 

The assessment of system properties illustrated here has been organ- 
ized so that a logical progression exists from simple to complex, pre- 
liminary design to detailed design, qualitative to quantitative. In the 
process a much clearer understanding of the essential nature and behavior 
of the optimal controller has been gained. In fact, there has been a 
quantum increase in our understanding of the optimal controller synthe- 
sized In Chapter IV,* there we had only four matrices that provided infor- 
mation about the design (Ap, Bp, Cp and Dp) . Now we are equipped with 
much more insight into the stability, response and sensitivity character- 
istics of the system. The remainder of this chapter is devoted to assess- 
ing certain properties that are peculiar to the controller. 

In the following sections we hope to better define and improve our 
understanding of the optimal controller itself. Attention will be focused 
first on the gain levels produced by the design. Our concern here lies 
both with too high and too low gain values — with control saturation as 
well as minimum increment of control. These extremes, taken together, 
provide measures of how well the controller uses the available dynamic 
range of a channel. The second property to be considered relates to the 
sensor array. In this connection we desire to exercise hardware and soft- 
ware economy in the sensor/equalization complex and to evolve a controller 
that will perform satisfactorily over a wide range of operating environ- 
ments. The discussion of these issues is not meant to be exhaustive, but 
it does cover those controller properties that should be evaluated in the 
preliminary design phase. 
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D. GAIN LEVELS 


The controller gain levels provide some important insight into the 
behavior of the FCS. Obviously, we are interested in both limiting 
cases : 

• What is the magnitude of each output that corre- 
sponds to saturation of each control point? 

• What is the magnitude of each output that corre- 
sponds to the minimum increment of control for 
each control point? 

Saturation levels permit us to identify boundaries in the control space 
which separate the space into regions where control can be exerted to 
modify the basic vehicle dynamics. Outside of these regions, the dynamic 
performance characteristics are those of the vehicle alone. On the other 
hand, the minimum increment of control establishes the control precision. 
The levels of the output variables which correspond to the minimum incre- 
ment define an inner region of the control space wherein control is 
effectively not present. Inside this region the system dynamics are, 
again, those of the vehicle alone. 

The assessment of gain levels involves a tradeoff between too high 
gains, which cause the control boundaries to be reached with small inputs 
(bang-bang control) and too low gains which require large inputs to even 
force entry into the available control space. The level required to satu- 
rate a control also depends on the type of actuators^ full authority paral- 
lel servos generally saturate at a larger percentage of maximum travel 
than restricted authority series servos. Regardless, the first step in 
determining the effective control space boundaries consists of defining 
the control point saturation levels. For the UH-1H, with its full author- 
ity parallel servos, the maximum longitudinal cyclic travel is 13-0 inches, 
while the maximum main rotor collective travel is 10.7 inches. 

Next we turn to the matrix, shown in Table 49, which contains the 
gains between each output and both control points. The units are also 
indicated. Since the elements of the control vector are control rates, 
and the saturation levels sire control positions, it is necessary to divide 
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each element of the Dp matrix by the appropriate gain between the control 
position and rate. For example, each gain between the outputs and 6g is 
divided by 29.9, the gain between Sg and 5g. The gains are thus expressed 
as inches of control per units of output. These normalized gains are 
listed in Table ^9* 

Dividing the saturation levels by the normalized gains gives the mag- 
nitude of each output which corresponds to saturation of the Sg or &c 
control point. The smaller of the two numbers for each output defines 
the control space boundaries within which linear control can be exerted 
to modify vehicle dynamics. These quantities may also be used to elimi- 
nate the ineffective gains from the Df matrix. 

We can also compute the output magnitude that corresponds to the mini- 
mum increment of control, by dividing the minimum increment of control by 
the normalized gains. The term "minimum increment of control" refers to 
the smallest amount a control effector or surface will move, regardless of 
input level. The minimum increment of control exists because of friction 
and other thresholds, and defines the precision with which a given control 
surface can be positioned. A typical power amplifier -actuator-load combi- 
nation has a dynamic range of 100-300^ thus, we have chosen to use l/l6 in. 
as the minimum increment of control for the computations in Table h9> 

The resulting output variable magnitudes that correspond to the mini- 
mum increments of control define the inner boundary of the linear control 
space. They can also be used to eliminate ineffective controller gains. 

For example, we see that using Sg to control the pitching velocity, q, 
provides an accuracy of approximately 0.2 deg/sec. If Sc were used, the 

accuracy is 78 deg/sec I The main rotor collective is obviously not useful 

Sq 

in commanding pitching velocity and thus the gain could be eliminated. 

Using the above reasoning, we can eliminate about half of the con- 
troller gains, which results in the simplified controller of Fig. 31 . 

One should be certain to determine that this simplified structure does 
not differ significantly from the one synthesized in Section IV by com- 
paring the stability, response and sensitivity characteristics of the two 
designs. The desire for a simple system leads naturally into the sensor 
array assessment of the next section. 
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Figure 31 . Simplified Controller Gains 
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E. SENSOR /EQUALIZATION COMPLEX 


This is the third look at assessment of the sensor complex. It was 
first considered "back in the filter -oh server synthesis, when the rms state 
estimation error was compared to the desired error budget. If the budget 
was exceeded, the sensor array was clearly not able to provide an accurate 
est ima tion of the states. Alternately, if the estimation error was well 
within the budget, sensors might be deleted in the name of cost effective- 
ness and other such considerations. 

The second pass at assessing the sensor array occurred during evalua- 
tion of the controller gain levels. Here it was seen that many ineffec- 
tive controller gains could be el imin ated, with no discernible performance 
degradation. The result was a further simplified controller. 

Now we turn to two further considerations for the sensor complex. The 
first is sensor /equalization economy, which refers to the desire for mini- 
mum compensation of the vehicle dynamics and a minimum array of sensor and 
equalization computation to achieve performance goals. This desire arises 
out of considerations such as cost, reliability, maintenance, and the like 
Referring back to Fig. 31, note that certain small gains in the Ap and Bp 
matrices have been set to zero. This is based on control realities, but 
is also consistent with economy of equalization. The importance of re- 
evaluating the stability, response and sensitivity of this further simpli- 
fied system has been emphasized because a design which does not meet per- 
formance standards is useless, regardless of how simple it may be. 

Examination of certain eigenvalue sensitivities can divulge the 
effects of deleting some controller gains upon the closed-loop root loca- 
tions. This was discussed previously regarding first-order parameter sen- 
sities of the eigenvalues to changes in the F and G matrices. Here we use 
first-order eigenvalue sensitivities to changes in the regulator matrix C. 
An example of the matrix sensitivities for the short-period mode eigen- 
value is given in Table 50 . The only significant sensitivities here indi- 
cate that the important controller gains are ; 
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v > 5 b 5 b 


w, 5c, / dt 5 q 

Two observations can be made. First, we presumably could have picked in 
advance which sensitivities would be important based on the controller 
gains that were eliminated. We see that their presence has negligible 
effect on the location of the short-period mode, based on the sensitivity 
magnitudes of the corresponding deleted terms in the controller gain 
matrices. 

The search for a practical optimal controller, which has guided our 
synthesis and assessment throughout, leads to one other goal for the 
design. This is the desire for satisfactory performance over a wide range 
of operating conditions. Suppose, for example, that it was necessary to 
store a redundant set of controller gains for ten different helicopter 
operating setpoints. The storage problems alone prevent such a design 
From ever being used in a real-world situation, not to mention the huge 
expense involved in assessing the performance of all ten designs. What 
we would prefer is a single design that is easily adapted to almost all 
operating modes. 

The concerns here are guided by questions like : 

Suppose we wish to create a command attitude (0 C ) system 
from a rate -of- climb command (h c ) design simply by open- 
ing the h loops. Does the resulting e c system exhibit 
satisfactory stability, response and sensitivity charac- 
teristics, and if not, how can they be improved? 

It is obvious that the desire for commonality of elements/gain settings 
is closely related to the desire for graceful degradation. Thus, the sce- 
nario posed above can be examined in parallel with the behavior of the 
9 loop when the h loop has degraded. These kinds of considerations insure 
that the optimal controller that is ultimately implemented is as simple 
and universal as possible. 
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F. 


SUMMARY 


In summary, we return to our original goals in assessing the design 
to see whether we have been successful. 

• D etermine whether design goals have been met . The primary de- 
sign goals in this work are the variable bandwidths that are set using 
the regulator cost function. An examination of the frequency response 
plots for the closed-loop system and various open-loop transfer functions 
assured that these bandwidths were attained. In addition, a good stretch 
of —20 dB /decade slope on either side of the open-loop frequency response 
crossover points was also created. This allows the bandwidth to be ap- 
proximated by the crossover frequency. 

Many other secondary design goals were also attained. These include: 
a filter -observer that has an rms state estimation error within the error 
budget^ zero steady-state response of the two primary error variables, 
accomplished by including the integrals of these variables in the regu- 
lator cost function^ and a controller that provides control decoupling 
for the longitudinal degrees of freedom. Each design goal was produced 
by a specific synthesis step (e.g., computing the regulator cost function 
to obtain the desired bandwidths) and the assessment was used to verify 
the desired results (and sometimes to guide iterations). 

• Define quantitatively all the important properties of the system . 
Three examples were cited to illustrate this goal. The first was to expand 
the scope of the results. When design synthesis was completed, the end 
product was a set of four controller gain matrices. It was practically 
impossible to predict system and controller properties and performances 
based only on looking at this rather compact collection of gains. Instead, 
the scope of the design was expanded to include frequency responses, tran- 
sient responses to commands and disturbances, root loci, various stability 
criteria, sensitivity criteria and many statistical metrics. All were 
essential in obtaining a quantitative design assessment. 

The second example given for this assessment goal was to highlight the 
system's dominant properties. In both the stability evaluation and the 
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response evaluation we were able to identify, for each variable, those 
modes that dominated the particular characteristics observed. In doing 
so, a fairly complicated transfer function was invariably reduced to a 
much simpler and more manageable one, making the stability and response 
behavior easier to comprehend. Often these dominant modes provided a 
convenient method for partitioning the response. The partitions corre- 
spond to certain well-known systems (such as zero position error, zero 
velocity error, etc.) and the behavior is thus more readily analyzed. 

A third example was to reveal the sensitivity of dominant properties 
to uncertainties. Having identified the dominant modes in connection 
with stability and response characteristics, we were then able to predict 
the effect upon system behavior of changes in these dominant modes. These 
"changes" ran the gamut from off-nominal conditions, controlled element 
equation uncertainies and simplifications, through thresholds, hysteresis 
and other parasitic nonlinearities, to component failure. The goal, of 
course, was to assure satisfactory FCS performance in the face of any or 
all such uncertainties. The sensitivity is variously quantified using 
first-order gain and parameter sensitivities, describing function analysis, 
and failure simulation. 

6 Broaden the understanding of the optimal system design results . 

This goal guided all of the techniques used. Too often, an optimal con- 
troller design is only a theoretical solution because the design results 
are not well understood. In the introduction we gave three examples of 
the breadth of understanding desired. The relative importance of the 
various controller partitions was illustrated In the controlled gain level 
considerations. On the basis of control saturation and minimum increment 
of control we were able to define control regions for each output varia- 
ble. The effectiveness of each control point in regulating each output 
quantity was also evaluated. 

Many aspects of the controller behavior could be attributed to speci- 
fic factors in the controlled element as a result of the assessment. For 
instance, the steady-state error response characteristics for the primary 
(or outer) loops were established. Dominant modes in various responses 
were identified. Conditions under which control poles are independent of 
the regulator solution were identified. 
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Implications for future design strategy appeared at each assessment 
step. In the synthesis of the controller little attention "was paid to 
questions of graceful degradation and commonality of elements. From our 
assessment of these design aspects it is apparent that more considera- 
tion should be given to systematic analysis of these features during 
design synthesis. 

• Explore simpli fication of the design •which may lead to a simpler, 
more practical controller . Throughout the assessment we have been con- 
cerned with design simplification. This began with identification of the 
dominant system properties and their sensitivity to uncertainty. Simpli- 
city was also paramount in considerations of controller gain levels, sen- 
sor array, equalization economy and commonality of elements. Although 
simplicity was continually emphasized, it was also stressed that the 
simplified systems should be rigorously analyzed to insure that the sim- 
plification was not accompanied by a loss in performance. 

• I nvestigate the implications of controller architecture on results 
Such an investigation was conducted during assessment of graceful degra- 
dation. This particular aspect of the evaluation was dependent on the 
controller architecture. We examined only two candidates — a single 
thread, multiply -redundant controller and an adaptive filter -observer 
controller — and saw that these two mechanization strategies required 
quite different assessment techniques. The resulting performance impli- 
cations were also different, as well as the implementation problems. From 
just two example architectures it was seen that the controller properties 
can be quite implementation-dependent. 

• Ge neralized design for other applications . The imposing task that 
characterizes the assessment of a new FCS design makes it desirable to 
design a system that is useful over many vehicles and many operating modes 
This aspect was touched upon when we looked at the suitability of a given 
design in several operating modes. Thus, it was desirable to permit an 
h c system to become a 8 C system simply by opening the h loops. We could 
also examine the performance of a candidate FCS design in various aircraft 
perhaps in a simulation environment. 
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So we see that the list of assessment tasks enumerated in Table 4k 
has allowed us to attain our assessment goals. Figure 32 depicts a flow 
chart that could he used in design assessment. If it is ever necessary 
to interrupt this flow, due to the need for redesign, the chart should 
he re-entered at the top. This is to insure that the design always meets 
all assessment goals. 

In final conclusion, it cannot he emphasized too strongly that the 
fundamental purpose of design assessment is to provide an in-depth under- 
standing of the candidate system properties and behavior. This enhanced 
comprehension is vital in an optimal controller design^ without it, such 
designs will he overly complex and the total system behavior will he 
poorly appreciated. 
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SYSTEM PROPERTIES 


Stability 


Response 


Stable closed-loop 
roots 


Adequate gain and 
phase margins 




Primary Variables 


Rapid, well -damped 
time histories 


Good output/command 
fol lowina 


Disturbance suppres- 
sion , 


Dominant Primary 
Modes 


Secondary Variables 

Motion harmony 

♦ 

Favorable couplinq 

I 

Disturbance sup- 
pression 


Dominant Modes of 
Secondary Variables 



Reasonable author- 
ity 1 imits 

* 

Motion harmony 

♦ 

Favorable/ 
negligible cross- 
coupling 


Key Response 
Modes 


Parasitic Non- 
1 i near i ties 



Relative insensi- 
tive to off-nominal 
conditions 


Design independent 
of threshold effects 


Component failure 
results in stable 
system 


CONTROLLER PROPERTIES 


Gain Levels 


Sensor Complex 


Good dynamic range 


Minimum sensors 


Reasonable sensor threshold 

♦ 

Ineffective pathways eliminated 


Equalization economy 

t 

Cormonalitv of elements 


Figure 32. Assessment Flow Diagram 
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APPENDIX A 


TRANSFORMATIONS FOR SINGULAR FILTER PROBLEMS 
ENABLING EIGENVECTOR DECOMPOSITION 

This appendix documents a transformation of the Euler-Lagrange equa- 
tions for the optimal state estimation or "filter" problem. This trans- 
formation is motivated by a desire to solve singular filter problems 
wherein some or all of the measurements are noise free. Attention is 
restricted to the steady-state solution for constant plants. Solution 
without transformation via existing computer codes (e.g., Ref. A-1 as modi- 
fied by R. Bach and G. Slater) is not possible because the inverse of the 
measurement noise spectral density matrix is required. Transformation 
for continued use of these codes is recommended since the codes make use 
of the robust and efficient OR algorithm for eigenanalysis . Furthermore, 
this approach satisfies our research needs for a minimum of additional 
software development and risk of encountering numerical difficulties. 

Consider the plant and measurement equations 


where 


x = Fx + Gu + rw 


x(0) = x Q 

z-j = H-|X + v-j 



z 2 = h 2 x + v 2 




Vector 

Dimension 

Vector Name 

X 

n 

State 

u 

P 

Control 

w 

q. 

Process noise 

Z 1 

mi 

Noisy measurements 

Z 2 

m 2 

Noise-free measurements 

Vi 

m-i 

Measurement noise 

v 2 

m 2 

Dummy measurement noise 
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G, r, H-| and H 2 matrices are assumed to have full rank. Process and meas- 
urement noise are Gaussian with 

E[w] = 0 E[v-,] = 0 E[v 2 ] = 0 

E[w(t 1 )w'(t 2 )] = QS(t-j - t 2 ) 

E[v-| (t-, )v' (t 2 )] = R-j S( t-] - t 2 ) 

E[v 2 (t-,)v 2 (t 2 )] = Q eRgSCt-, - t 0 ) 

w, v-| and v 2 are mutually statistically independent random processes. 

Their power spectral densities, Q, R-j and R 2 , are diagonal matrices having 
full rank. 

It is always possible to order the measurements so that noisy measure- 
ments are in z-j and noise-free measurements are in z 2 as reflected in the 
above equations. 

Also, it is always possible to reorder the state vector elements by 
means of elementary transformations to obtain H 2 = [H 2 i H 22 ] wherein H 22 
full rank. 

Furthermore, the transformation matrix, T, defined in Table A-1 can 
always be used to obtain H 2 in the form: 

H 2 = [H 21 H 22 ] = [0 I] 

Next, consider r in partitioned form wherein r 2 has rows. 



Form A 1 = P 2 Qr 2 - Let the rank of A 1 be defined as r. Assume for the 
present that A -1 has full rank, r = nip. This assumption is equivalent 
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TABLE A-1 . TRANSFORMATION MATRIX 


Given : 


x = F x + G u + r w 
Zj = H-j x + v -j 
x + v 0 


x(o) = x 0 


z 2 ~ **2 
**2 = [H 2 i Hgg] 

and. H22 is m 2 m 2 nonsingular . The nonsingular transformation, T, 
x = r x where 


results in 



' I 0 1 


I 

0 

T = 

ffil 

ro 

ro 

1 . 

, T 1 = 

_-h 2 - 2 1 h 2 i 

H 2 2 _ 

) . 


1 



1 ' 

, Z 2 ) 

( **21 X 1 + H 22^2 ) 






F = TFT 


[Fll Fi 2 l 

L F 21 F 22 J 


Fll 

F 2 i 

F 2 2 


F 1 1 - F 12 h 22 1h 21 


Fl2 = F 12 H 22 


**22 F 21 — H 2 2 F 22 **22 **2 1 + **21 F 11 — **21 F 21 **22 **21 


**22 F 22 **22^ + **21 f 12**22^ 


G = TG = 


T = TT = 


[:;] • 


Gl = Gl 


F i = r ^ 


Gq = **2i Gi - 


r 2 “ **21 r 1 


H = HT 1 

"Hn 

**12 " 




H21 

H22J 



Hn = **n - 

_ _-i_ 

**12 1*22 **21 


7 

**12 = H12 **22 

H21 = 0 



9 

**22 = 1 
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**22 **2 
**22 F 2 



ii i mi 


mil hi 


i i 


to the requirement that the derivative of every noise-free measurement 
contain an independent white noise component. All subsequent develop- 
ments either require this assumption he satisfied, or that a single para- 
meter relaxation of this assumption he used (refer to Appendix B) . 

Euler-Lagrange equations for the filter problem are (e.g., Ref. A- 2 , 

pp. 395-396): 



These equations may he partitioned more finely and Laplace transformed. 


(si-f-h) 

" F 12 

r q QTi 

-f 2 i 

(si — F 22 ) 

r 2 or q 

H'l'jR'j H-j 1 

» -1 

H-j -j R 'j H*| g 

(sI + F qi ) 

1 Hl 1 

eH] 2 Ri ^H-| 2 

eF 12 
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The last Euler -Lagrange equation has been multiplied through by the para- 
meter e- If e is allowed to approach zero from above, the result is that 
z 2 = x 2 since R 2 is diagonal and nonsingular by assumption. z 2 = x 2 re- 
places the last Euler-Lagrange equation. 

Let us examine this result from another viewpoint. Assume a solution 
for the optimal state estimate, x, of the form (e.g.. Ref. A-2, page 396) 



However, in order that the error in a log-likelihood cost function remain 
bounded for the noise-free measurements [recall these are weighted by 
(l/e)R 2 ^ in that cost function], z 2 = x 2 = x 2 is required. Consider the 
assumed solution. Since X-j and X 2 are not necessarily zero, P 12 and P 22 
must be zero because P 2 2* the error covariance for the error in the esti- 
mate of x 2 , is zero. 

The second Euler-Lagrange equation may be solved for X 2 . Making use 
of z 2 = x 2 and A 1 = r 2 QT 2j > the result is: 

X2 = A[F 2 iX-i — (si — F 22 )z 2 + z 2 q — r 2 QTlll + G2 u ^ 

The first and third Euler-Lagrange equations can be rewritten as 


(si— f -] 1 + r 1 Qr 2 AF 21 ) 

(rw! -r-iQT2Ar2 ( ri) 

(-) 

( H 1 1 R 1 1 + F 2lAF21 ) 

(si +F^-| — F 2 iAr 2 QrT ) 

U) 


! (x 10 -r 1 'T2Az 2 o) 
(Mo -F 21 Az 20) 


(G-j -riQT2 A G2) 
(-F 21 AG 2 ) 


0 (f-| 2 + r-) qr 2 A[ si — f 22 ] ) 

(Hi 1 ) (— i^i 1 Hi 2 + F 2 -|A[ sl — F 22^ ) 



Specialization of the assumed solution for the optimal state estimate 

x-j = x + PX-) 
x 2 = x 2 = z 2 

was justified previously. This assumed solution may he substituted in 
the first and third Euler-Lagrange equations to obtain them in terms of 
P, x-j and X-j . 


(sI-Fn + r QT2-AF2 1 + P[HtiRi 1 H-|i -t-F^iAFgi]) 


(HiiR^Hn +F21AF21) 


(r-|Qr-| — riQr^A^Qri + t F i i — r-iQT2AF2i]P 
+ P[F^ 1 — F^-jA^OTi 3 -PCH^r/h-ji + F2-jAF2t ]P 

(sI + Fii — F 2iAT2QTi — [h] iRi^Hi 1 + F2iAF2i3P 



( ( x 10 — Ft GF2 Az 20 + P[^10 “ f 21 Az 2o3) N \ 

- 

(F12 + Ti QT2 a ( s I — f 22) + FF2iA(sI — F 22) 



\ (^10~ f 21 Az 20) / 

( f 21 A ( s I —f 2 2 ) ~ h 11 r T 1h 12) 



(g 1 — r -j Qr^AG^ — PF 21 AG 2 ) 


(ph^Rt 1 ) 

+ 


u + 



(-f^ag 2 ) 


(h^rT 1 ) 


Elimination of the costate-to-state coupling in the above equation results 
in an algebraic Riccati equation. 

riOTi -riqrpArpQn + [Fn -norg^ilP 


+ P[F* 1 — F 2 lAT 2 , T'l 1 ” P[Hi 1R1 ^Hj j + F21AF2I ]P 


0 



Assume the Riccati equation is satisfied, and rewrite the x-j equation in 
terms of the filter gain matrices 


Kn = EH^E, 1 

and 

k 12 = ( pi 21 + IW2) A 

This results in 

(si — F i -j + K-| 2 F 2 1 K-| i Hi i )x-j 

= © ~ 2 Zg©) K-j -j z -| 4- "”^11^12^ z 2 (^1 ~ ’ K-j 2 ^ 2 ^ ^ 

Occurrence of [s(x^ — K-| 2 Z 2 ) ~ (^1 0 “ K 12 z 20^ in atov ' e equation suggests 
a more appropriate choice of filter state coordinates is the (n— m 2 ) vec- 
tor, y. 

y = x-, - K 12 z 2 

Then the differential equation for the estimate becomes; 

y = (F-| •] — K12F — K-j iHt i )y + ( [F-] i — K12F21 — K-j iH-] •] ]K-j 2 

+ F 1 2 ” Ki 2-^22 — ^1 1 Hi 2 ) z 2 + ^- 11 z 1 + (^1 — 
y(0) = x 10 -K-12 z 20 


x-| is obtained from y by means of 

x 1 = y + K -| 2 z 2 

Figure A-1 shows the structure of the resulting reduced order filter. Re 
call, however, that this solution is in terms of a transformation of the 
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Figure A-1 . Optimal Filter-Observer Structure (r20r2 Full Hank) 






state vector for the original problem. Appropriate inverse transforma- 
tion must he applied before interfacing this solution with the regulator 
portion of the optimal stochastic controller. 

Consider the error in the estimate, y. The quantity corresponding to 
y is y where.- 


y = x-j - K 12 z 2 


X 1 ~ K 12 x 2 


The differential equation for y is; 

y = (Fn-Ki^i-K^Hi^y 

+ ( [F-j i K-| 2 F 2 -j — K-| -[H-| i ]K-j2 + F "j 2 — 2^22 ) ^2 

+ K'l'jH'i^x-] + (Gi — K-\^G 2 )m + (pi — K , \^' 2 )'W 

y(o) = x 10 -k 12 z 20 

The differential equation for the error in the estimate, y, is: 

y = y-y = (Fn -K i 2 f 2 i -K iiHn)y 

— (ri -K 12 r 2 )w + K-|iV-| 

A 

The error in y is also the error in x-j because 

y = (x-, — K -[2 z 2 ) - ( Xl -K 12 z 2 ) = x-| 

Consider the properties of the covariance of x-| . The covariance of 
x-| is denoted by X. X is the solution of the equation: 
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(Fn — ^■i2^ , 2i K n H n) x + x(Fi 1 .-f 21 k 12 -h 11 Ki 1 ) 


+ (r-j -K 12 r 2 )Q(r 1 -r 2 K 12 ) + k-j -j Ri Ki -] = o 


Necessary and sufficient conditions for a minimum of the trace of X with 
respect to K-jp and K 22 are 

K* q = XH^rY 1 

K 12 = (XF 2 -j + Tv^) (r 2 QT 2 ) 

provided R-| and (r 2 qr 2 ) are positive definite. R-j and (r 2 QT 2 ) are always 
positive definite here hy definition and assumption. 

Substitution of K* q and K* 2 into the covariance equation results in the 
Riccati equation obtained earlier with P replaced by X. This confirms the 
optimality of the singular filter solution and shows that P is the covari- 
ance for x-| and y. 

Several concluding comments are in order. 

The solution given here is a special case of the solution given in 
Ref. A-3, pages 327-332. It is a special case in that here we have trans- 
formed the state vector at the outset in order to avoid introducing (n— m 2 )^ 
parameters unnecessarily. (Those parameters are in the first (n~m 2 ) columns 
of the T matrix in Ref. A~3> page 329- ) 

In view of the above, one might ask, "Why bother with this development 
if the results of Ref. A -3 may be specialized?" The answer has several 
points : 

• The development in Ref. A-3 is a "heuristic" one, 
particularly with respect no incorporating the 
"free gain matrix, B-] . " 

O There is no demonstration that the Riccati equation 
in Ref. 2 results from the Euler- Lagrange equations 
or their equivalent. This is necessary in order that: 
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— The resulting filter he optimal. 

— The eigenvalue decomposition method of 
solution he applicable. 

• There is a desire to eliminate the unnecessary para- 
meters in the Ref. A- 3 formulation for the sake of 
simplicity. 
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APPENDIX B 


UH-1H EQUATIONS OP MOTION 
AND ACTUATION SYSTEM 


The first quantitative step in design is formulation of the equations 
of motion in state vector form: 

x = Ex + Gu + rw 

z = Hx + v 

y = HpX 


where 

x is the state vector, composed of shaping 
filter states (x s ) and controlled element 
states (x c ) 

u is the control vector 

w is the process noise vector 

z is the measurement vector 

v is the measurement noise vector 

y is the output vector, composed of the 
quantities being measured (Hx) and the 
other output quantities (y 0 ) 

The specific selection of the vector components and their corresponding 
equations (kinematic, perturbation, shaping filter and auxiliary) is the 
subject of the next two subsections. The longitudinal vehicle dynamics 
are given first, then the dynamics of the shaping filters and, finally, 
the actuator model. 

LONGITUDINAL EQUATIONS OF MOTION 

For the longitudinal problem, the vector components indicated in 
Table B-1 are used. The equations of motion can be stated in terms of 
the variables comprising those vectors. 
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TABLE B-1. VECTOR COMPONENTS FOR LONGITUDINAL PROBLEM 


Vector 

Elements 

Computer 

Mnemonic 

Description 

X S 

u g 

xoi a 

longitudinal gust velocity (ft/sec) 


w g 

X02 

Vertical gust velocity (ft/sec) 



X03 

Commanded altitude rate (ft/sec) 

*c 

u 

X06 

Longitudinal velocity (ft/sec) 


w 

xoT 

Vertical velocity (ft/sec) 


q 

X08 

Pitching velocity (rad/sec) 


3 

X09 

Pitch angle (rad) 


S B 

X10 

Longitudinal cyclic (in.) 


Sc 

XII 

Main rotor collective (in.) 


/ h e dt 

X12 

Integral of altitude rate error (ft) 


X BAR 

X13 

Horizontal stabilizer bar (ft) 


/ u ASg dt 

Xl4 

Integral of airspeed error (ft) 3 * 


/ x dt 

XI h 

Integral of inertial speed error (ft) c 

u 

&B 

DBD 

Longitudinal cyclic rate (in. /sec) 



DCD 

Main rotor collective rate (in. /sec) 

w 

Wug 

PNU 

Longitudinal gust velocity process 
noise (ft/sec) 


Ww s 

PNW 

Vertical gust velocity process noise 
(ft/sec) 


V h c 

PNC 

Commanded altitude rate process noise 
(ft/sec) 

z 

h 

HD 

Altitude rate (ft/sec) 


q 

Q 

Pitching velocity (rad/sec) 


e 

TH 

Pitch angle (rad) 


5b 

DB 

Longitudinal cyclic (in.) 


s c 

DC 

Main rotor collective (in.) 


/ h e dt 

HDI 

Integral of altitude rate error (ft) 


fee 

KDE 

Altitude rate error (ft/sec) 


/ u AS e dt 

ASI 

Integral of airspeed error (ft) 13 


u AS e 

ASE 

Airspeed error (ft/sec) 33 


/ x dt 

XDI 

Integral of inertial speed error (ft) c 


X 

XD 

Inertial speed error (ft/sec) c 


ax 

AXP 

Longitudinal acceleration independent 
of pitch (a x — g sin P) (ft/sec2) d 


a z 

AZP 

Vertical acceleration measured £ x ft 
forward of c.g. (ft/sec 2 )* 3 - 

y 0 

a 

AOA 

Angle of attack (rad) 


X 

XD 

Inertial speed (ft/sec) 35 


u AS e 

ASE 

Airspeed error ( ft/sec ) c 


a State vector mnemonics are automatically generated by software. 

^Used in 100 kt case. c Used in 1 kt (hover) case. d Not used in all cases. 
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Perturbation Equations of Motion 


These equations are derived in Ref. B-1, with the exception of the Xg^ 
equation, derived in Fig. B-1 . 

u = (“X u ) u g + (-Xw)^g + X u u + X w w + (X a -W 0 )q 
+ (— g cos 0 O )9 + XggSg + 

w = (~Z u )u g + (-Z w )wg + Z u u + 7 w w + (Zq + U Q )q 

+ (-g sin 0 o )e + Z 5 b 5 B + Z- 5 C & C + ^5b x BAR 

q = (-M u )u g + (-Mfcr)wg + M u u + MwW + Mqq 
+ M6b s B + m 5c 5 C + M &B X BAR 

0 = q x BAR = -0.333 x BAR + 4.97 q 



Figure B-1 . Pitch Axis Stabilizer Bar (from Ref. B-2) 
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Kinematic Equations 


These .equations are derived in Ref. B-1 . 

(/ h e dt) = h c + (-sin 0 o )u + (cos 0 o )-w + (-Vt 0 )9 + (-e)/h e dt 

(/u ASE dt) = (-cos 0 o )u g + (-sin 0 o )w + (cos 9 0 )u 

+ (sin 9o) v + (~e)f U ASE ^ 

( / x dt ) = (cos 0 o )u + (sin 0 o )w + (-e)/xdt 

The e factors are used to obtain bounded rms integral responses when inte- 
gral loops are not yet closed. Otherwise, e = 0. 

Shaping Filter Equations 

These equations are derived using a simplified version of the model 
for random turbulence. Gradient effects associated with the normal turbu- 
lence component are neglected. The random turbulence components have 
Gaussian probability density functions with zero means. The standard 
deviation should be chosen from a Rayleigh probability density func- 

O 

tion* having a characteristic speed of ob^ ft/sec. However, for the sake 
of simplicity, the mean value of aug* which is ob^, is used. 

a u = a CT = 2.79 - 0.2^5 log 10 h ft/sec h > 100 ft 
& 

= 2.3 ft /sec h <. 100 ft 


The standard deviation cfwg is a function of Oug* The frequency content 
of the random turbulence and are functions of altitude. 


The power spectral densities for the longitudinal and normal random 
turbulence components at a given altitude are respectively: 


*The Rayleigh probability density function is for Ug. 



where 


°u 

- g o 

Ug oi 2 + (V Ao /La ) 2 

gg^d-s^vAo/jy) 

$W S cc? + (1.59^/Lw) 2 


o 2 



<J) dm 


<IVg is a lower-order approximation to the Dryden power spectral density 
such that the mean- square level and half -power frequency are preserved. 

The differential equations for unit-white-noise shaping filters pro- 
ducing output variables Ug and Wg having power spectral densities and 
^Vg respectively are; 

Ug = - | V Ao | /Lu u g + crug \/2 |Va 0 |/Lu Wug 


w g = - 1 . 59^ |v Ao I A wWg + (% g ^2(1 .^9i)|VA 0 | Aw 


S 


where w u and w w are independent, zero-mean, unit white noises. Va q i s 
the trim airspeed. 

The integral scale lengths and 1^ are given as functions of alti- 
tude h by 


= 1 h-5[ 1 oo] 1 = 673 


1 00 < h < 1 750 ft 
h < 100 ft 


iw - h 


h < 1750 ft 
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The standard, deviation for the normal turbulence component tJvrg is related 
to the standard deviation for the longitudinal turbulence component crug 
through the integral scale lengths. 


c% g 


^/Lw/Lu. hUg 


The rate -of -climb command equation is given by: 


(h c ) = h c h c + ffh c N/ 2 aj h c ™h c 

where w^ c is an independent, zero-mean, unit white noise. 

Auxiliary Equations 

These are derived in Ref. B-l . 

h = (sin e 0 )u + (—cos 0 o )w + Vt o 0 

h e = h c + (-sin 0 o )u + (cos 0 o )w + (-Vij o )0 

u AS e = (-cos 0 o )u g + (-sin 0 o )wg + (cos 0 o )u + (sin 0 o )w 

x = (cos 0 o )u + (sin 0 o )w 

a x = (— X u )ug + (— X w )wg + X u u + X-wW + Xqq 

+ (— g cos 0 O )S + Xggbg + X§q 5 q + ^ 5 b x BAR 

a-z ~ — (Zu“ - 0 xMu) u g ~ ( Zw — ix^w^g + ( Zu _ -£ x ^u) u g 

+ (Zyr ~ i x M w) w + (Zq-ix^q)^ + 

+ (Zs c — i X M5c)5 c + ( Z§g — ^ x M6b ) X BAR 

a = (-l/V To )w g + (l/V To )w 
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Table B-2 presents the longitudinal matrix equations in literal .form. 
Table B- 3 lists values for the parameters found in the longitudinal state 
vector equations, for two key flight conditions. All parameter values 
axe from Ref. B-2 except those for the shaping filter equations. These 
latter were derived above. 

ACTUATION SYSTEM 

The actuation system has identical dynamics for all four control axes. 

A first-order approximation to the third-order combination of series and 
boost actuators is used. The first-order model includes an approximation 
of the series actuator phase lag at low frequencies. Dynamics in the 
parallel trim actuator path are neglected. The bases for these approxi- 
mations are explained elsewhere in this report. 

Let DX represent the output displacement of any boost actuator (DX = DB, 
DC) and DXD represent d(DX)/dt. The differential equation included in the 
plant model for DX is: 


DX = DX 

Unity gain feedback around the integration element of the actuator model 
is provided by the optimal regulator solution. The Qr and Rp cost func- 
tion weighting coefficients are selected to obtain each effective actuator 
break frequency, co^, which approximates the combination of the series and 
boost actuator dynamics. The method for selecting the cost function coef- 
ficients to accomplish this is described in Sections II and IV. 

^A = 1 / to boost + ^series/^series^ = 25*73 rad/sec 

“>boost = 50 . rad/sec 

series = ^ ( From Ref * B ~3) 

•^series = 75- rad/sec 
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TABLE B-2. LITERAL MATRICES FOR LONGITUDINAL INPUT DATA 
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a XM used in hover example; 
ASI used in 100 kt example 


TABLE B-2. (Continued) 



(concluded on following page) 


TABLE B-2. (Concluded) 
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a XDI used in hover example ; ASI used in 100 kt cruise. 

t>XD used as measurement in hover, output at 100 kt,- ASE used as measurement at 100 kt, output in hover. 
c Used to investigate acceleration measurements in 100 kt case. 



TABLE B-3. VALUES FOR LONGITUDINAL EQUATIONS 
OF MOTION PARAMETERS 


Parameter 


YA 

kt 

<J 

t-3 

0 

II 

& 

o 

ft /sec 

Lu 

ft 

Lw 

ft 


rad/sec 

a h c 

ft/sec 

a Ug 

ft /sec 

CTwg 

ft/sec 

ho 

ft 

Wo 

ft /sec 

u 0 

ft/sec 

S 

ft /sec 2 

So 

deg 

X u 

1 /sec 

Zu 

1 /sec 

Mu 

1 /sec-ft 

X-w 

1 /sec 

Zw 

1 /sec 

1% 

1 /sec-ft 

x q 

ft/sec 

Z q 

ft /sec 

Mq 

1 /sec 

X&B 

ft/in. -sec 2 

z 5b 

ft /in. -sec 2 

MSB 

1 /in. -sec 2 

X&c 

ft/in. -sec 2 

Z&c 

ft/in. -sec 2 

MS C 

1 /in. -sec 2 

-^x 

ft 

h 

ft 


100 kt 

1 kt 

(Case 128) 

(Case 122) 

100.0 

1.0 

168.78 

1.69 

1^50. o 

1450.0 

1000.0 

1000.0 

0. 1 

0. 1 

1 .0 

1 .0 

2.06 

2.06 

1.71 

1.71 

1000.0 

1000.0 

8.91 

0. 12 

168.54 

1.68 

32.2 

32.2 

3.02 

4 .o 4 

-0.0451 

-0.0034 

0.0888 

-0.0991 

0.0050 

0.0019 

0.0950 

0.0250 

-0.9963 

-0.3850 

-0.0066 

-0.0038 

1 .7727 

0.5797 

-3.4493 

0.2913 

-0.7012 

-0. 1900 

0.5568 

1 .o 4 o 6 

5.0738 

0.3214 

-0. 1728 

-0.1691 

1 .1402 

0 . 6806 

-13. 1274 

- 9.7745 

—0.0223 

-0 . 0033 

7.30 

7-30 

1000.0 

1000.0 


Refer- 

ence 
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CONTROL AUTHORITY MDDEL 


The model for determining control displacement and rate requirements 
assumes mechanization in terms of series, parallel and boost actuators as 
shown below. Mechanization is the same for all control axes. Dynamics 
of the series and parallel servos are neglected. This block diagram 
enables us to write equations for the series actuator displacement 

XS = -KXS + (1 - K/o> A ) DXD 
DXS = XS + DXD/cjo a 

and the parallel actuator rate 

DXPD = K DXS 

where K is the gain in the integral path determined in the course of design. 



Parallel 
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APPENDIX C 


UH-1H, HOVER, EXAMPLE APPLICATION 

This appendix presents the highlights of the synthesis and assessment 
procedures used in the design of a longitudinal flight control system for 
the UH-1H at hover. None of the design iteration steps are included, 
only the final results. Also, just a sampling of the assessment data is 
presented. The next several paragraphs describe the data which follow. 

Table C-1 provides an overview of the primary design goals in this 
synthesis. They involve identifying the flight condition and the corre- 
sponding functions to be performed by the FCS, which in turn specify the 
bandwidths. The actuator bandwidths were computed using an approximation 
to the actual UH-1H hardware actuators (see Appendix B). 

Tables C-2 and C-3 list the matrices used to define the problem. The 
shaping filter and controlled element equations are provided in Table C-2, 
while the measurement and output equations can be found in Table C-3- 

Table C~4 contains the elements of the final filter -observer synthe- 
sis, discussed previously in Section IV (see Table 16). 

•Table C-5 begins the regulator synthesis portion of this example with 
a list of the relevant controlled element transfer functions. Figure C-1 
presents an example of the portion of the regulator synthesis procedure 
that is concerned with making approximations to the transfer functions, 
computing the associated cost function weightings for each surrogate con- 
trol, and choosing the most effective control point by selecting the 
smallest cost function weighting. The 0 closure is computed here, using 
cyclic and collective as the candidate control points. Notice that for 
the TH/DB transfer function both the asymptotic and actual Bode plots are 
used in computing K and n. Compared to the cost function weighting when 
collective is used as the control point there is little difference between 
the asymptote and the actual transfer function approximations . Based on 
the cost function weightings, longitudinal cyclic is chosen as the most 
effective 0 control point. 
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As was mentioned in Section IV, it is often possible in a particular 
application to predict the effect of various closures on subsequent 
closed-loop transfer functions and thereby eliminate the iterative sec- 
tion of the regulator synthesis by effectively closing all of the loops 
simultaneously. This was possible for the example at hand. Table C-6 
presents the design of the regulator by indicating the process of choos- 
ing the most effective control point at each step. Notice that the cost 
function weighting on collective rate, is an order of magnitude 

larger than the cost function weighting on cyclic rate, Rjjbd* This 
reflects the results of the 100 kt design study wherein collective acti- 
vity was judged to be too severe when the cost function weightings on the 
actuator rates were equal (see Section IV). 

Table C-7 presents the final regulator design for this example. 

Once the filter -observer and regulator syntheses are complete, all 
that remains in the optimal controller design is to compute the controller 
coefficient matrices, shown in Table C-8, and the controller transfer func- 
tions, shown in Table C-9. 

The design assessment begins in Fig. C-2 with closed-loop transfer 
functions and transient responses to step inputs of Ug, Wg, and h c . This 
provides information on the gust suppression characteristics of the design, 
as well as the response to command inputs. 

Next we examine Bode plots of various closed-loop transfer functions 
to see if the design goals were attained. Figure C-3 shows the rate-of- 
cliirib response to rate-of-climb command^ the —^5 deg phase angle falls 
at approximately 1.2 rad/sec. Figure C-4 presents the error response for 
the same variable — the —3 dB point is at 1.0 rad/sec. 

Figure C-5 is a Bode plot for the approximation of the e/e c transfer 
function with all loops closed. Figures C-6 and C-7 plot the propor- 
tional plus integral and proportional plus integral error transfer func- 
tions, respectively, for the groundspeed loop. Again we see a close 
agreement between the design goals and the actual design. 

The final steps presented here for assessing the design via classi- 
cal control metrics make use of the 1^ matrix (see Fig. 6) to compute 
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various open-loop transfer functions. These open-loop transfer functions 
provide a means for determining the effect of the outer-loop closures on 
the inner-loop bandwidths. 

Figure C-8 is a plot of the rate -of- climb proportional-plus-integral 
transfer function with the proportional and integral loops opened. The 
crossover frequency is about 1.2 rad/sec. Figure C-9 examines the rate- 
of-climb integral bandwidth by plotting the HDl/VHI transfer function 
with the HD I loop opened. Figures C-10 and C-11 show similar plots for 
the proportional -plus-integral groundspeed loops. Figure C-12 plots a 
more inner loop, the 0 loop, with the 0, x, and fxdt loops opened. Here 
we see a crossover frequency of approximately 2.3 rad/sec. 

Finally, an assessment of this FCS design based on statistical met- 
rics is presented in Table C-10. The process noise augmentation has been 
removed here. We can compare the rms response of the plant states and 
filter states to the combined process and measurement noise with any 
secondary design goals specified at the outset. We also have at hand 
the rms response of the outputs and controls, which enables us to compute 
measures of control activity such as the positive-going zero crossing rate, 
also shown in the table. 

This appendix has presented just a sample of the design synthesis and 
assessment tools available. Increased familiarity with the design pro- 
cess will enable a designer to feel comfortable with this method and the 
resulting data analysis. 
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TABLE C-l 


EXAMPLE APPLICATION 


UH-1H Longitudinal Flight Control System 


Flight Condition; Hover 


Functions ; 


Rate -of -climb command 


Groundspeed hold 


Bandwidths ; 

(rad/sec) 

Cyclic (DB) ; 

25-73 

Collective (DC) : 

25.73 

Pitch (TH) : 

2.0 

Rate-of-climb error (HDl) : 

1 .0 

Integral of HDE (HDl) : 

0.82 

Groundspeed error (XD) ; 

0.5 

Integral of XD (XDl) ; 

0.1 




TABLE C-2. PLANT.- CONTROLLED ELEMENT AND SHAPING FILTERS 
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TABLE C-2. (Continued) 
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e- i-h q 
< m 

£ Q 

U 
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TABLE C-2. (Concluded) 


I 
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000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0. 000 1.00 1 XI 4 



TABLE C-3- PLANT,* OUTPUT AND MEASUREMENT EQUATIONS 


D > ^ 


% ^ 

-> % 

s , 


% 

*&> 


-b%, 



I 
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941E-08 -0 . 679E-09 -0.159E-07 0.000 0.000 -0.286E-16 0.000 0.195E-08 0.000 ! E13 



TABLE C-5. CONTROLLED ELEMENT TRANSFER FUNCTIONS 


CASE: UH1H HOVER 122LONG 31-JAN-79 CONTROLLED ELEMENT TF'S 


DENOMINATOR: 


1. 0000 

{ .00000 ) 

( .38494 ) 

(( .10339 
(( .26279 , 

< . 11863E-01 > 


( .00000 

.18934 

.92717 


) ( .00000 ) ( .00000 ) 

. 19576E-01 , .18833 )) 

.24365 , .89459 )) 


NUMERATOR: TH/DBD 


16910 



. 00000 ) . 

( .00000 ) 

( .00000 

.33300 ) 

( .39184 ) 


. 17445E-03> 




(- . 790 65E-0 2 ) 


NUMERATOR: TH/DCD 

33000E-02 

( .00000 ) ( .00000 ) ( .00000 ) ( . 108S4E-01 ) 

( .33300 ) (-11.270 ) 

< ,13442E-03> 


NUMERATOR: HDE/DBD 


. 24729 


( .00000 

) ( .00000 

(-1.4789 

) 

(( .35965 

, 1.3542 

<-.22334 

> 


) ( .00000 ) ( .33300 ) 

.48705 , 1.2636 )) 


HDE/DCD 


NUMERATOR: 

-9.7980 
( .00000 
(( .18545 
(( .25760 
<-.17544 


) ( .00000 
, .14356 
, .93210 


) ( .00000 
, . 26623E-01 , 
, .24011 , 


) 

.14107 )) 

.90064 )) 


NUMERATOR: 

-9.7980 
( .00000 
(( .18545 
(( .25760 
<-. 17544 


/DCD 


) ( .00000 
, .14356 
, .93210 


) 

, . 26623E-01 , 
, .24011 , 


.14107 )) 

.90064 )) 


NUMERATOR: XD/DBD 

1.0606 

( .00000 ) ( .00000 ) ( .00000 ) ( .33300 ) 

( .38296 ) 

(( .19712E-01, 2.2663 , .44674E-01, 2.2659 )) 

< .69473 > 


NUMERATOR: XDI/DBD 

1.0606 

( .00000 ) ( .00000 ) ( .33300 ) ( .38296 ) 

(( . 19712E-01 , 2.2663 , .44674E-01, 2.2659 )) 

< .69473 > 
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u) l rad /sec) 


K/s n+1 

( 2 n <x§ TH ) 1 Qdb 

APPROXIMATION 

KIM 

Asymptotic 
K = 0 . 1 69 ^ n = 1.0 

183 

G(jcd) 

K = 0 . 316 , n = 1.5 

2,831,632 


Figure C- 1 . Approximating Transfer Functions 
Selecting Value for Qo 














TABLE C-6. SUMMARY OF Q R AND R R SELECTION 
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Q MATRIX DIAGONAL, REGULATOR CLOSED LOOP EIGENVALUES, REGULATOR 


cm o* roX *r co mu vd cj r* i-* ooD n 'sj < 


c-i^ 



TABLE C-8. CONTROLLER COEFFICIENT MATRICES 


HI*) NH fON 

r-i «s> 

CO u cu 


N 

N 



PQ 
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<>> 
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II 

II 
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(M O 
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538. 0.U 9E+04 -29. 9 0.518E-01 -0.717 -0. 252 -26. 6 -1.89 1 DBD 



TABLE C-9 


CONTROLLER TRANSFER FUNCTIONS 


CASE: UHlH HOVER 122LONG 2-FEB-79 CONTROLLER TRANSFER FUNCTIONS 


DENOMINATOR: 


1. 0000 



( .31507 

) i 

( .33300 

< 1.6383 

> 



( 15.614 ) 


NUMERATOR: 

DBD/ HD 


21. 541 

( .11140 

) ( .31682 ) 

( .33300 

< .25317 

> 



NUMERATOR: 

DCD/ HD 


-42.506 

( .11589 

) ( .31500 ) 

( .33300 

<-. 51672 

> 



NUMERATOR: DBD/ Q 


NUMERATOR: DCD/ Q 


538.16 

( 15.694 ) 

(< .79748 , . 24 4 2(3 , .19475 , .14734 )) 

< 503.67 > 


1.3019 

(-. 28591E-01 ) ( 1.2724 ) ('53.279 ) 

< 2.5324 > 


NUMERATOR: DBD/ TH 

1485.9 

( .24018 ) ( .33300 ) < 15.988 

<1901.3 > 


NUMERATOR: DCD/ TH 

79. 797 

( .33300 ) (-.45581 ) ( 2.4905 ) 

<-30.165 > 

FLOATING UNDERFLOW PC=036035 


NUMERATOR: 

DBD/ DB 



NUMERATOR: 

DCD/ DB 

-29.879 




.51744E- 

■02 

< .18186 

) < 

.33300 

) ( 15.423 ) 

(-.11856 

) { .33300 ) (-2084.8 

<-27.917 

> 



< .41471 

> 





FLOATING UNDERFLOW PC=036035 

NUMERATOR: 

DBD/ DC 





. 517B1E- 

01 



NUMERATOR : 

DCD/ DC 

( .29977 
<-21. 332 

) ( 
> 

.33300 

) (-4127.1 ) 

-27. 008 
(-. 94888E 

-01) ( .33300 ) ( .33773 





< .28821 

> 

NUMERATOR: 

DBD/HDE 



NUMERATOR: 

DCD/HDE 

-. 71669 
( .31507 
<-1.1741 

> ( • 
> 

. 33300 

) ( 15.614 ) 

4.1570 
{ .31507 
< 6.8103 

) ( .33300 ) ( 15.614 

> 

NUMERATOR: 

DBD/HDI 



NUMERATOR: 

DCD/HDI 

-. 25178 
< .31507 
<-.41249 

) ( . 
> 

,33300 

) ( 15.614 ) 

2. 1563 
( .31507 
< 3.5326 

) ( .33300 ) ( 15.614 

> 

NUMERATOR: 

DBD/ XD 



NUMERATOR: 

DCD/ XD 

-26. 641 
( .27441 
<-37.964 

> ( • 
> 

,33300 

) ( 15.595 ) 

-.43622 
( .15039 
<-.39283 

) ( .33300 ) ( 17.982 

> 

NUMERATOR: 

DBD/XDI 



NUMERATOR: 

DCD/XDI 

-1.8892 
( .31537 
<-3. 0950 

) ( ■ 
> 

,33300 

) ( 15.614 ) 

-. 22059E- 
( .31507 
<-. 36139E 

01 

) ( .33300 ) ( 15.614 

-01> 
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Longitudinal Gust Transfer Functions 


C.*»SE : UH1H HOVER 122L0NG 31-JAN-7? ALL LOOPS CLOSED 

FLOATING UNDERFLOW PC=036035 


denominator: 


1.0000 

< .10175 

< . 50520 
<< .92313 
<( .67380 

118.7? 


< .31507 

( 15.61-1 
.70723 
2.77-11 


) ( .33300 ) < 

) < 25.716 ) ( 

.837-18 r .3-1882 

1.8672 » 2 . 0-498 


.33595 

25.726 


Transient Response to -10.0 ft /sec ug Step Input 


numerator: hd/ 

NEW FILE 

-.98613E-01 
( . 00000 > 

( .48737 ) 

<C .60809 
(( .*9639 

<-.71328E-02> 


FILE NAME? HDUG.2CL 


< .10125 

( 25.726 
.40520E-01 . 
2.7215 . 


) ( .32805 1 ( .33300 

1 ( 12.347 ) 

. 24640E-01 . .32167E-01)) 

1.8966 . 1.9518 )) 


h 

(ft/see) 


L 


numerator: 

NEW FILE 


FILE NAME? BBUG . 2CL 


.5 r- 

& o [_ 


3 . 4937 
. 10110 
3.7979 
( .94433 

< .39054 
-1 .3571 


( .33300 

( 25.724 
.29608 
.90773 


) ( .49451 


) ( 2.3-103 


.7 S370E-01 ) ) 
.39836 )) 


e 

(rad) 


.005 


numerator: 

NEW KILE 


FILE NAME? TMUU.2CL 


-. 19090E-02 
( .891 19C-01 ) 

( 23.721 ) 

(< .50319 i 

<( .92523 
- . 67934E-02. 


( .33300 

( 33.150 
.21630 
.09364 


). ( .33300 


. 12142 
.33146 


) ( 7.5257 



jumfrator: X 

NEW FILE 

. 10372E-01 
( , 00000 
( 6.8075 
(< .93505 
{ (-.12977 
. 599 72 E- 01 


FILE NAME? XDUfi.ZCL 


.91390 
1 .8327 


I) ( . 33300 

) ( 39.714 

.85445 r 
-.34779 » 


5 0 

(in.) 


-2 sec - 


NUMERATOR ! 


FILE NAME? DCUG . 2CL 


1 .59 46 
. i l i6 

.-ie 28 -i 

.?'.<076 
. 5/479 
1.1518 


.84304 

2.7655 


. 757.58 

1.8667 


) < .34074 


36614 

.0404 
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Figure C-2. Closed-Loop System Transfer Functions, 
Transient Response to Commands, Transient 
Response to Disturbances 
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Normal Gust Transfer Functions 


CASE: UH1H MOVER 122LONG 31-JAN-79 ALL LOOPS CLOSED 

iriNG UNDERFLOW PC^036035 


^nominator: 


1.0000 

( .10195 > 

( .31507 

) ( .33300 ) ( 

* 33595 

< .50520 ) 

( 15.614 

) ( 25.716 ) ( 

25 . 726 

(( .92313 r 

.90723 

f .83748 f .34882 

) ) 

(( .=7330 f 

2.7741 

f 1.8692 f 2.0498 

) ) 

' 119.79 > 

ATING UNDERFLOW 

FC=036035 



JMERATOR: HD/ UG FILE NAME? HDWG.2CL 

NEW FILE 


-.39580 

( .00000 ) 

( • 33300 > 

< <12.790 > 

<( .63102 
<-.96393E-01> 


( .15150E-01) < .10199 ) ( .32014 

( .33892 ) ( .51520 > ( 25.726 


2.8095 i 1.7722 


2.1787 


) > 


Transient Response to -1.0 ft/sec w g Step Input 


f, 

(ft/sec) 


.5 r- - — 

0 L. 


umerator: db/ uig file name? dbug.gcl 

NEW FILE 


-6.6001 
< .10211 
(( .99749 

(( .95584 

(( .90319 
1 2.5732 


( .33300 

* 33096 
.72377 
1.0757 


> (-6.7123 


.33013 

.69182 

.97155 


, 23482E-01 ) ) 
.21270 )) 
,46175 )) 


•iumerator: th/ wg file name? thug.ecl 

NEW FILE 


Sc 

(in.) 


'[ 


.005 

8 

(rad) 

O 



.33000E-02 
( .10290 ) 

(( .99403 

(( .92193 

(( .77974 


< .?1509E-01> 


( .33300 
.34557 
.91517 
23.639 


) ( .33300 


.34351 

.94369 

22.302 


1 < 25.726 

.37710E-01 ) ) 
,35460 )> 

17.966 >> 


X 

(ft/sec) 



NUMERATOR: t XD/ UG FILE NAME? XDUG.2CL 
NEW FILE 


. 21957E-02 

( .00000 ) 

( .33468 > 

< 91.909 > 

(( .90770 

< .34359 > 


< . 14809E-01 ) ( .27605 

(-7.0035 ) ( 25.726 


( .33300 

(-29.403 


.90729 


.92355 r .33073 >) 


NUMERATOR: DC/ UG FILE NAME? DCWG.2CL 

NEW FILE 



[-*- 2 sec -— ( 


13,007 
( .10188 
( .51001 
(( .91195 

(( .67067 

< 4.7364 


> ( .31464 

) ( 25.726 

r .84947 
r 2.7835 


) ( .33300 > ( .33590 

) 

,7746a f .34955 )> 

1.3668 f 2.0647 )) 


(concluded on following page) 


Figure C-2. (Continued) 
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Rale -of- Climb Command Transfer Functions 


CASE! UH1H HOUER 122L0NG 31-JAN-79 ALL LOOPS CLOSED 


denominator: 


1.0000 
t .10195 
( • 50520 

(( .92313 
(< *67300 
< 119.79 


< .31507 

< 15.614 
.90723 
2.7741 


) ( .33300 > ( 

> ( 25.716 > < 

. 03748 r .34882 

1.0692 r 2.0498 


.33595 
25.726 
) ) 

) ) 


) 

) 


Transient Response to 1.0 ft/sec lie Step Input 


NUMERATOR: HD/HDC FILE NAME? HDH0C.2CL 

NEW FILE 


40.900 
< .10195 
( .50240 

<( .67251 
1 118.79 


( .31507 
< .52197 
2.7694 


) ( .33300 ) ( 

) ( 15.614 ) < 

1.0625 » 2.0496 


.33601 
25.726 
) ) 



numerator: dd/hoc file name? dbhdc.2CL 

NEW FILE 



-.71669 
( .10135 

(-2.2043 
(( .95577 
<( .73732 

< 2.6625 


( .31506 
( 15.615 
. 49323 
2.0271 


) ( .33300 

> C 25.741 
.47290 r . 
1.4956 r 1 


) ( .39446 

) 

14017 )) 

.3603 )) 


.005 

e 

(rod) 

0 



numerator: TH/HI'C FILE NAME? TUH0C.2CL 

NEW FILE XDHOC . 2CL 


-.71 391E-03 
( .31506 > 

<( .58193 

(( . 097B9 

(< .99980 


< .33300 

.10137 
6.9191 
16.332 


> < .33300 
•58992E-01 , 
6.2125 r 
16.330 r 


> < .37296 

. 82440E— 01 ) ) 
3.0439 >> 

.23708 )) 


X 

(ft/sec) 


.05 r 
0 


numerator: xd/hpc file name? xdhdc.dcl 

NEW FILE 


.5 r 

S 9 

(in.) „ 


-.80053 
< .00000 
( ,33773 
<< .59773 
; 7.7291 


> ( .31503 

) <-.57722 

f 1.7890 


) ( .31968 > < .33300 ) 
> < 15.614 . ) < 25.741 ) 
1.0514 f 1.4474 )> 


K 2 sec —H 


numerator: 

NEW FILE 

DC/HDC FILE NAME? DCHUC.2CL 


4.1570 

< .10235 

) < .31507 ) ( .33300 > < 

.33684 

< .38233 

< 25.726 

) < .47213 > < .55313 > < 

> 

15.614 

<C .67330 
' 4.6420 

> 2.7743 * 1.3680 r 2.0513 

) ) 


Figure. C-2. (Concluded) 
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Phase (deg) ” Amplitude ( dB ) 



CASE: UH1H HOVER 


ALL LOOPS CLOSED 


denominator: 


1.0000 
l .10195 
( .50520 
{ ( .92313 
<( .67390 
< 118.79 


( .31507 
C 15.614 
.90723 
2.7741 


) ( .33300 ) < .33595 

) ( 25.716 ) < 25.726 

,83748 , .34092 >> 

1.8692 ? 2.0498 >> 


NUMERATOR: HD/HDC FILE NAME? HDHDC.2CL 

NEU FILE 


40.909 
' .10195 

' .50240 
< ( .67251 
< 118.79 


( .31507 
( .52197 
2.7694 


) ( .33300 ) < .33601 > 
) ' 15.614 ) i 25.726 ) 
1.8625 t 2.0496 >> 


j}_ ± ^ 1 . 0(8 + .^ 22 ) 

h 0 (s + 25.T)[s 2 +2(.923)(.91)s+ (-91) 2 ] 


Figure C-5. 


Closed- 


Loop Frequency Response: 


h/h c 
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case: UH1H HOVER 122L0NG 31-JAN-79 ALL LOOPS CLOSED 


DENOMINATORS 


1 .0000 

( .10195 > < .31507 ) ( .33300 > ( .33595 ) 

( .50520 ) ( 15.614 ) ( 25.716 ) ( 25.726 ) 

(( .92313 t .90723 . .33748 , .34882 )) 

(( .67380 r 2.7741 r 1.8692 r 2.0498 )) 

< 118.79 > 


NUMERATORS TH/VTH FILE NAME? THVTH.2CL 
OLD FILE 

-251.70 

< .00000 ) ( .51016E-03) < .24046 > C .33300 ) 

( .33300 ) < 15-968 ) < 25.716 > 

<< .92283 r .90699 t .83699 » .34939 >> 

<-1.1565 > 


-251.7s 2 

(s+ .l)(s+ .51)(s + 25.7)[s2 + 2 (.67)(2.77)s+ (2.77) 2 ] 


Figure C-5* Closed-Loop Frequency Response: e/v g = e/e c 




CASE 


UHIH l-IOVER 122LQNG 31-JAN-79 ALL LOOPS CLOSEP 


rEMOMrNAroR: 


1 .0000 

( .10195 
( .50520 
(( .92313 
<< . 37380 
■: 118.79 


< .31507 
( 15.614 
.90723 
2.7741 


1 < .33300 ) < 
) < 25.716 ) t 
.83748 . .34832 
1.8692 , 2,0498 


.33595 
25.726 
) ) 

> > 


numerator: xd/vxo 


-23.252 

( .ooooo > 

< 15.594 ) 

<< .92390 t 

<< . 20097E-01 r 

'.-1456.9 > 


< .27433 

( 25.716 
.90551 
2.2706 


) < .33296 ) < 

> 

.83660 r . 3464B 
. 45632E-01, 2.2702 


33300 

) ) 

) ) 


) 


numerator: xdixvxi 


-2.0036 

< .31507 ) < .33295 

< 25.716 > 

<( .92325 » .90595 

<< • 20122E-01 r 2.2693 
■1-118.79 > 


> < .33300 ) < 15.614 ) 

.83642 r .34806 )) 

. 45663E-01 t 2.2689 >) 


-28.3[s 2 + 2(.02)(2.27)s + (2.27) 2 ] 

(s + .5l)(s + 25.7)[s 2 +2(.67)(2.77)s + (2.77) 2 ] 


Figure C-6. 

(*/▼*) 


Closed- Loop Frequency Response; 
+ (/ *dt/vr ) = i/x c 




Phase (deg) "" "" Amplitude (dB) 



CASE! UH1H HOVER 122L0NG 31-JAN-79 ALL LOOPS CLOSED 


denominator: 

1.0000 


( 

. 10195 

) ( .31507 

) ( .33300 > ( 

♦33575 

< , 

.50520 

) < 15.614 

> ( 25.716 ) ' 

25 • 726 

< ( 

.92313 

, .90723 

» .83748 , .34832 

) ) 

< ( 

.67380 

, 2.7741 

7 1,8692 7 2.0499 

>> 


113.7? 


numerator: xd/vxd 


-2S.252 

( .00000 > ( .27433 

( 15.594 > < 25.716 

(< .92390 r ,90551 

!< .20097E-01, 2.2706 
<-1456.9 > 


} ( .33296 ) < .33300 J 

) 

.83660 i .34648 )) 

. 45632E-01 r 2.2702 )> 


numerator: xdi/vxi 

-2.0036 

( ,31507 ) ( .33295 

( 25.716 ) 

<( .92325 , .90595 

(( .201 22E -01 , 2.2693 

<-118.79 > 


) ( .33300 ) < 15.614 ) 

.83642 r • 34S06 )) 

■45663E-01, 2. 2639 >) 


fe ^ (s — . 001 ) ( s + .001) 
Xc (s + 1 ) ( s + • 51 ) 


Figure C-7* 
(x/ v x) + 


Closed-Loop Frequency Response 

(fxdt/vr + 1) = x e /£ c 


Phose l deg) Amplitude ( d8 ) 


&□ 


4-D 


2D 


□ 


D 


-XDO 


-2DD 


Figure C-8. 




</> u = 74 deg 


CASE: UH1H HOVER 122LONG 2-FEB-79 HDE AND HDI OPENED LOOP 


DENOMINATOR: 


1. 0080 
( .OOeOO 
( .33300 
( 25.726 
(( .67207 
< 17.547 


) ( .10199 

) ( .33575 

) { 27.232 

, 2. 7736 


) ( .11459 ) { . 31507 > 

) ( . 50556 ) ( 1 5. 614 ) 

) 

1.8663 , 2.0518 )) 


NUMERATOR: HDE/VHE 


-40.908 
( .00000 
( .33609 
(( .67265 
<-229. 16 


) ( .10195 

) ( .50555 

, 2.7692 


) ( .31507 ) ( .33300 

) ( 15. 614 ) ( 25. 726 

1.8627 , 2.0491 )) 


) 

) 


NUMERATOR: HDI/VHI 


-21. 190 
( .10195 

) ( .31507 

) < .33300 > ( 

.33506 

( . 50574 

) ( 15.614 

) ( 25.726 ) 


(( .67328 

, 2.7706 

, 1. 3654 , 2.0406 

)) 

<-118.79 

> 




) 


_h_ ^ -40.91 (s + .51) 

/ s( s + . 11 5 ) ( 3 + 27*37 

hc OL 


Open-Loop Frequency Response : h/h c ^ ^ dt Qpen 
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CASE: UH1H HOVEH 122LONG 2-FEB-79 HDE AND HDI OPENED LOOP 


DENOM INATOR: 
1.0000 


( .00000 
( .33300 
( 25. 726 
(( .67287 
< 17.547 

) ( .10199 

) ( .33575 

) ( 27.282 

, 2. 7736 

> 

) ( .11459 
) ( .50556 

) 

, 1. 8663 , 

) ( .31507 
) ( 15.614 

2.0518 )) 

NUMERATOR: HDI/VHI 



-21.190 
( . 10195 
( .50574 
(( .67328 
<-118.79 

) ( .31507 

) ( 15.614 

, 2.7706 

> 

) ( .33300 

) ( 25.726 

, 1. 8554 

) ( .33506 

) 

2. 0485 )) 

fh e dt 


~21 . 

2 

he 

s( 

OL 

>+ .51)(s + 25.7) 


Figure C-9* Open-Loop Frequency Response; 

I h e dt/h c . # 

/ h e dt open 





Phase (deg) Amplitude IdB) 



CASE: UH1H HOVER 122LONG 2-FEB-79 XD AND XD I OPENED LOOP 

DENOMINATOR: 

1 . 0000 


( .00000 ) 

( . 36768E-04 ) ( .23482 

) ( 

. 33283 

( .33300 ) 

( 15.611 ) ( 25.716 

) < 

26.961 

(( .9236S 

.90563 , .83649 

. 34707 

)) 

(( .50066 
< .84033E-01> 

3. 1452 , 1. 5999 

2.7079 

)) 


NUMERATOR: 

XD/VXD 



-28.252 

( , 00000 

) ( .27433 

) ( .33296 ) ( 

. 33300 

( 15.594 

) ( 25.716 

) 


(( .92390 

, .93551 , 

, .83650 , .34648 

) ) 

(< . 20097E 

<-1456.9 

-01, 2.2706 , 

> 

, . 45632E-01 , 2. 2702 

) 1 


NUMERATOR: 

XDI/VXI 



-2. 0036 




( .31507 
( 25.716 

} ( .33295 

) 

) ( .33300 ) < 

15.614 

(( .92325 

, .90595 , 

, .83642 , .34806 

) ) 

(( . 20122E 

-01, 2.2693 

r . 45563E-0 1 , 2. 2689 

)) 


<-118.79 > 


X 


OL 


-28.3[s 2 + 2(.Q2)(2.27)s+ (2.27) 2 ] 
s(s + 27.0)[s2 + 2(.51)(3.1)s+ (3.T) 2 ] 


Figure C-10. 


Open-Loop 

=A C | . r. 

|x, J X 


Frequency Response; 
dt open 
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uj( rad /sec) 


20dB/decade 

i * ,02 rad /sec 


<£ m - 83 deg 


CASE: UH1H HOVER 122L0JIG 

i 31-JAN-79 XDI OPENED 

LOOP 

DENOhlNATORJ 




1.0000 




< .00000 

) I .31307 

) C .33300 > 

.33412 ) 

( ,64743 

> < 15.614 

> < 25.714 > < 

25.723 ) 

(( . 92324 

, .90740 

f .83776 f .34864 

) > 

<< .67371 

r 2.7471 

r 1.8515 f 2.0293 

) ) 

< 1456.9 

> 

• 


numerator: 

XDI/VXI FILE 

NAME? XDIVXI 


NEW FILE 




-2 . 0036 




( .31507 

) < .33295 

) ( .33300 ) < 

15.614 ) 

( 25.716 

) 



(« .92325 

t « 90575 

f .33442 f .31906 

> > 

(( . 20 122E 

-01* 2.2693 

f . 45663E-0 1 f 2.2609 

> ) 

<-118.79 





-2.0[s 2 + 

2(.02)(2.27)s + 

(2.27) 2 ] 


s(s+ .65)(s + 25-7) [s 2 + 2(. 67) ( 2 . 75 ) S + (2-75) 2 ] 


Figure C-11. Open-Loop Frequency Response: 


f xd t/i 


f xdt open 






CASE: UH1H HOVER 12ZLONG S-FEB-79 XD,XDI AND TH LOOPS OPENED 


DENOMINATOR: 


1.0000 
( .00000 

) 

{ . 30483 

) ( .33300 

) ( .33863 

( 3.5187 

) 

( 15.568 

) ( 25.714 

) ( 26. 577 

{{ .77420 

, 

. 16484E-01, 

-12762E-01, 

. 10433E-01 ) ) 

(l .92994 


.90422 , 

.84086 , 

.33250 )) 

< . 28597 

> 





) 

) 


NUMERATOR: TH/VTH FILE NAME? THVTH 

NEW FILE 


-251. 70 
( .130000 
( .33300 
(( .92498 
< .20183 


) (-. 89352E-04 ) < .23999 ) { .33300 

) ( 15.961 ) ( 2S.721 ) 

, .90635 , .83836 , .34442 )) 


) 


_ Q _ ± , -251. 7 s 

0 C (s + 3.52)(s + 26.6)[s 2 + 2(.T7)(.oi6)s+ (.016)2] 


Figure C-12. Open-Loop Frequency Response: 
0 / 0c lx, f x dt, 0 open 
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TABLE C-10. CLOSED-LOOP SYSTEM RMS PERFORMANCE 
(PROCESS NOISE AUGMENTATION REMOVED) 
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APPENDIX D 


UH-1H, 100 KT, EXAMPLE APPLICATION 


This appendix provides an overview of the design synthesis and assess- 
ment of a longitudinal flight control system for the UH-1H at cruise 
(100 kt). No design iteration steps are included^ just the final design 
data are given. Only a small sample of the design assessment data is 
included. The next paragraphs explain the data which follow. 

Table D-1 lists the primary design goals for this application. These 
design goals are the various bandwidths, which are determined from the 
functions intended to be accomplished by this FCS. 

Tables D-2 and D-3 present the input data used to define the system. 

The shaping filter and controlled element state equations are given in 
Table D-2, and the measurement and output equations follow in Table D-3. 

Table D-U lists the elements of the final filter -observer design. 

Table D-5 lists the various controlled element transfer functions of 
interest, while Table D-6 presents the procedure used to meet the band- 
width requirements by choosing the elements of % and R R for the regulator. 
As in the hover example, we were able to make use of the relatively un- 
coupled nature of this problem to design all of the loop closures simul- 
taneously. Table D-7 provides a summary of the final regulator design. 

With the filter-observer and the regulator synthesized, the software 
is used to obtain the controller gain matrices and the controller transfer 
functions, given in Tables D-8 and D-9, respectively. 

Design assessment begins with the closed-loop transfer functions, 
shown in Fig. D-1, and the corresponding transient responses to Ug, Wg 
and £l c step inputs, also shown in Fig. D-1. These transient response 
plots enable us to determine the behavior of the closed-loop system in 
the presence of gust disturbances, as well as the behavior in response 
to command inputs. 


D-1 



Subsequent figures provide data on the frequency response character- 
istics of the closed-loop system. Figure D-2 presents the rate -of -climb 
error response to rate -of -climb command. Figure D-3 plots the pitch 
attitude response to pitch attitude command. The frequency response of 
longitudinal velocity to a longitudinal gust disturbance input is plotted 
in Fig. D-4, while the Bode plot of airspeed error to longitudinal gust 
disturbance is shown in Fig. D-5* These are only a few of the many inter- 
esting frequency response plots for this design. 

An example of the statistical data available in the design assessment 
stage is shown in Table D-10. Here the process noise augmentation has been 
removed and we can examine the rms response of the plant and filter states 
to the combined process and measurement noise. We can also look at the 
rms response of the outputs and controls. These rms values are used to 
determine control activity levels, including positive -going zero crossing 
rate, which is shown to the right of the rms data. The control activity 
data are used in determining whether appropriate cost function weightings 
were chosen for the control rates (the matrix) in the regulator design. 
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TABLE D-l 


EXAMPLE APPLICATION 


UH-1H Longitudinal Flight Control System 
Flight Condition: 100 kt (Cruise) 

Functions : 

Rate-of-climb command 
Airspeed hold 

Band-widths: (rad/sec) 

Cyclic (DB): 25-73 

Collective (DC): 25-73 

Pitch (TH): 1-5 

Rate-of-climb error (HDE) : 0-7 

Integral HDE (HDl): 0-32 

Airspeed Error (ASE): 0.1 

Integral ASE (ASl) : 0.02 
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TABLE D-2. PLA NT; CONTROLLED ELEMENT AND SHAPING FILTERS 


CO ZZ <J\ ^ CM 


V-i, H 

i*9 . ^ 0 

% Q % 

°-s> 

^ ^ PQ 

7s. ® 0 « 

a 

% H 

•^b/V § 

n 




V g 


V^O ^ ^ 








\\ ^ 

^ % 

■&> ^ .O 


(continued on following page) 


TABLE D-2. (Continued) 


<D 

no 

nJ 

P< 


a* 

•rl 

§ 

t — i 

H 

O 

£ 

0 

CD 

1 
O 
£ 
o 
o 



8 

Q 


Q 

PQ 


G Q 

CNi O IS. G 

(J G ® 

Q 

G G 


Ci G G 
G Q G 
C3 Q Cx 

G G G 


G 

G 

Ct 

G 


<3 

G 

Ci 

G 


Q 

G 

G 


G 


G C 

G G 

G G 

G G 



G G 

G G 

G G 


ea G 

G> G 

G G 


G 


G 


G 

Q 

G 


I 


G 

G 

G 


G B G G 
G G G G 
G .GO 


G G 

Q G 

G G 


G G 


G C 


D G 
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TABLE D-2. (Concluded) 




TABLE D-3- PLANT,- OUTPUT AND MEASUREMENT EQUATIONS 
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TABLE D-4. FILTER -OBSERVER SYNTHESIS 



m • —* • r-* . 
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TABLE D-5. CONTROLLED ELEMENT TRANSEER FUNCTIONS 

CASE: UH1H 100KT 128LONG 30-JAN-79 CONTROLLED ELEMENT TF ' S 


DENOMINATOR: 

1. 0000 

( .00000 ) ( .00000 } ( .00000 ) ( .00000 ) 

( .69142 ) 

(( .22354 , .17435 , .38973E-01, .16993 )) 

(( .43625 , 1.4971 , .65312 , 1.3472 )) 

< .4710SE-01> 


NUMERATOR: TH/DBD 

-.17280 

( .00000 ) ( .00000 ) ( .00000 ) ( . 14710E-01) 

( .33300 ) ( 1.2044 ) 

<-. 10194E-02> 


NUMERATOR: TH/DC D 

22300E-01 

( .00000 ) ( .00000 ) ( .00000 ) ( . 29081E-01 ) 

( .33300 ) (-3.1286 ) 

< . G7562E-03> 


NUMERATOR: HDE/DBD 

5. 0374 

( .00000 ) ( .00000 ) ( .00000 ) ( .42303E-01) 

( .33300 ) 

(( .15642 , 2.6370 , .41248 f 2.6046 )) 

< .49346 > 


NUMERATOR: H DE/DC D 

-13.169 

( .00000 ) ( .00000 ) ( .00000 ) 

(( .70493 , .17927 , .12637 , .12715 )) 

(( .30308 , 1.3401 , .40618 , 1.2771 )) 

<-.76012 > 


NUMERATOR: HDI/DCD 

-1 3. 1 G 9 

( .00000 ) ( .00000 

(( . 70493 , .17927 

{( .30303 , 1.3401 

<-.76012 > 


) 

, . 1 2G37 , .12715 )) 

, .40618 , 1.2771 )) 


NUMERATOR: ASE/DBD 


. 82341 

( .00000 ) ( .00000 

( 1.2596 ) 

(( . 73104E-01 , 2.5391 
< 2.2266 > 


) ( .00000 ) ( .33300 ) 

.18562 , 2.5323 )) 


NUMERATOR: ASI/DBD 

. 82341 

( .00000 - ) ( .00000 ) ( .33300 ) ( 1.2596 ) 

(( .73104E-01, 2.5391 , .18562 , 2.5323 )) 

< 2.2265 > 
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TABLE D-6. SUMMARY OF Q R AMD R R SELECTION 
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CLOSED LOOP EIGENVALUES, REGULATOR 


D-1 1 



TABLE D-8. CONTROLLER COEFFICIENT MATRICES 


AF MATRIX 


1 

2 

3 


EI3 

E01 

E02 




i 

1 

! -0. 333 

-0. 614E-08 

0. 130E-06 J 

El 3 



i 

2 

! 8. 06 

-0.183 

1. 52 ! 

E01 



r 

3 

j -125. 

1.17 

-25.0 ! 

i 

E02 

*_ 





BF MATRIX 
1 

2 

3 

4 

5 

6 

7 

8 

9 

HD 

Q 

TH 

DB 

DC 

HDI 

HDE 

AS I 

ASE 

-0. 314E-05 

4.97 

0. 530E-03 

0. 657E-06 

-0. 171E-05 

0. 263E-12 

-0. 180E-12 

-0. 908E-11 

0. 162E 

-36. 3 

261. 

0. 613E+04 

a) 

Si 

-20. 1 

-0. 145E-04 

-0. 157E-03 

0. 478E-03 

1.87 

594 . 

-0. 419E+04 

-0 . 100E + 06 

-125. 

325. 

-0. 422E-04 

0. 546E-04 

0. 147E-02 

-30. 6 


CF MATRIX 


1 

2 

3 


E 13 

E01 

E02 

1 1 

-23. 9 

-0. 558 

-1 . 22 

i DBD 




i 

J 2 

0.934 

-0. 313 

0. 479 

! DCD 


I 


DF MATRIX 


l 

2 

3 

4 

5 

6 

7 

8 

9 

HD 

Q 

TH 

DB 

DC 

HDI 

HDE 

AS I 

ASE 

28.1 

271. 

-0. 434E+04 

-27.9 

0.627 

-1.34 

-2.99 

-0. 193E-01 

-2.15 

-11.6 

-2. 14 

0 • 20 6E +0 4 

0. 627E-01 

-25. 9 

0. 324 

0.981 

-0. 786E-02 

0.281 


1 

E 1 3 
2 

E0 1 
3 

E02 


1 

DBD 

2 

DCD 
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TABLE D-9. CONTROLLER TRANSFER FUNCTIONS 


CASE: UH1H 100KT 128LONG 2-FE8-79 CONTROLLER TRANSFER FUNCTIONS 


DENOMINATOR: 

I . 0000 

{ .11136 ) ( .33300 ) ( 25.074 ) 

< .92977 > 


NUMERATOR: DBD/ HD NUMERATOR: DCD/ HD 


28. 057 



-11. 590 




(-.817G1E-01) 
<-. 990 96E-01> 

( .12973 ) 

( .33300 } 

( .10443 
< .19091 

) ' 
> 

[ .33300 

) (-.47364 


NUMERATOR: DBD/ Q 


NUMERATOR: 

DCD/ Q 


270.94 

( .57922E-01) ( .19074 

< 129.19 > 

) ( 43.160 ) 

-2. 1404 
( .14888 
<-128.03 

) ( .40243 

> 

1 ( 998.36 


NUMERATOR: 

DOD/ TH 


NUMERATOR: 

DCD/ TH- 


-4344.4 
( .11516 
< 386.40 

) ( .33300 

> 

) (-2.3193 ) 

2059.9 
( .12469 
< 68.546 

) ( .33300 ) 

> 

( .80143 


NUMERATOR: 

DBD/ DB 


NUMERATOR: DCD/ DB 


-27.876 
( .12362 
<-22. 684 

} ( .33300 ] 

> 

1 ( 19.767 ) 

. 62709E-01 

( .16364 ) ( .33300 

<-3.3050 > 

) (-967.21 


NUMERATOR: 

DBD/ DC 


NUMERATOR : 

DCD/ DC 

. 42710 
( .10160 
<-12. 500 

) { .33300 

> 

) (-589.14 ) 

—25. 870 
( .10784 
<-17. 484 

) ( 
> 


NUMERATOR: DBD/1IDE 


NUMERATOR : 

DC D/HDE 


-2. 9S85 

( .11133 ) { .33300 ] 

<-2.7778 > 

1 ( 25.074 ) 

. 98066 
( .11140 
< .91216 

) ( .33300 : 

> 

1 ( 25.074 


NUMERATOR : 

DBD/HDI 


NUMERATOR: 

DCD/HDI 


-1.3402 
( .11135 
<-1.2460 

) ( .33300 ] 

> 

> ( 25. 074 ) 

. 32357 
( .11137 
< .30089 

) ( .33300 1 

> 

1 ( 25.074 


NUMERATOR: 

DBD/ASE 


NUMERATOR: 

DC D/AS E 


-2.1499 
( .11284 
<-.66249 

) ( .33300 ] 

> 

) ( 8.2005 ) 

. 281 38 
( .11284 
<-.30735 

) ( .33300 

> 

) (-29. 069 


NUMERATOR: D8D/ASI 

-. 19271E-01 

( .12935 ) ( .33300 ) ( 25.162 ) 

<-. 20387E-01> 


NUMERATOR: DCD/ASI 

-. 78565E-32 

( .13233 ) ( .33300 ) ( 24.982 ) 

<-. 864 90S -0 2 > 
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Transient Response to -1.0 ft/sec Wn Step Input 
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(continued on following page) 
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Figure D -1 . (Concluded.) 


Phase (deql Amplitude IdB) 


u ( rad /sec ) 




• 'jq w ( rod/sec) ^ 


CASEs UH1H 100KT 12SL0NG 31-JAN-79 ALL LOOPS CLOSED 
DENOMINATOR: 

I. 0000 

( . I5598E-0 1 ) ( . 626S9E-01 ) ( .11136 ) ( .33300 

( .33514 ) { 25.074 ) ( 25.736 ) ( 25.736 

(( .99458 , .72781 , .72387 , .75647E-01)) 

(( .56441 , 2.2044 , 1.2442 r 1.8197 )) 

< .51948 > 


NUMERATOR: TH/VTH FILE NAME? THVTH.CL 

NEW FILE 

704.78 

( .00000 ) ( . 1181 6E-01 ) ( .12216 ) ( .21162 

( .33300 ) ( .33300 ) ( 1.2135 ) (-2.7713 

( 25.736 ) 

<-2.0661 > 


705- Os (s + .212) (s + 1.2) (s -2.77 


(s+ .063) (s + 25- 1 )(s + 25 . 7)[s^ + 2( -99)( . 75)s + ( . 73) 2 ][s 2 + 2 ( .56)(2.2)s + (2.2) 2 ] 


Figure D-5* Closed-Loop Frequency Response; e/v 9 = 0/0 c 
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CASE: OH 1H 100KT 128LONG 31-JAN-79 ALL LOOPS CLOSED 


DENOMINATOR: 


1.0000 

( . 1 5598E-0 1 ) 

( .33514 ) 

(( .99458 , 

(( .56441 
< .51948 > 


( . 62 689E-01 
( 25. 074 
. 7 2 7 1; 1 
2.2044 , 


) ( .11135 

) ( 25.736 

.72387 , 

1.2442 


) ( .33300 

) ( 25.736 

. 7564 7E-0 1 ) ) 
1.8197 )) 


) 

) 


NUMERATOR: X06/ UG FILE NAME? UUG.CL 

NEW FILE 


.451 00E-01 
( . 11484E-01 ) 
( .48113 ) 

( 46.923 ) 

(( .54154 
< .51818 > 


( .11322 
( 1.2669 

2. 1330 


( .33297 
( 23.801 


( .33301 
( 25.732 


1. 1551 


1.7931 


)) 


u ^ .04^(s+ .48)(s + 1 .3) 

”g (s + .063)[s 2 +2(.99)(-73 )s + (.73) 2 ] 


Figure D-4. Closed- Loop Frequency Response: u/u g 


D-19 



Phose (deg) Amplitude ( dB ) 


,□□1 ■ i ;□! • u ( rod /sec ) . 



^200 _ 


CASE: UH1H 1 0 i3KT 128LONG 31-JAN-79 ALL LOOPS CLOSED 

DENOMINATOR: 

1.0000 


( . 15598E-01 ) 

( . 62 6 89E -0 1 

) ( .11136 

) ( .33300 

( .33514 ) 

( 25.074 

) ( 25.736 

) ( 25.736 

(( .99458 

.72781 

.72387 , 

. 75647E-01 ) ) 

(( .56441 
< .51948 > 

2.2044 

1.2442 , 

1.8197 )) 


NUMERATOR: 

ASE/ UG FILE 

NAME? ASEUG.CL 


NEW FILE 




-.99860 




( .00000 

) (-. 64568E- 

-02) ( .11254 ) ( 

.33300 

( .33490 

) ( .65545 

) ( .79111 ) ( 

24.383 

( 25.736 

) ( 26.458 

) 


(( .55909 

, 2.2381 

, 1.2513 , 1.8556 

)) 


< 3.4897 > 

UAS e ^ -1 .Qs(s + .6^5) (a + .79) 

u g (s + .02)[s2 + 2(.99)(.T3)s + (.T3) 2 ] 

Figure D-5. Closed-Loop Frequency Response: 


TABLE D-10. CLOSED -LOOP SYSTEM RMS PERFORMANCE 
(PROCESS NOISE AUGMENTATION REMOVED) 
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APPENDIX E 


METHOD EXTENSION WHEN (^QT-g) IS SINGULAR 


APPROACH 


Augment the rg and. Eg matrices so that the rank of T2 is n W2 


r = 



and the first nw-j columns of r represent the process noise effects which 
are actually part of the first principles plant model. The remaining 
U-w — n w-| = nw2 components of process noise are added for convenience in 
obtaining alternative limiting forms for the singular filter. The inten- 
sities of these n^g components are allowed to approach zero according to 
some externally specified ratio relationship chosen by the designer. (This 
is akin to the arbitrariness in selection of observer poles by the designer.) 

This means that the diagonal Q matrix has the form 

q = diagjq!, q^ , \i, dqnw-,+2* ( 


where ^ = 1 and jjl approaches zero from above. Notice that there is 

also some arbitrariness in the elements of the n W g columns used to augment 
the Tg and Tq matrices. The effective maximum number of parameters which 
may be arbitrarily selected is (nnwg — 1 ) . The only restriction on para- 
meter choices is that (r2QT2) have full rank, r = mg, for q > 0 . Notice 
also that the columns augmenting r-| may be identically zero by choice. 

This in turn results in A-] 1 and (see below) being null matrices. 

The terms involved in the Euler-Lagrange equations for the optimal 

I *1 

filter are r-iQTi, r-|QP2 and • 



Separate the terms riQTi, r-|Qr2 (^ 2 ^ 2 ) ^ into components which 
are independent of and dependent on p. That is, let 

p 1 Qp -| = H-] 1 + M-A-| 1 

r 1^2 = H 12 + MA-,2 = (r 2 qr^) 

r 2 ^2 ~ ^22 + ^22 ~ A 


Next, consider the behavior of various factors occurring in the Euler- 
Lagrange equations in the limit as p approaches zero from above. 


lim 
P -i-o 


r lQTi’ 


1 


S (r 2^ r ’ 


(A (l > + A^/l-O 


P -► 00 


r i«'2(r 2 <3’a) 1 


lim 

e-*-0 


|h 12 A (1 ) + A 12 a(2) + E 12 A (2) /^ 


lim 


+T 0 r ]Qr 2 (r 2 Qr 2 ) 1 r 2 on' 


lim 

+ 


(D W 


p^O 1*12^ H 12 + A ^ 


(2) ^12 + S 12 A ^ A 12 + -12 a ^ 2 ^ S J 2 /^ 


The procedure continues by substituting the separated quantities into 
the Euler-Lagrange equation: 


E-2 


(si-F n ) r-|Qri 
H^RT’Hn (sI + Fn) 



0 

Hi'iRT 1 


f 12 

-HiiRiVig 



r 1 Qr2(r 2 QT2) 1 f 21 
F 21^ F 2^2^ F 21 


r^^^QTo) 1 g 2 

F2i(r 2 QT2) 1 G 2 


-riQT2(r2QT2) 1 ^qt] 
-F^CrgQT^) -1 r 2 Qri 



0 -riCT2(r2Qr 2 ) 1 (si-f 22 ) 

O -F2i(r2QT2) 1 (sI-F 2 2) 


r 1 qr2(r 2 Qr2) 

F2i(r2QT2) 1 


The result is; 

(si — F-| 1 + [E-| 2 A^ ^ 

+ A 12 A (2) + Hi2A^ 2 Vm- 1 f 21 ) 

(Ri'iRT’hh 

+ F 2 i [A^ 1 ) + A^ 2 V^1 f 21 ) 


(^11 ~^ 2 k ^' > ~\2 + A 12 A ^~12 
+ H 12 a ^ a 12 +H 12 A^ 2 ^Hj 2 /^]) 


(si + F-Ji — F 21 [A^ 1 

+ A^ 2 ^Ai2 + A^Hi2/V] ) 



11 

- 

(Hi2A^ 1 ‘ l +A 12 A^ 2 ^ +Hi2A l ' 2 V^) 

(x 10 ) 

1 

CM 

+ 

-CM 

h 

1 


(G-l -[H 12 A (1) 

0 

+ A 12 A ^ ^ + Hi2A^ ^ / M-3 Cx2 


( f 12 + [H^2 A *' 1 ^ + A 12 A^ 2 ) 

+ H^ 2 a ^ ^ /m-] ( S I — F 22^ ) 


+ 


-1 


H n Ri 


(-HnR/H^ + FgiU^ 1 ) 

+ A^ 2 ^/h][sI — F 22] ) 


E-3 


-Fg^A^ +A^ 2 V^]G 2 



The revised steps in solution are as follows. 

• Evaluate each element of each matrix in terms of 
numerical coefficient values for s^ s^ s^/p and 
s yp in each matrix. 

• Examine the resulting matrix equations via inspec- 
tion of the numerical coefficients. If any equation 
has one or more non-zero coefficients for s°/m. and/or 
s/p, then multiply that equation through by p and let 
u. -^-0. (This is equivalent to simply retaining only 
the coefficients for s°/p and sVp in that equation.) 

• Eliminate a number of variables corresponding to the 
number of equations which are rendered algebraic in 
xi and >1 as the result of the previous step, via 
substitution , 

• Apply eigenvalue decomposition (or matrix Riccati 
equation solution) technique to the resulting reduced 
Euler -Lagrange equations (refer to pages A-5 to A-7 
of Appendix A) . Notice that the form of the filter- 
observer solution must be modified to allow for dif- 
ferentiation of the measurements. Multiple differ- 
entiation may be required. 

• Use the algebraic equations identified in the next -to - 
last step to obtain estimates of the remaining states. 



APPENDIX F 


EFIECT OF AUGMENTING PROCESS NOISE INTENSITIES ON 
FILTER EIGENVALUES AND RMS ESTIMATION ERROR 


Data for three filter design passes is summarized in this appendix. 
The filter is for the longitudinal dynamics of the UH-1H helicopter in 
100 kt cruise. All measurements are noise free. The nine measurements 
include airspeed error and integral airspeed error as appropriate for 
cruise. The third and final filter design has acceptable rms estimation 
error and closed-loop eigenvalues which are within the acceptable range. 


F-l 



TABLE F-1. INITIAL FILTER DESIGN, 100 KT 



TABLE F-2. SECOND FILTER DESIGN, 100 KT 

0 MATRIX DIAGONAL, FILTER fWS STATE EST ERROR, FILTER 



TABLE F-3. FINAL FILTER DESIGN, 100 KT 




CC 

CO 



CO 

E-» 

in 


to 

E- 

< 

£- 

cn — < 

cn 

s: 

c: 


,-H rH cn CN m ro 

G G ca 

XXX 


CL CN 
H H 

X 


»-i m cm 


\ 


ro 

G 

CN 


CN 

O 


C2 

CL 

C£ 

G 


m 

G 

CN 


r- 

G 
*— i 

Q 


G 

G 

G2 

G 


CN 

r-i 

I 

CO 




G 

G 

G 


Gi 


G 

G 

G 




O 

CL 

G 


O 


V Z) 


a 


m 2 vdO r- — 

CL Oj Cl 


CT* CJ G CO -H 03 

q h r «“i x 

a a cl. 


cn to 
»-» cn 
o- 


N* 

CN 


N* 


CN 

C\ 


CN 


G 


I 




f-4 


L 



1. Report No. 

NASA CR-3275 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 

PRACTICAL OPTIMAL FLIGHT CONTROL SYSTEM DESIGN FOR HELICOPTER 
AIRCRAFT. Volume I - Technical Report 


5. Report Date 
May 1980 


6. Performing Organization Code 


7. Author(s) 

L. G. Hofmann, Susan A. Riedel, and Duane McRuer 


8. Performing Organization Report No. 
TR-1127-1-I 


10. Work Unit No. 


9. Performing Organization Name and Address 

Systems Technology, Inc. 

13766 South Hawthorne Boulevard 
Hawthorne, California 90250 


11. Contract or Grant No. 

NAS2-9946 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 


13. Type of Report and Period Covered 
Contractor Report 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Ames Technical Monitor: Heinz Erzberger 


16. Abstract 


1 


A method by which modem and classical control theory techniques may be integrated in a synergistic] 
fashion and used in the design of practical flight control systems is presented here. A general 
procedure is developed, and several illustrative examples are included. Emphasis is placed not 
only on the synthesis of the design, but on the assessment of the results as well. 

The first step here is to establish the differences, distinguishing characteristics and connec- 
tions between the modem and classical control theory approaches. Ultimately, this uncovers a 
relationship between bandwidth goals familiar in classical control and cost function weights 
in the equivalent optimal system. In order to obtain a practical optimal solution, it is also 
necessary to formulate the problem very carefully, and each choice of state, measurement and 
output variable must be judiciously considered. These so-called "engineering art" matters 
allow us to bridge the gap between the optimal control theory and its practical application. 

Once design goals have been established and problem formulation completed, the control system 
is synthesized in a straightforward manner. Three steps are involved: filter-observer 

solution, regulator solution and the combination of those two into the controller. Assessment 
of the controller, which is obten the bulk of the task at hand, permits, an examination and 
expansion of the synthesis results. Often, the composite picture which results may lead to a 
revised design which is simpler and more practical . 

Key contributions in this work include the solution of the singular Kalman filter problem, 
presented in an appendix, and the development of a user-oriented computer software package for 
the design and assessment of optimal flight control systems. Volume II of this report contains 
the software details in the form of a user's guide. Volume II is available as a high contractor 
report number 152306 on microfiche. 


17. Key Words (Suggested by Author (s) ) 

Optimal Control -Theory, Helicopter Flight 
Control Systems, Singular Kalman Filter, 
Optimal Control -Software, Optimal Control- 
Application, Cost Function Weights, Classical 
Control Analysis 


19. Security Classif. (of this report) 
Unclassified 


18. Distribution Statement 

Unclassified - Unlimited 


STAR CATEGORY - 08 


20. Security Classif. (of this page) 
Unclassified 


21. No. of Pages 
271 


22. Price 

$10.75 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1980 




